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1. Introduction 

Departure from thermal and chemical equilibrium occurs frequently in hypersonic flow. 
When a fluid element travels through a flowfield it will stay in equilibrium with its surroundings 
if the molecular collision rate is sufficiently high. On the other hand, if the fluid element 
experiences sudden changes in conditions and the collision rate is relatively low, then departure 
from thermodynamic equilibrium, or nonequilibrium, occurs. In hypersonic flight 
nonequilibrium can occur, for example, as gas passes through shock waves, suddenly expands, or 
undergoes combustion. The resulting nonequilibrium can affect the population of molecules in 
translational, rotational, and vibrational internal energy modes as well as the chemical 
composition. Such nonequilibrium phenomena, if ignored, can lead to erroneous predictions of 
lift, drag, aeroheating, and engine performance. 

In a gas in thermal equilibrium the distribution of population among its various thermal 
energy levels (translational, rotational, and vibrational) is governed by the Boltzmann 
distribution law, which is expressed in terms of a temperature. If the gas is not in thermal 
equilibrium, then no unique temperature can be assigned. In moderate density (non-rarified) 
hypersonic flows it is often useful to define separate vibrational and translational-rotational 
temperatures (where translation and rotational modes are assumed to be in equilibrium among 
themselves and with each other). For example, the translational-rotational temperature will rise 
abruptly across a shock wave, but the distribution of energy in the vibrational levels will take 
some time (distance) before it comes into equilibrium. It has been shown that under certain 
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conditions an intermediate vibrational temperature may be defined (i.e., a Boltzmann distribution 
of population among the vibrational levels) at various times in the relaxation. 1 In rarefied 
hypersonic flows, the rotational and translation temperatures can be unequal; furthermore they 
can have non-Boltzmann distributions.” 

Computation of fluid dynamical flowfields is greatly complicated when nonequilibrium must 
be considered. Instead of computing chemical equilibrium, finite rate chemistry must be 
included, resulting in the calculation of tens to thousands of equations instead of just a few. 
Similarly, thermal nonequilibrium conditions require bookkeeping of different temperatures for 
different energy modes such as translational, rotational and vibrational temperatures, T trans , T rot 
and Tvw, respectively, assuming that each of those energy modes can be described by a 
Boltzmann distribution. Further complicating matters, different species in the same flow can 
have different vibrational temperatures, or may even have nonequilibrium distributions that 
cannot be described by a temperature. For example, the products of some exothermic chemical 
reactions are formed with a non-Boltzmann distribution of vibrational energy. 1 

While thermodynamic and chemical nonequilibrium occur in many different fields of 
chemistry, physics and engineering, hypersonic flows are particularly susceptible to these 
complicating effects. Hypersonic vehicles (for example, entry vehicles or hypersonic cruise 
vehicles) typically fly at high altitude where the pressure is low. Because of the high speed, 
kinetic energy from the vehicle is converted to thermal energy in the gas flowing over the 
vehicle. Thus hypersonic gas flow typically encounters low pressure, P, and high temperature, 
T, conditions. Since gaseous collision rates scale as P/VT, low-pressure high-temperature flows 
can have relatively low collision rates. These low collision rates, coupled with flow phenomena 
like shock waves, expansions and boundary layers which have steep gradients, cause sudden 
changes in flow conditions, resulting in nonequilibrium. 

The different thermal and chemical phenomena mentioned above require different numbers 
of collisions to maintain equilibrium. Translational and rotational modes of a molecule can be 
brought into equilibrium with 10 or fewer collisions, occurring in less than 1 ns at P = 1 atm and 
T = 273 K. For comparison, it takes about 20,000 collisions to vibrationally excite an CF 
molecule. 4 Chemical changes can occur even more slowly. Dissociation of the CF molecule by 
collisions requires on the order of 200,000 collisions. 4 Thus vibrational and chemical changes in 
a gas occur much more slowly than translational and rotational changes. Consequently, 
vibrational and chemical nonequilibrium are more commonly observed in hypersonic flow than 
translational and rotational nonequilibrium, and measurements of these phenomena will be 
emphasized in this chapter. The complexity of modeling these nonequilibrium phenomena, and 
the strong effect that they can have on the heating of aerospace vehicles, makes it essential to 
develop experimental techniques capable of quantifying the degree of nonequilibrium in the 
flowfield without interfering with the distributions that they are measuring. The non-intrusive 
and quantum-state-specific characteristics of spectroscopic techniques make them uniquely 
suitable for this purpose. Spectroscopy, then, can provide direct experimental evidence of these 
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equilibrium or nonequilibrium thermal distributions and chemical compositions. Nonintrusive 
spectroscopic techniques can be used in ground-based flow simulation facilities and sometimes 
even in flight to probe individual atomic and molecular states, thereby measuring population 
distributions, translational, rotational and vibrational temperatures, chemical compositions, and 
other data that can help modelers compute these challenging flowfields. 

A variety of different molecular-based optical measurement techniques will be detailed in 
this chapter. Optical emission spectroscopy (OES) is the only technique detailed herein that does 
not require a light source (e.g. a laser). OES instead disperses naturally emitted light to provide 
quantitative flow field information. However OES can only be used to study luminous 
environments. Also it is difficult to determine spatially-specific information using OES since the 
collection optics average all the light along its optical path. These limitations can be overcome 
by using a light source such as a lamp or a laser. If a light source is incident upon a gaseous 
medium the incident photons can either be absorbed, scattered, or can transmit through the 
medium. These different possible light paths present opportunities for different methods to probe 
and measure the gas. Absorption spectroscopy uses the transmitted light to diagnose the sample: 
the method compares the transmitted intensity with a sample present to a reference thereby 
determining the sample’s absorption integrated along the optical path. Rayleigh scattering, 
Raman scattering and laser induced fluorescence (LIF) collect scattered light at an angle to the 
incoming light. This localizes the measurement point at the intersection between the incoming 
light and the optical path of the detection optics. Measurements can be made at a single spatial 
location or, using a camera, measurements can be made along a line or in an image plane. 
Whereas Rayleigh scattering is an elastic scattering process, Raman scattering is inelastic, with a 
small amount of energy either gained from or lost to the surrounding medium. LIF is a two-step 
process wherein a photon from the incident beam is absorbed and then is subsequently emitted. 
These scattering processes emit light approximately isotropically in three dimensions (into 4n 
steradians) so that the intensity of the light drops off as one over the distance from the 
measurement point squared ( inverse square law). Consequently, using conventional lenses, only 
a small fraction (typically ~1%) of the scattered light can be collected. 5 The abovementioned 
scattering techniques do not require the coherence property of a laser and are known as 
incoherent measurement techniques. 

Coherent measurement approaches involve the use of multiple laser beams which are 
directed into and interact with the flowfield. These techniques make use of the coherence of 
laser light to perform nonlinear wave-mixing which generates coherent signal beams. The 
resulting laser-like beam does not follow the inverse square law. Instead, it can propagate long 
distances without loss. It can also be spatially filtered to reject flow luminosity (which does 
follow the inverse square law) without losing signal intensity. For coherent measurement 
techniques, the measurement volume is determined by the location of the crossing of the laser 
beams. Unfortunately, these coherent techniques are much more complicated to set up and 
harder to keep aligned than incoherent techniques. Also, these methods are nonlinear and are 
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therefore generally more difficult to analyze and interpret than incoherent techniques. However, 
coherent techniques, such as coherent anti-Stokes Raman spectroscopy (CARS), provide high 
signal-to-noise ratio data in applications with limited optical access and/or flow luminosity. 

This chapter mainly focuses on spectroscopic measurement techniques for studying gaseous 
hypersonic nonequilibrium flow. In particular, we will focus on techniques that have either been 
used to study hypersonic nonequilibrium flowfields or techniques that show promise for studying 
these flows. Several classes of measurement techniques that are, or could be, used to study 
hypersonic nonequilibrium flows are not described herein. For example, surface measurement 
techniques such as pressure and temperature sensitive paint, phosphor thermography, skin 
friction measurements and photogrammetry (for model attitude and deformation measurement) 
are excluded to limit the scope of the chapter. Physical probes such as pitot probes, static 
pressure probes, heat flux gauges, total temperature probes, and others are excluded to limit 
scope. We further exclude measurement techniques that require particle seeding because of 
practical difficulties in seeding particles into hypersonic nonequilibrium flow facilities and 
because the particles may not faithfully follow the flow through steep gradients 6 often associated 
with hypersonic nonequilibrium flow. However, particle based methods may still be useful in 
certain hypersonic nonequilibrium flow applications. The chapter emphasizes quantitative 
techniques, so excluded are techniques that typically provide only flow visualization, such as 
schlieren. 

Hypersonic facilities present many challenges to the implementation of optical measurement 
techniques. Facilities capable of producing high-enthalpy flow often have hazards that must be 
accommodated in order to make measurements. Such hazards include those associated with high 
pressure, explosives, high voltage, high acoustic noise and toxic gas generated by the facility and 
the measurement system. Consequently, in many facilities the operator cannot be near the 
facility during operation so the equipment must be operated remotely or must operate 
autonomously. To reduce costs and increase test section size and operating pressure range many 
hypersonic facilities operate for very short duration, for example 100 msec or less. These are 
known as impulse facilities. For laser systems operating at 10 Hz only one measurement can be 
acquired in these facilities per tunnel run. Since some impulse facilities can operate only once 
per day it is difficult to obtain a statistically significant number of measurements with some 
measurement techniques. In such facilities higher-speed measurement systems that can obtain 
data at kHz or MHz rates are desirable. Large impulse facilities can move a few centimeters 
during operation, complicating measurement technique alignment. The flow through longer-run- 
duration facilities induces vibrations that can misalign optical measurement systems and changes 
in the model position may occur via thermal expansion of the support system. Furthermore, 
many large hypersonic facilities have poor or no optical access; if there is optical access the 
window material may not be suitable for the desired measurement technique. Hypersonic 
facilities are often located in dirty operating environments and for practical considerations the 
lasers often are placed far from the test section requiring long path lengths resulting in 
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misalignment and reduced laser energy delivered to the test section. Finally, some hypersonic 
facilities produce very luminous flow, making difficult the detection of weak light signals. In 
such flowfields, coherent optical measurement techniques are desirable. 

In some instances, it is desired to study nonequilibrium effects. However, nonequilibrium 
also frequently occurs as a byproduct of heating and accelerating the test gas in hypersonic 
facilities. For example, in order to achieve the required high flow enthalpy and Mach number, 
gas is often heated to thousands of Kelvin and then expanded through a nozzle. Energy can be 
added to the flow using electricity (via convection over a heating element or by an electric arc), 
by driving shockwaves through the gas, or by combustion. ' Upstream of the nozzle, the gas can 
be vibrationally excited and may reach chemical equilibrium in which minor species such as NO, 
O, and N are present for the case of air. When the gas expands through the nozzle, the collision 
rate is not sufficient to keep the chemical composition in equilibrium and the chemical 
composition effectively freezes at an intermediate composition. Similarly, vibrational energy 
modes do not stay in equilibrium with rotational and translational modes, which can result in T vi b 
> T rot ~ T i runs- This test gas in a state of chemical and/or thermal nonequilibrium, expanded to 
the desired Mach number, then passes over the test article being studied. An accurate 
understanding of the nozzle exit, ox free stream conditions, is required to properly characterize 
the test stream, calibrate the facility, interpret the test results and to scale the test results to flight. 
For a detailed description of the different types of hypersonic test facilities see the chapter titled 
“ High Enthalpy Facilities and Plasma Wind Tunnels for Aerothermodynamics Ground Testing ” 
in this manuscript. 

The next section of this chapter reviews some of the equipment commonly used to perform 
the different measurement techniques. The following sections introduce and develop the theory 
behind some of the most important measurement techniques used to study hypersonic 
nonequilibrium flows. Theoretical developments are limited to introductory material with 
references provided for more detailed analysis. Then, selected measurement examples are 
described. These examples do not provide an exhaustive review but are meant to illustrate 
typical applications. When possible, applications to hypersonic nonequilibrium flows are shown. 
However, examples of other flowfields are provided if the measurements demonstrate the 
potential for application to hypersonic nonequilibrium flows. 

2. Equipment Used for Molecular-Based Optical Measurement 

There are similarities in the hardware used in most optical spectroscopic measurement 
systems. Some of the more common lasers and detection systems will be summarized here. 
References for commercially available equipment described are not cited to avoid preferential 
treatment to certain vendors. It is difficult to capture all possible types of lasers and optical 
detection systems in a few paragraphs. This section is intended to provide an introduction and 
overview to the uninitiated and is not a comprehensive discussion. The various specifications 
described are approximate and are presented only to provide rough orders of magnitude. 
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Laser systems can be classified as continuous wave or pulsed, though hybrid systems also 
exist. Continuous wave (CW) lasers operate continuously and are specified by their operating 
power in Watts. CW lasers can be made to operate with very narrow spectral resolution suitable 
for high-resolution spectroscopy where the spectral width of the laser is negligible compared to 
the spectrum being measured. Inexpensive, low powered diode lasers having power on the order 
of 1 mW are used for absorption measurements. Higher powered CW lasers on the order of 10 W 
are used for Rayleigh scattering measurements. These higher powered CW lasers are generally 
not broadly tunable, so for resonant techniques like laser induced fluorescence only coincidental 
spectral overlaps between available laser frequencies and atomic or molecular transitions can be 
exploited. To perform spectroscopic measurements at arbitrary wavelengths, high-powered CW 
lasers can be converted to tunable light using dye lasers and/or other energy conversion systems 
such as doubling or mixing crystals, sometimes in combination with other CW lasers. 
Currently, high-resolution (i.e. laser linewidth much narrower than the typical absorption 
linewidths) CW lasers are available commercially from the ultraviolet (UV) to the infrared (IR). 
The main limitations of CW lasers are their relatively low power output which may require long 
integration times to acquire data with sufficient signal-to-noise ratio in some applications. 
Whereas absorption spectroscopy, which detects the full laser beam, can be performed at 
hundreds of kHz, scattering approaches like the Rayleigh, Raman and LIF that use CW lasers 
typically acquire data several orders of magnitude slower because of low signal levels. For 
comparison, a typical hypersonic flow with a gas velocity of 2000 m/s, a 1 mm resolution would 
require 500 ns time resolution, corresponding to 2 MHz data acquisition rate. 

Diode lasers provide CW light in most wavelength ranges of spectroscopic interest from the 
ultraviolet to the far infrared. Wavelengths can be varied over a limited range by modulating the 
temperature or injection current of the diode, at up to MHz scan rates, 9 providing extremely high 
time resolution for wavelength- scanning measurements. Because the lasers and detectors are 
small, have no moving parts, and the lasers are often coupled to optical fibers, these devices can 
be used to produce relatively robust measurement systems that are particularly well suited to the 
often harsh environments associated with nonequilibrium hypersonic flows. These 
characteristics have made diode laser absorption spectroscopy a particularly desirable technology 
for hypersonic flight testing programs. 10 

o 

Pulsed lasers typically producing much higher power (on the order of 10 Watts) for short 
durations, typically over a few nanoseconds, have been commercially available for about three 
decades. Such powers are usually required to perform nonlinear spectroscopic measurement 
techniques like CARS. These lasers are usually specified in terms of the energy per pulse. For 
example, frequency doubled Nd:YAG laser systems can produce 1 J per pulse at 532 nm. Again, 
this type of laser is not broadly tunable but it can be used to pump a dye laser or optical 
parametric amplifier and together with other conversion devices (doubling and/or mixing 
crystals) can produce high-energy tunable output from the UV to the IR. Ten ns duration pulsed 
lasers allow measurements to be made with time resolution that “freezes” the flow (prevents 
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blurring). Raman, Rayleigh, PLIF and CARS systems based on Nd:YAG lasers typically 
integrate signals over 10 to 500 ns time durations. Unfortunately, commercially available high 
energy pulsed laser systems typically operate at repetition rates of 10 to 100 Hz, which according 
to the calculations above, are too slow to resolve unsteady flow fluctuations. Therefore, 
individual measurements can be considered instantaneous, but successive measurements are not 
correlated in time. The linewidth of pulsed lasers is typically broader than high resolution 
narrowband lasers because of a number of factors. A consequence of the Heisenberg uncertainty 
principle is that the shorter the duration of a laser pulse the spectrally more broad the laser pulse 
will be. Consequently, the spectrally narrowest 10 ns pulses typically have line widths equal to 
0.004 cm" 1 . These pulses are routinely obtained by injection seeding pulsed Nd:YAG lasers so 
that only a single longitudinal mode is amplified. These injection-seeded pulses are much 
broader than commercially available CW Nd:YAG lasers. Without seeding, the same Nd:YAG 
laser would have an even larger linewidth of about 1 cm" 1 . Nd:YAG-pumped narrowband dye 
lasers typically have line widths within the range of 0.03 to 0.5 cm" , though modifications such 
as injection seeding or intercavity etalons can be used to reduce the linewidth of a pulsed dye 
laser. Conversely, spectrally selective elements can be removed to create a broadband (red) dye 
laser with a linewidth of 20 nm (~ 500 cm" 1 ) or more. 11 

There are several new laser technologies that have been developed in the past decade with 
high output power and fast repetition rates. This allows for flow freezing, temporally resolved 
measurements to be obtained with good signal-to-noise levels. Pulse burst laser technology uses 
Nd:YAG lasers that have been optimized to be pulsed 25-100’s of times in succession at 
frequencies up to 1 MHz. Pulse burst laser technology is not yet widely commercially available 
although it is being brought to market. Another relatively new technology that produces high 
power and high repetition rates is known as ultrafast or femtosecond lasers. As the name implies, 
the duration of the laser pulses can be nearly a million times shorter than a typical nanosecond 
laser. While the energy per pulse (typically less than 10 mJ, for commercially available systems) 
is much lower than ns-based systems, the output power can be on the order of 10 10 Watts, 
opening up new opportunities for coherent measurement approaches. Such high powers can be 
efficiently converted using optical parametric amplifiers and doubling or mixing crystals to 
extend the tuning range of femtosecond lasers to the UV and IR. Femtosecond lasers are 
commercially available at repetition rates up to 10 kHz - still two orders of magnitude too slow 
to time-resolve unsteady hypersonic flows. Even though the resulting measurements would not 
be time correlated, such high data rates are especially useful for acquiring statistics in short 
duration or expensive-to-operate test facilities. The higher data rate allowed by femtosecond 
lasers would be beneficial in most experiments. Picosecond-duration lasers are also available for 
spectroscopic measurement applications though picosecond lasers have neither the high energy 
per pulse of nanosecond lasers nor the high peak power of femtosecond lasers. The linewidth of 
a femtosecond laser is closely related to the duration of the pulse. Typical off-the-shelf 
femtosecond lasers with pulse duration on the order of 100 fs have linewidths of ~10 
nanometers. While this relatively large linewidth prevents conventional high-resolution spectral- 
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domain spectroscopy, the time domain can instead be used to obtain detailed spectral 
information. 

Most measurement techniques use some type of spectrally selective optic or optical system to 
collect or process signal. The simplest spectrally selective optic consists of filters that are placed 
in front of detectors or cameras. Absorptive and reflective filters are commercially available (or 
can be custom manufactured at a reasonable cost) at many wavelengths, for example to block the 
laser’s wavelength and collect the signal in another wavelength. In some applications, the 
narrower the spectral filter the better, since the signal is often very close to the laser’s 
wavelength such as in Rayleigh and rotational Raman scattering. Spectral filters blocking a 10- 
20 nm full width half maximum (FWHM) wide band near the laser’s wavelength with greater 
than 6 orders of magnitude of attenuation while transmitting >95% of the Raman shifted signal 
can be purchased. However, such filters do not provide sufficient resolution to observe Doppler 
broadening and Doppler shifts required for gas temperature and velocity measurement from 
Rayleigh scattering, for example. Gaseous filters, such as h vapor contained in a heated cell, can 
provide spectrally narrow (0.0015 nm or 0.05 cm 1 FWHM) and high (6 or more orders of 
magnitude) absorption which is useful in many applications. Etalons are spectrally dispersive 
optics or optical systems that are made of a pair of partially reflective mirrors. An etalon’s 
resolution depends on its flatness and reflectivity; the resolution can be high enough (0.01 cm" 1 ) 
to resolve the Doppler broadening (0.05 - 0.5 cm" 1 ) in typical hypersonic flow experiments and 
have an advantage over gas filters in that they can operate at any wavelength. Spectrometers, 
containing a dispersive optic such as a prism or grating, are typically used to disperse a signal 
beam into a spectrum that can be acquired by a detector or camera. Small (100 mm focal 
length), easy-to-use fiber-optic coupled spectrometers are commercially available for 
applications in which spectral resolution on the order of 0.02 nm (5 cm 1 ) is acceptable. Larger 
spectrometers (for example with 1 meter or longer focal length) provide higher resolution spectra 
- the resolution generally increasing linearly with the focal length. Coupled with a CCD camera, 
such spectrometers can typically achieve resolutions of 0.15 cm" /pixel. Resolution can 
sometimes be improved by using a camera with smaller pixels or by using a telescope at the exit 
of the spectrometer to enlarge the spectrum. 14 

For spectroscopic measurement techniques, light is detected with a variety of different single 
point, line, or imaging devices. Among single point detectors, photodiodes are typically used in 
applications where the intensity of the collected light is high, for example in absorption 
spectroscopy. Photodiodes are used because they are linear over a wide range of light intensities, 
have fast time response, require relatively simple electronics and are inexpensive to produce. 
Photomultiplier tubes (PMTs) are usually single -point detectors but they allow variable gain 
useful for amplifying light at low signal levels, while also providing fast (sub-ns) time response. 
While photodiodes and PMTs are both available as linear and even two-dimensional arrays, other 
camera technologies are usually used in spectroscopic applications where high resolution spatial 
or spectral information is required. Charge coupled device (CCD) cameras provide the linear 
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response (less than 1% nonlinearity), high dynamic range (12-16 bits) and high sensitivity 
required by many spectroscopic imaging or spectral applications. For applications requiring short 
exposures, for example to freeze the flow or to reject unwanted luminosity while capturing a 
transient signal, CCD cameras can be fitted with an image intensifier allowing exposure times 
down to hundreds of picoseconds. CCD cameras having 1 million pixels or more typically frame 
at tens of Hz consistent with high-powered YAG laser systems. CCD camera technology has 
been extended to much higher speeds (up to MHz rates) albeit with lower spatial resolution (for 
example 160x160 pixels). Recently, CMOS camera technology has allowed megapixel cameras 
to frame at tens of thousands of Hz allowing new opportunities for optical measurement 
techniques. 

3. Optical Emission Spectroscopy 

Optical emission spectroscopy (OES) is one of the oldest diagnostic techniques and has been 
applied to various fields of research. Traditionally, astronomy was one of the first fields to apply 
this passive method. Due to the large distances to even the nearest stars, emission spectroscopy 
continues to be one of the main techniques for understanding the stars. Despite its disadvantages, 
OES presents an often simple experimental method to access the thermochemical state of 
emitting gases. The experimental set-ups are relatively simple (not requiring a laser, for 
example) and the underlying optical principles well understood. However, the scope of this 
chapter does not allow for a complete description of emission spectroscopy methods applied to 
plasma diagnostics, even if restricted to nonequilibrium applications. Instead, selected examples 
of successful applications to nonequilibrium environments will be shown to demonstrate 
fundamental approaches to access these complicated problems. 

In general, emission spectroscopy (if properly calibrated) can provide absolute, quantitative 
information on excited states (electronic, vibrational, or rotational). However, even in strongly 
radiating environments, the majority of particles will be in the ground state or in weakly excited 
levels. Therefore, an extrapolation of the measured quantities to the low lying states is necessary 
to obtain information about the gas chemistry. In Thermodynamic Equilibrium , the upper state 
population distributions are linked to the ground state through known relations, given by 
Boltzmann distributions. Even if temporal or spatial changes occur in a medium, an equilibrium 
assumption may be justified if the variations are slow enough that, for any point, one can assume 
thermodynamic equilibrium in some neighborhood about that point. This condition is widely 
called Local Thermodynamic Equilibrium (LTE). When nonequilibrium processes are involved, 
these relations may no longer hold. Hence, a combination with other methods (LIF, CARS, 
absorption spectroscopy, etc.), which can access the ground and exited states not covered by 
emission spectroscopy, may better describe the thermodynamic state of the plasma. However, 
even without knowledge of the relation to the total particle density of the species under 
consideration or under nonequilibrium conditions, basic information may be extracted from 
spectral features. In particular, the line shape of individual lines carries important information on 
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translational temperature which often can be determined from Doppler broadening or the 
electron density can be inferred from the line’s Stark broadening. Line ratios may be used to 
determine rotational, vibrational, or electronic excitation temperatures as opposed to the absolute 
intensity of the detected emission. Ratios between lines from inert elements (e.g. argon) and 
reacting components such as oxygen may carry information on chemical processes, even in 
nonequilibrium situations, without the need to compute total species densities. 

One major difficulty of emission spectroscopy measurements is that each measurement 
provides integrated intensities along the entire line of sight. Furthermore, absorption along this 
line cannot be separated from single measurements. In case of an optically-thin medium, local 
distributions can be obtained through tomographic methods or through deconvolution of vertical 
profiles (e.g. through Abel-inversion or related methods if rotational symmetry of the emitting 
volume can be assumed). In the general case of an absorbing medium, only indirect methods 
made through a comparison with computed radiation based on numerical simulation are feasible. 

Emission spectroscopy is typically performed in atmospheric entry wind tunnel tests for at 
least two different reasons: 1) to determine radiation as a load factor for design purposes (typical 
examples are given by recent testing in the NASA Ames EAST facility), 15 and 2) to use radiation 
as a tool for diagnostic purposes, e.g. for facility characterization 16 or to gather information for 
comparisons to numerical simulation (e.g. material response). One further important application 
of emission spectroscopy is flight measurements which can provide validation data for 
theoretical simulation and extrapolation to flight of ground test data, either by flight 
experiments ' or by airborne or ground based observation. ’ ’ In industrial plasma 
applications (e.g. coating or edging), OES may be used for process control through empirical 
correlations or even to provide information on the underlying physics. 

The more recent development of off-the-shelf miniaturized spectrometers and the availability 
of fiber-coupled instrumentation have simplified emission instruments. The hardware 
requirements can be roughly grouped into systems with 1) low resolution for survey and 
overview applications, 2) more sophisticated instruments (usually with larger focal length) for 
higher spectral resolutions of specific spectral features such as molecular band structures, and 3) 
high resolution systems (e.g. Fabry-Perot interferometers) for detailed investigation of single 
spectral lines in particular for the investigation of line shapes and broadening mechanisms. 

The following sections will give a short introduction to basic physics for emission of 
radiation of atoms and molecules followed by examples for plasma diagnostics to determine 
temperatures, velocities or even particle densities. A description of some simulation codes with 
particular respect to hypersonic flows and examples for successful simulation of equilibrium 
plasmas with these codes will also be given. In a second part, the effect of nonequilibrium on 
plasma radiation will be discussed with examples showing emission spectroscopy applied under 
nonequilibrium conditions. 
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3.1 Basic Physics Describing the Emission of Atoms and Molecules 

Unlike the emission of thermal radiation of solid bodies, which yields continuum radiation 
following Planck’s formula, the emission of flames and plasmas usually consists of discrete 
emission lines. If measured in lower spectral resolution these lines appear as bands. These lines 
and bands originate from bound-bound transitions (i.e. in both upper and lower states of the 
transition, the participating electron is bound to an atom or molecule). These line spectra may 
overlap with continuum radiation produced through free-free or bound-free radiation generated 
through interaction of electrons with each other or heavy-particles (non-electrons) which, in the 
latter case, involves ionization and recombination processes. In general, the prediction of single 
emission lines can be divided into three problems: 

a) determination of the spectral position of the line using empirical or quantum mechanical 
information on the governing energy difference, 

b) determination of the line strength using thermodynamic information about the plasma 
state (i.e. the upper state population), and 

c) determination of the spectral shape of the line determined by broadening processes, again 
from the thermodynamic state (e.g. the gas temperature for Doppler broadening and 
pressure/species information for pressure or/and Stark broadening). 

For bound-bound emission processes, each emission line will be produced through a 
radiative transition from an upper energetic state £3 to an energetically lower state Ej. The 
resulting photon will carry the energy difference AE=E 2 -Ej which determines its frequency 

v=AE/h (3.1) 

with Planck’s constant h. This frequency is equivalent to the speed of light divided by the 
wavelength of the emitted radiation, in a vacuum: 

v=c<A (3.2) 

Often, the energy is used and tabulated normalized by the product of (h Co), given in the unit 
[cm 1 ]. Consequently, the spectral position of an emission line might be characterized by its 
wavenumber, v , which is the inverse of its wavelength, 
again often given in [cm 1 ]. 

The energy emitted by one emission line, often called 
the volumetric emission coefficient, can be written as the 
energy of one photon (hv) times the number of atoms 
which undergo this transition per unit time, defined as a 
probability (Einstein coefficient of spontaneous emission, 

A 21 ) times the population density in the upper energy state, 
emitted into full space of 4n: 



Figure 3.1. Energy level diagram for 
spontaneous emission. 
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We recognize that the emitted power depends on the product of the transition probability A 21 , 
which only depends on the quantum mechanical characteristics of the energy level £ 2 , and the 
upper state particle density m which depends on the thermodynamic state of the emitting 
medium. For a computation of the emitted power, we therefore have to determine the quantum 
mechanical relations for the species and the energy level of concern and the thermodynamic state 
of the emitting gas. 


First, let us consider the emission of atoms, which then will be extended to molecule 
emission. For atoms, the only degree of freedom for excited states is given by electronically 
excited levels. These excited levels are not arbitrary but characteristic for each species. In the 
simplest case of a hydrogen atom with only one electron, the energy levels can be calculated 
through the Rydberg formula, resulting in a series of spectral lines emitting at wavelengths, A: 


~ = R h 

A h 


1 1 


v n i 


n 


2 J 


(3.4) 


where R H is the Rydberg empirical constant. The energy levels can also be expressed in terms of 
the angular momentum of the electrons (by Planck) as: 


A E 




(3.5) 


who also expressed Rh is given by universal constants: 


Rh 


2n 2 k 2 e 4 m 
ch 3 


(3.6) 


For more complex atoms with more electrons, such laws cannot be easily derived and laws 
for allowed radiative transitions between electronically excited states must be extracted from 
quantum mechanics. Detailed descriptions of the nomenclature and the basic governing 
processes are found for example in Refs. 22 and 23, tables for transition probabilities for the 
different atoms can be found in Ref. 24 or in online databases, for example published and 
maintained by the National Institute of Standards and Technology (NIST). These transitions 
can be summarized in so called Grotrian diagrams. 26 

To determine the final emission coefficient, the number of particles in the upper excited state 
must be known. In general, this could be done by balancing all populating and depopulating 
processes, such as collisions (excitation and quenching), radiative transitions (absorption and 
emission), and chemical reactions. In thermodynamic equilibrium, the forward and backward 


12 



processes equal each other and the plasma will achieve a stationary state in which the 
populations of the excited states, n 2 , are related to the total density of the species of concern, n (h 
by a Boltzmann distribution: 


o 

"2 ~ " : eX P 


Q(Tj 


kT 


ex J 


(3.7) 


with the excited state energy, E 2 , the excitation temperature, T ex , the multiplicity (sometimes 
called degeneracy) of the upper state, g 2 , and the temperature dependent partition function, Q: 


Q(T,J=Y,ki ex P 


A. 

v kT exJ 


(3.8) 


Under nonequilibrium conditions, the population balance for each individual state has to be 
determined separately by considering all individual processes. Common approaches for such 
situations will be discussed in the second part of this chapter. 

For molecular emission, the same principles used for atoms apply but additional degrees of 
freedom of rotation and vibration have to be considered. Internal energy now can be stored in all 
three degrees of freedom and a radiative transition may contain contributions of all three 
energies. Therefore, the relation for one photon becomes: 

AE = AE el +AE, lb +AE m =hv = j . (3.9) 

The following considerations are restricted to diatomic molecules, which are relatively well 
understood and also are the major molecular constituents of high-temperature air (N 2 , 0 2 and 
NO). The governing processes for poly atomic molecules become rapidly too complicated to be 
considered in this work. Although simplified models for diatomic (considering the molecule as 
two point masses with a rigid massless connection (rigid rotator) or as two point masses 
connected by a perfectly elastic massless spring (harmonic oscillator)) are sometimes used in 
computational fluid dynamics to account for the internal energies, these models are not sufficient 
for spectroscopic considerations, in particular at high temperatures. The molecular potential are 
usually described by an anharmonic oscillator where the energy approaches infinity as the 
nuclear distance decreases and approaches the dissociation energy with increasing intemuclear 
distance. The solutions of the Schrodinger equation yield discrete energy levels, characterized by 
the vibrational and rotational quantum numbers V and J. Figure 3.2 illustrates a molecular 
transition. 
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Figure 3.2. Illustrations of a) vibrational potentials of a harmonic oscillator and a realistic molecular 
potential, and b) a ro-vibrational transition between two electronic levels. 


Recently, ab-initio computations for molecular potentials were performed using quantum 
chemistry.” However, these models require substantial computational effort both during the 
initial calculations and in post processing. More commonly, approximations as proposed by 
Dunham ” are used where vibrational and rotational energies are approximated through 
polynomial series depending on the rotational and vibrational quantum numbers J and V. 
Usually, vibrational and rotational energies, respectively, are computed in the normalized form 
with the units [cm" ] with the spectroscopic constants derived from empirical observation:” 
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( 3 . 11 ) 


( 3 . 12 ) 

( 3 . 13 ) 


The coefficients for these series can be found in the literature. Detailed computation of 
different molecules may take additional effects into account such as spin splitting, lambda 
doubling, etc. A rather complete description was given by Herzberg,” while more current 

9Q 

molecular constants for air were published by Gilmore.” 
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Figure 3.3. Deformation of the NO ground state 
potential with rotational quantum number. ' 0 Modified 
and reprinted with permission of the author. 


The dependence of rotational energies 
on vibrational quantum numbers shown 
in Eq. (3.11) yields a deformation of the 
molecular potential with rotational 
excitation rather than just a vertical shift 
in energy as suggested by the illustration 
in Fig. 3.2. 

The resulting potentials show a potential 
well above the dissociation limit (and 
therefore stable excited molecules with 
energies higher than the dissociation 
energy)’ 1 as shown in Fig. 3.3. 


If a truncated Dunham series (3.10 
and 3.11) after second-order terms is used, Morse showed that an exact solution for the 
Schrodinger equation representing the motions of the nuclei in a diatomic molecule can be 
found. He proposed a molecular potential in the form: 


U{x) = D e { \-e- px ) (3.14) 

With the well depth of the potential D e , the difference x to the equilibrium distance r e , and the 
parameter [i that depends on the Eigen-frequency of the oscillator coo, the reduced mass of the 
molecule ju A and D e \ 32 


J3 = 1.2177 • 10 7 m — 

1*>. 


(3.15) 


This representation, however, becomes inaccurate for high quantum numbers and more 
sophisticated potentials may be used, e.g. the three parameter Fippincott potential, 33 the five 
parameter Hulbert-Hirschfelder potential, or Rydberg-Klein-Rees (RKR) potentials.' 

For the computation of line strengths, the definition of the volumetric emission coefficient 
remains but the transition probability becomes a function of the relation between the upper and 
lower state rotational and vibrational configuration, the indices ‘ and ‘ ‘ indicating the upper and 
lower state, respectively: 
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J’ is the rotational quantum number of the upper state, the electronic transition moment 
(R e (r vV , )) 2 q vV . describes the probability of a particular electronic transition depending on the 

s J 

inter-nuclear distance, the Honl-London factor — indicates how the total intensity of a 

2/'+l 

transition is distributed among the different branches (in the easiest case the P-, R-, and Q- 
branches defined by a change in rotational quantum number from the upper state to the lower 
state of -1, 1, or 0, respectively). 


The electronic transition moment is based on the Franck-Condon principle which postulates 
that the internuclear distance does not change during the transition. Therefore, the wave 
functions of the upper and lower level will yield a higher transition probability when the overlap 
of the upper and lower molecular potentials is larger as illustrated in Fig. 3.4. The transitions 
indicated by solid lines are probable because of the overlap of wave functions. The hypothetical 
transition indicated by the dashed line on the right side of Fig. 3.4 is not possible because of non- 
overlap of the wave functions caused by the different internuclear spacing in ground and 
electronic states. Both the electronic transition moment and the Honl-London factor depend on 
the molecular species of concern. 


a) b) 



Figure 3.4. Illustration of the Franck-Condon 
principle. Reprinted from Ref. 61 with permission 
of Springer Science and Business Media. 


In equilibrium, population of each energy 
level is distributed through a Boltzmann 
distribution. Although the definition of 
thermodynamic equilibrium postulates that 
only one temperature governs all distribution 
functions, separate temperatures for rotational, 
vibrational, and electronic excitation are 
already introduced at this point to be able to 
extend these considerations to different degrees 
of nonequilibrium. Therefore, the upper- state 
density now depends on three temperatures and 
can be described through a combination of 
three Boltzmann distributions of rotational, 
vibrational, and electronic energy: 
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with the partition function: 
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3.2 Line shape of emission lines and line broadening 

So far, each emission line was considered to be emitted at one discrete wavelength 
determined by the energy difference of the upper and lower different states. However, emission 
and absorption lines can be broadened through different processes (some of which are detailed in 
the next section on Absorption Spectroscopy) such as: 

• natural broadening due to the Heisenberg uncertainty principle (Lorentzian) 

• Doppler broadening due to the thermal motion of the emitting particles (Gaussian) 

• Collision broadening (Lorentzian) 

- Van-der-Waals due to collisions with neutrals 

- Resonance broadening due to collisions with the like species (perturber’s state 

connected by an allowed transition to the upper or lower state of the transition) 

- Stark for interactions with ions or electrons, specific for each transition/line 

• Zeeman effect due to interaction with magnetic fields (Lorentzian) 

Each of these processes actually causes a line shift for each photon. The integration over 
many photons with different shifts produces a broadening by distributing the transition energy 
over a certain wavelength (or wavenumber) range. Depending on the physical effect, the final 
line shape will follow a Gaussian or a Lorentzian profile. Usually, the broadening is described by 
the linewidth at its half maximum (full width, FWHM, or half width, HWHM). Half widths of 
Lorentz profiles can be summed up linearly, half widths of Gaussian profiles can be summed in 

quadrature: fwhm g = J FWHM ^ + FWHM 2 G , • If several processes occur simultaneously, the 

resulting broadening may be described by a Voigt profile. Spectral lineshapes are detailed in 
Section 4 of this manuscript, which focuses on Absorption Spectroscopy. 

3.3 Spatial resolution 

One major disadvantage of emission spectroscopy (shared with absorption spectroscopy) is 
the integrating character of the measurements. Even if the set-up is focused on one particular 
point in the field of view (e.g. the center line of a flow field), the detected signal is not limited to 
that point. In a radiating volume, the measured signal will result from an integration along the 
line of sight with all cross sections along this line contributing with equal efficiency as long as 
the local cross section defined by the limiting rays does not exceed the volume dimensions. The 
measured integrated intensity I(z) is coupled to the local emission coefficient e(r) by the Abel 
equation: 35 






dr 


(3.19) 
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where R is the outer radius of the measured volume, z the off-axis position, and r the local radius. 
Solving for s(r) yields 



(3.20) 


In the special case of axi-symmetry and an optically thin volume radiator, local emission 
values can be determined by analyzing a series of off-axis lines of sight. For the example of a 
homogeneously emitting medium, an elliptic distribution of the measured emission versus off- 
axis position will be obtained. Local values may be obtained by applying a so-called Abel- 
inversion. However, an analytical solution can only be found if the measured profile versus off- 
axis position can be described by particular functions (such as a Gaussian distribution, a pair 
polynomial function, or elliptic functions). In the general case, i.e. if the measured profile I(z) 
cannot be described by one of the above mentioned functions, approximate methods can be used. 
The most common one is the “onion peeling” method: 

Each measurement is interpreted as an average over a ring with constant emission coefficient 
measured at an off-set value zj. The emission along the lines of sight is then the sum of the 
individual emission coefficients <?, times the length Ljj with which the line of sight cuts each 
individual ring r, , yielding a rather simple system of linear equations. The geometric principle is 
illustrated in Fig. 3.5(a) for 4 measurement positions. The system for 4 lines of sight reads as: 
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and can even be solved recursively for small numbers of off-axis positions or through matrix 
inversion for a larger set of measurements. For this method, the individual lines of sight do not 
have to be parallel to each other, e.g. if the variation of z is realized through tilting the optical 
axis rather than moving the set-up, as illustrated in Fig. 3.5(b). To reach higher accuracy, radial 
distributions over each ring element might be assumed and included in the set of equations. 


18 





Figure 3.5. Illustration of the Abel-inversion principle for a set of off-axis measurements a) parallel to 
the rotational symmetric radiation field, and b) realized through tilting the detection set-up (from 
reference 36). Reprinted with permission of the authors. 

3.4 Analysis of emission spectroscopy data: 

The most traditional application (e.g. in astronomy, or process control in industrial 
applications) of emission spectroscopy is probably the identification of radiating species in the 
emitting medium. Even for an unknown thermodynamic state of the emitting medium, the simple 
occurrence of particular emission lines is proof for the existence of the element of concern. If a 
number of species contribute to the unknown spectrum, the occurrence of multiple lines of the 
same element can provide a conclusive identification. Since the emitted spectrum will change 
with temperature, some basic knowledge of the thermodynamic state might be required already 
at that stage of interpretation. For industrial applications, an engineering approach might be 
suitable where spectra are recorded at the desired condition, and only deviations from these 
spectra are used as indication for a possible failure of the process without even knowing the 
cause. 

If quantitative measurements of the chemical composition are intended, however, the 
thermodynamic state of the emitting medium has to be known to a significant accuracy. All 
atomic emission lines (and most molecular spectra in the visible to near-infrared region) will 
originate from electronically excited states. Quantitative interpretation of these spectra requires 
detailed knowledge of the population of these excited states and their link to the ground state of 
the species of concern which, even at high temperatures, usually contains the highest density. 
Without detailed collisional models, quantitative interpretation usually requires the validity of 
Boltzmann distributions, and therefore an assumption of equilibrium. In general, equilibration 
will be reached through collisional energy exchange. As described in the introduction, the 
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various forms of internal energies (rotational, vibrational, and electronic excitation) require 
different numbers of collisions to equilibrate. While some 10s of collisions are needed to 
equilibrate rotational excitation, vibrational equilibration already requires some 100s of 
collisions, and electronic excitation may be reached after several 1000s of collisions during the 
life time of the excited states. Therefore, equilibration of electronic excitation with the ground 
states is most likely achieved at high pressures. 

3.4. 1 Temperature determination methods 

While number densities usually require quantitative knowledge of the absolute measured 
emission, temperatures can be obtained from the qualitative spectral appearance, e.g. from 
relative line ratios or line shapes. At low pressures, the dominant broadening mechanism of 
emission lines is the Doppler broadening which is related to the translational temperature 
through equation 3.21. If the emission lines are measured with sufficiently high spectral 
resolution, a Gaussian profile can be fitted to the measured line shape and the extracted half 
width can be evaluated with respect to the governing temperature. The basic assumption is only 
that the concept of temperature is defined which implies that the thermal velocities are 
distributed according to a Maxwell distribution. No information about the strength of the 
measured line is needed. Typical half widths for temperatures up to 10s of thousands of Kelvins 
are on the order of 10s of picometers (~0.5 cm) at maximum. Such spectral resolution is usually 
obtained through interferometers (e.g. Fabry-Perot interferometers) or Echelle spectrometers. 
The main influences on the measured data are usually the instrument broadening of the set-up 
used and line broadening caused by other physical effects that are potential error sources that 
need to be modeled and corrected for. The higher the pressure, the more dominant collisional 
broadening becomes but collision broadening typically decreases with T 05 while Doppler 
broadening increases with T (from Eq. 3.11). For the individual lines of carbon dioxide 
(CO 2 ) or water (H 2 O), for example, the Doppler broadening will be dominant at 0.1 atm or below 
for temperatures higher than 1000 K, while at 1 atm Doppler and collisional broadening are on 
the same order of magnitude for temperatures of 3000 K (compare Fig. 10-7 in Ref. 37). 

If a significant amount of ionization occurs (e.g. in electrically-heated plasmas), Stark 
broadening through interaction of excited electronic states with the electrons may contribute to 
the line broadening. This can potentially interfere with temperature measurements or alternately, 
allow electron number density measurements. However, Stark broadening factors are different 

o o 

for each individual transition, so an appropriate choice of emission lines under investigation 
may provide useful temperature measurements, even at significant degrees of ionization. On the 
other hand, the broadening of selected lines (e.g. the hydrogen lines H a or Hp) may provide 
measurement of electron density. 15 ' 39 If no hydrogen is available in the original plasma 
composition (e.g. for air plasmas) it has been shown that even intentional seeding with small 
amounts of hydrogen enables electron density measurements at hypersonic conditions while not 
significantly altering the plasma state. 
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If lines originating from the same atom or molecule but different excitation energies are 
measured simultaneously, governing temperatures can be determined from ratios of the total 
emission of the individual transitions. This data is obtained by integrating the measured spectrum 
over the line shape (in the following referred to as line emission). If a Boltzmann distribution of 
the excited states can be assumed, the ratio of different emission lines can be analyzed by 
evaluating the distribution function with respect to temperature. If the thermodynamic state of 
the emitting medium is not known to be in equilibrium it is appropriate to analyze more than two 
emission lines to determine the temperature and, at the same time, prove the validity of the 
Boltzmann distribution using a so called Boltzmann plot. If the equation for the electronic 
excitation of atoms is written as 
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with I 2-1 being measured line intensities for upper state excited energies Ej. A plot of the 
logarithmic terms vs. excitation energy follows a straight line if the upper state densities are 
Boltzmann distributed. The agreement with this straight line shows the validity of a Boltzmann 
distribution, the slope contains the desired excitation temperature (slope = -1 /(kT ex )). The higher 
the difference in excitation energy, the better the accuracy of the measurement, a rule that will be 
discussed in more detail in the context of PLIF in section 4.5. 



Figure 3.6. Effects of the plasma state on electronic 
excitation temperature in comparison to equilibrium 
temperature, based approximately on data in Ref. 40. 


However, even if the Boltzmann 
plot indicates a Boltzmann distribution 
of the excited levels under 
consideration, the resulting temperature 
might not necessarily equal the 
equilibrium temperature. As illustrated 
in Fig. 3.6, in a recombining or 
ionizing plasma, the population of high 
electronic levels might be modified 
with respect to the equilibrium 
population through the chemical 
processes (overpopulation of highly 
excited levels in a recombining plasma, 
underpopulation in the ionizing case) 
yielding a Boltzmann-like distribution 
within the measurement accuracy with 
a governing temperature different from 
the equilibrium temperature. 
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If the same procedure as described above is applied to molecules, similar equations can be 
derived for rotational, vibrational, and electronic excitation. The individual temperatures can be 
determined by appropriate selection of evaluated emission lines (for T rot , all lines under 
consideration should originate from the same upper and lower vibrational and electronic levels, 
etc.). Figure 3.7 gives an example of a Boltzmann plot to determine the rotational temperature of 
Ni + in an air plasma from different rotational lines of the 0-0 transition of the No" 1 " 1 st Neg. 
System. 



0 5,000 10,000 15,000 

rotational energy, cnr 1 


Figure 3.7. Boltzmann-Plot for the determination of 
rotational temperature from rotational lines of the 0-0 
transition of the N 2 + 1 st Neg. system in an air plasma based 
on data from Ref. 30. Reprinted with permission of the 
American Institute of Aeronautics and Astronautics. 


A more common approach in the 
last few decades, however, to 
determine vibrational or rotational 
temperatures is a comparison of a low 
resolution molecular spectrum with 
theoretical simulation. The band shapes 
are usually characteristic for the 
rotational temperature, band ratios of 
different upper vibrational excitation 
levels may be used to determine 
vibrational temperatures. Again, this 
method does not necessarily require 
intensities calibrated to absolute values 
since only band ratios are used. It is, 
however, mandatory to at least perform 
a relative calibration within the 
wavelength range of concern, even if 
this range is rather narrow. 


3.4.2 Simulation Codes 

Over the years, several codes for the simulation of plasma emission with particular respect to 
the application to hypersonic flows have been developed worldwide. In the following, only line- 
by-line approaches are summarized. In these codes, emission, absorption, and broadening 
properties for each individual transition of the species of concern are calculated, and a final 
spectrum is composed through a superposition of all lines under consideration. In comparison 
with other methods (smeared band, narrowband, or wideband models) that are usually faster, line 
by line models still yield the most accurate simulation and offer the advantage of covering all 
ranges of spectral resolution (e.g. through simulation in high spectral resolution, followed by a 
spectral integration). Typically, lines of sight through a radiating volume with known 
composition are constructed and the emission is propagated through this volume. Typical 
applications are found in rebuilding experimental conditions in ground experiments or flight for 
a comparison of measured data. Another task is to determine radiation loads in ground test 
facilities or, even more important, in flight situations to determine final loads and support the 
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design process. A list of some important codes is given in Table 3.1. Again, it is not possible to 
cover all existing codes, so the list was compiled with the intention of being representative but 
not comprehensive. Typically, these codes are under continuous development to account for 
recent simulation demands (e.g. adding radiating species for atmospheric entries into other 
celestial bodies’ atmospheres, the most recent examples being Mars, Titan, or Venus) and new 
physics models being developed. 

Table 3.1. Representative (but not comprehensive) list of codes for the simulation of plasma radiation. 


Code 

(Organization, Country) 

Simulated species 

Radiation Transport 
capabilities 

Restrictions/availability 

Ref. 

Neqair 

(NASA Ames, USA) 

Air, H, PL, C, 
CN, CO z 

Line of sight, spherical cap, 
infinite slab 

EAR restricted 

41 

Hara 

(NASA Langley, USA) 

Air, H, PL, C, 
CN, C0 2 , Ar 

Line of sight, spherical cap, 
infinite slab 

EAR restricted 

42 

Spradian/Radipac 
(JAXA, Japan/South 
Korea) 

Air, H, PL, C, 
CN, C0 2 , Ar 

Line of sight, spherical cap, 
infinite slab 

none, 

available through JAXA 

43 

Parade 
(ESA, Europe) 

Air, NH, C, CN, 
CH, C 2 

Ray tracing through additional 
Fortran code (Herta) 

none, 

available through ESA 

44 

Specair 

(Ecole Centrale, France) 

Air, H, PL, C, 
CN, CO z , Ar 

None 

reduced version available 
for download 45 

39 

Lifbase 
(SRI, USA) 

OH, OD, NO, H, 
CN, CF, SiH, N, + 

None 

none, available for 
download 46 

47 

Spartan ( ipfh , Instituto 
de Plasmas e Fusao 
Nuclear, Portugal) 

Air 


freely distributed for 
academic and research 

48 

Hyperrad , (NASA 
Ames, USA) 

Air 


EAR restricted 

27 

CNRS-Rad Database, 
(Ecole Centrale, France) 

Air, CN, CO, 
CO + , C 2 

None 

Available as 
spectroscopic database 

49,50 


Note: “Air” is: N 2 , O z , NO. N, O. N 2 + 


3.5 Application examples of emission spectroscopy 

Below we present some typical examples of the application of emission spectroscopy to 
hypersonic flows in ground test facilities. There is a wide spectrum of facilities operational 
around the world which may be divided up into impulse and stationary facilities. Shock tubes 
belong to the first category. They typically recreate the hypersonic flight regime accurately but 
provide only short test times. Stationary facilities provide sufficient test time to yield a realistic 
material response for thermal protection system (TPS) qualification but only partially rebuild a 
true flight environment. Both types of facilities need plasma diagnostic methods to characterize 
the flow, emission spectroscopy being widely applied. 

In impulse facilities, imaging set-ups have been used in recent years to take a snapshot of the 
flow propagating through the test section while gaining simultaneous spectral information from 
the entire shock and post-shock system. 15 ' 51 Figure 3.8 shows such an image with extracted 
spectra and spatial intensity distributions across the shock and the post-shock region. 
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Shot 35 (Blue) - 8.53 km/s. 0.098 torr 


The purpose of recent work was to 
experimentally quantify radiation 
loads to spacecraft during 
atmospheric entry to reduce 
margins in the design process, 
mainly due to equilibrium 
radiation generated in the post- 
shock region. However, increased 
interest in the nonequilibrium 
radiation inside the shock and in 
the immediate relaxation region 
has been expressed recently. 52 
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A simple to interpret case of 
emission spectroscopy is plasma 
states in thermodynamic 
equilibrium. For most high- 
temperature air species, the 

simulation capabilities are 

sufficient to reproduce measured spectra. A good example is given by data taken at Stanford 

on 

University in an atmospheric, inductively coupled plasma torch by Laux et al. Figure 3.9 
shows measured and simulated emission spectra in that equilibrium case that are in excellent 
agreement. 


Figure 3.8. Example of data obtained from a single 
spectrograph in the EAST facility. Vertical and horizontal cross 
sections of the data shown depict the spatial and wavelength 
resolved data, respectively. 51 Reprinted with permission of the 
American Institute of Aeronautics and Astronautics. 
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Figure 3.9. SPECAIR: NO p, y, 8, s, P’, y’, 0 2 SR and atomic lines, including self-absorption within 
plasma and absorption in 5 m of room air in the low UV (left), and measured and computed LTE spectra 
of the nitrogen/argon plasma at the generator nozzle exit from the near-UV to the near infrared.’ 9 
Credentials to von Karman Institute for Fluid Dynamics, in Radiation and Nonequilibrium Collisional- 
Radiative Models by C.O. Laux, VKI LS 2002-07. 
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These measurements demonstrate the ability to make absolute concentration measurements using 
OES; in this case the absolute intensity of No emission matches the predicted emission assuming 
equilibrium chemistry. 

Sometimes changes in the gas 
dynamic state happen too fast for the 
plasma to reach equilibrium (both in 
the thermal and chemical sense). In 
this case, populations of the excited 
states may no longer follow well 
defined (Maxwellian) distribution 
functions. An analysis requires a 
detailed balance of population and 
depopulation for each excited state, 
yielding the necessity of collisional 
radiative models. The assumption of 
a quasi-steady-state, where the time 
change of the excited states is 
considered small against the 
population processes themselves, is 
often applied as described by Park 31 
and may hold over a wide regime of 
conditions. However, additional population processes may become dominant, in particular those 
through chemical reactions, requiring knowledge of excited-state chemistry. Evidence of excited 
state chemistry was found during the space shuttle development in the so called “shuttle glow”. 
This basic process was thoroughly described by Golde and Thrush. 53 Excited states are 
populated through chemical reactions that end up directly in these states without the need of 
collisional excitation from the ground state. The final radiating level might even be populated 
through overlapping potentials (usually of molecules) through level crossing processes. 54 An 
experimental realization was given by Laux et al. by rapidly cooling down an equilibrium 
plasma. The molecular emission of No from the 1 st Pos. System (N? B 3 n g -> N 2 A 3 Z U + ) was 
monitored and showed strong deviations from the simulation based on thermal excitation with 
significant overpopulations of high vibrational quantum numbers, as illustrated in Fig. 3.10. 

Indeed, these overpopulations could be related to the aforementioned level crossing processes 
(also known as inverse pre-dissociation) and the spectrum was rebuilt applying collisional 
radiative models. 54 Recently, the same processes have been shown to be dominant for the 
emitted radiation of No in highly expanding high-enthalpy flows used for TPS material testing at 
NASA by Winter et al , 36 as demonstrated in Fig. 3.11. The measured spectrum can only be 
rebuilt by introducing significant overpopulation factors for high upper vibrational quantum 
numbers (V u =8 to V u =13) of the emitting electronic state of No. It was shown that the molecular 



Figure 3.10. Measured spectra and computation assuming 
LTE and non-LTE simulation of a nitrogen/argon plasma 
after rapid cooling from an equilibrium state from Ref 39. 
Credentials to von Karman Institute for Fluid Dynamics, in 
Radiation and Nonequilibrium Collisional-Radiative Models 
by C.O. Laux, VKI LS 2002-07. 
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radiation (at least of these high vibrational levels) is not produced by excitation of No ground 
state molecules but by recombination of N atoms into electronically excited No levels. This 
provides a potential method for measuring ground state atom concentrations in nonequilibrium 
flows without the need of much more complicated techniques such as two-photon laser-induced 
fluorescence (TALIF). 


b) 


£ 4.0E-04 


2.0E-04 

I 


vibrational overpopulation 



Exp. Acton-Pimax 

— Exp. Ocean Optics 

sum over scaled vib. levels 

— 5.5E+06x Vu 13 
2.3E+06 x Vu 12 
1.3E+06 x Vu 11 
4.0E+03 x Vu 2 
2.3E+03 x Vu 1 


tjjj 


\ ■ 

HUH 

AM 


Hr- 

/ 


H 1 i\fl 
/( 1.4 

/ 

Wfl ■ 


650 700 750 

wavelength, nm 



Figure 3.11. Measurements of the N 2 1 st Pos System in the free stream of the IHF arc -jet facility in 
comparison to LTE simulation with NEQAIR and superposition of individually scaled spectra simulated 
separately for each upper vibrational level (left), 36 and potential crossing between the N 2 A' 5 I g + , and N 2 
B 3 n g states 54 (right), which is reproduced with permission from the authors. Right hand figure is 
Copyright 1988, AIP Publishing LLC. 


4. Absorption Spectroscopy 

Absorption spectroscopy has been used to study nonequilibrium flows since at least 1961. 55 
In their textbook on Gas Dynamics, Vincenti and Kruger describe how Camac (1961) used the 
technique to measure vibrational relaxation times behind normal shock waves in oxygen-argon 
mixtures. 56 Camac used a xenon-filled microwave-excited resonance lamp, a hot cathode 
hydrogen discharge tube, and a high-voltage Lyman discharge lamp to produce the ultraviolet 
light needed for the absorption measurements. Absorption spectroscopy has been since used to 
extract other quantitative thermodynamic and fluid dynamic properties flows such as 
temperature, pressure, species concentration, and velocity. 

4.1 Absorption Theory 

Quantitative measurements are achieved with absorption through the analysis of the fluid’s 
interaction with light. As light, with spectral irradiance I v [W/cm /cm’ ] propagates through an 
infinitesimal thickness dx [cm] of an absorbing medium of concentration C [molecules/cm ] with 

2 57 

a molar absorption coefficient 8 [cnr/molecule], its intensity decreases according to: 
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til,, = I „sCdx 


(4.1) 


Using a spectral absorption coefficient k v [cm ], which is equivalent to the product of e and 
C, Eq. 4.1 can also be written as: 58 ' 59 


- dl r = k v I v dx . 


(4.2) 


Integration of this equation gives the well-known Beer-Lambert law: 


- [ k v dx 

= e J 


(4.3) 


In this equation, 7 0 is the incident irradiance prior to absorption. If the medium can be specified 
by a single temperature, velocity, and absorbing species concentration and the pressure is 
uniform and absorbs linearly over the path length L, the relationship can be simplified to: 
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(4.4) 


The parameter k v contains key information about the absorbing medium, including the 
population of potential absorbers and the probability of absorption. Strategies that exploit the 
relationship between this information and thermodynamic properties lead to various quantitative 
measurement techniques. Contributions to k v can be better explained with the aid of an energy 
diagram. Figure 4.1 shows incident light being absorbed, which causes molecules to transition 
from an energy ground state (E\) to an energy excited state 
(£ 2 ). The energy of the incident photon, which is the product 
of Planck’s constant h and the frequency of the light v, is 
transferred to the excited molecule. 


Not all potential absorbers in the ground state are excited. 
The rate of absorption is governed by the transition line-shape 

function Y v [cm] and the transition line- strength factor 3 ] 2 [cm" 

2 

]. These two components contribute to the spectral absorption 
coefficient: 60 



Figure 4.1. Two-level energy 
diagram of absorption process. 


K=s 12 y v . 


(4.5) 


If we assume that the population in the upper state is negligible, S 12 is proportional to the 

population of absorbers in the ground state Ni, Einstein’s coefficient of stimulated absorption Bn 
1 2 - 1-1 

[s' (W/cm'/cm" )" ] and the energy of the incident photon. Y v describes the shape of the transition, 
spectrally, and is related to the probability that molecules will absorb photons at various 
frequencies. By definition, the transition line- shape function is integrated to give: 
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(4.6) 


\Y v dv = 1 • 

Light sources typically have a discrete spectral width, which further complicates the process. 
The spectral irradiance of light can be calculated from: 

I V =IL V (4.7) 

where / [W/cm ] is irradiance and L v [cm] is the line-shape function of the light source. Similar 
to Y v , L v is defined as: 


J L v d v = 1 . (4.8) 

To account for these two different line-shape functions, the integrated absorption coefficient 
is used: 


k m{ =S n G, (4.9) 

where G is the spectral overlap integral defined as: 

G = \Y v L v dv (4.10) 

The spectral overlap integral represents the amount of intersection between the light source 
line-shape and the transition line-shape. There are two extreme cases to consider. 59 The first case 
assumes a monochromatic light source and a relatively broad absorption transition. Here, G 
approaches the value of Y v at the light source’s frequency. A spectral scan of a relatively broad 
transition using a narrow-band laser to determine the absorption line-shape is an example of this 
case. The second case assumes a monochromatic absorption transition and a broadband light 
source. Here, G approaches the value of L v at the absorption transition’s frequency. Note that Eq. 
4.10 can be used to calculate G for all cases. For single-mode diode lasers the approximation of a 
monochromatic laser is frequently used. This is because the laser line-width is of the order of 10s 
of MHz while transitions are of the order of a few GHz at standard atmospheric conditions. 
However, at low-density conditions, such as those appropriate to flight in the upper atmosphere, 
the monochromatic approximation may no longer be valid. In these cases, corrections such as 
Eq. 4.10 need to be used. 

4.2 Broadening Mechanisms 

The transition line-shape is influenced by broadening mechanisms, which have various 
thermodynamic dependencies. Typically, the spectral width and shape of the transition are 
affected by natural (Av N ), collision (Av c ), and Doppler-based (Av^) broadening mechanisms. 
Natural and collision broadening are homogenous where each absorber has the same transition 
energy and frequency. Natural broadening is related to the Heisenberg uncertainty principle, as 
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the frequency of the transition cannot be perfectly known. Natural broadening is usually small in 
comparison to the other broadening mechanisms. 61 Collision or pressure broadening occurs from 
physical collisions with other molecules. Not only do collisions increase the spectral width 
(Av c ), they also cause a shift (Av v ) in the center frequency and asymmetry in the line-shape. For 
a multi-component gas mixture, these values can be calculated from: 62 


a L = P X(^ 2 ^) 

i 

(4.11) 

A ^, = p Yj(x< d i) 

(4.12) 


l 


where P, /■,, 2y h and 3/ are the pressure, mole fraction, collision width per unit pressure, and 
collision shift per unit pressure, respectively, of the absorbing molecules. The terms y,, and 3/ 
have temperature dependencies which, for many species of interest to absorption spectroscopy 
can be evaluated with the High Resolution Transmission (HITRAN) database. 63 As we shall see, 
the temperature can be determined through the area ratio of two transitions. Because the effects 
of pressure are the same for both transitions, these effects cancel in the integrated signal ratio 
allowing an unambiguous temperature measurement. If the temperature is known, then the 
analysis of the collision shift and/or collision broadening can be used to calculate the gas 
pressure. Homogeneous broadened transitions have a Lorentzian line-shape function: 


A ^ H 1 

ln (v-v 0 f +(A v h / 2f 


(4.13) 


where A vh is the spectral width of a homogeneous broadened transition (full width half 
maximum). 

Doppler broadening is inhomogenous, where each absorbing molecule has a transition 
frequency that is Doppler shifted according to its velocity component in the laser propagation 
direction. This is because of thermal motion, which induces Doppler shifts. Doppler broadened 
transitions have a Gaussian line-shape function: 
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(4.14) 


where vo is the center frequency corresponding to a stationary absorber. 

When both homogeneous and inhomogeneous broadening mechanisms contribute to the line- 

co 

shape function, the Voigt function is used. The Voigt function is a convolution of the 
Lorentzian and Gaussian profiles: 64 
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VM=T , f v dy 

n J -°°<r +(x-y) 

(4.15) 

where 


a = Vln2 A Vh 
Av d 

(4.16) 

x = 2 Vln2 ~ F ° f 
A V D 

(4.17) 

Physically, the use of the Voigt profile can be justified by each collection of molecules or 
atoms in a particular Doppler subgroup being subject to the effects of collisional broadening, and 
this is expressed through the convolution operation. The Voigt function V(a,x) cannot be solved 
analytically but is usually calculated using numerical methods. Figure 4.2 illustrates the 
different line shapes computed using an approximation to Eq. 4. 15. 65 Figure 4.3 shows the 
effect of Lorentzian and Gaussian broadening on the spectral shape of a line using the fraction of 
the total transition energy contained within a given frequency difference from the line center. 
With Gaussian broadening dominant, 99% of the transition strength is contained within two half 


widths at half maximum (Av HWHM ), whereas for a Lorentzian profile, this value is not reached 

for ~50 half widths. In other words, a strong Lorentzian-broadened line might significantly 
influence distant spectral regions if Lorentz broadening is dominant. 



(v-v 0 )/Av hwhm 

Figure 4.2. Spectral line shape for Lorentz 
(collision), Gaussian (Doppler) and Voigt broadening 
for equal transition strength and half-width. 



Figure 4.3 Fraction of transition energy 
contained between (v-v 0 ) for a Lorentzian 
(collision) and a Gaussian (Doppler) line profile. 


30 






4.3 Translational Temperature 

If A vd can be measured, then the translational temperature T trans can be determined by the 
following relationship: 59 


= 2v 


2kT 


(MW )c 


-In 2 


(4.18) 


where k is Boltzmann’s constant and MW is the molecular weight of the absorber. If other 
broadening mechanisms are small in comparison and a spectrally narrow light source is used to 
spectrally scan across the transition, then A vd can be extracted with relative ease to determine the 
translational temperature. Even in flows with considerable homogeneous broadening the 
translational temperature can be determined from the Doppler component of the Voigt line- 
shape. 66 Using linewidth to determine translational temperature has the advantage of only 
requiring a single absorption line. This is particularly useful for absorption measurements of 
atomic species where there may only be one transition available within a laser’s output range, or 
for lasers that have a very limited tuning range. It is also necessary for measurements in cases 
where there is nonequilibrium between rotational and translational energies. Measurements of 
translational temperature based upon transition width may have higher uncertainties compared to 
measurements of rotational or vibrational temperature based upon a Boltzmann plot. 67 This is 
because of the square-root dependence of the Doppler width on temperature and the additional 
challenge of separating the Doppler width from other broadening mechanisms, especially at high 
pressures, when measuring translational temperature. 68 For example, an uncertainty of 1.5% in 
the measurement of the Doppler width in A-band Doppler-broadened oxygen lines at 296 K leads 
to an uncertainty in temperature of 3%. Zhang et al. 69 measured the translational temperature in 
an arcjet facility using absorption of Argon at 81 1 nm. The Doppler width was measured with an 
uncertainty of 12% after the effects of homogeneous broadening, Stark broadening and noise 
were accounted for, leading to an uncertainty in the translational temperature of less than 30% 
over temperatures from 2000 K to 10,000 K. The difference in sensitivity between line-strength- 
based temperature measurements and Doppler width measurement depends strongly on the 
difference in ground-state energies between the available transitions. 

Another example of a translational temperature measurement is described in the work of 
Parker et al where absorption experiments were performed in the Large Energy Shock Tunnel 
(LENS I) at the Calspan-University of Buffalo Research Center (CUBRC). Quantum cascade 
lasers were used to obtain absorption measurements of NO deep into the infrared spectrum. 
Figure 4.4 shows a sample absorption spectra from that work near 5447 nm (1835.56 cm 1 ), 
which corresponds to one of six fundamental NO rovibrational transitions used in that study. 
Using the Doppler width contribution based on the fitted Voigt line-shape, a translational 
temperature of 190 K was computed (Eq. 4.18). This temperature, measured in the test section, 
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corresponds to a lower enthalpy condition (5 MJ/kg) and using a helium driver gas. At high- 
enthalpy conditions, the LENS I facility produces free-stream flows that are characterized by 
chemical and thermal nonequilibrium. In the work, translational temperatures were compared to 
rotational temperatures and to CFD predictions at these conditions. Although the measurement 
uncertainty was not reported, the authors predict a temperature resolution of 50 K. 



Figure 4.4. Absorption spectra measured in the CUBRC LENS 1 facility used to calculate translational 
temperature. Reproduced from Ref. [70] with permission of the authors. 

4.4 Rotational Temperature 

Rotational temperatures are typically measured by analyzing multiple absorption features that 
have line- strengths with different rotational temperature dependencies. Usually either multiple 
light sources are used to excite separate isolated transitions or a single light source is used to 
spectrally scan over adjacent absorption features. This technique is referred to as multi-line or 
ratio thermometry. The spectrally integrated absorbance that is measured experimentally is 
defined as: 71 
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The subscript j refers to the absorption feature that is probed. Often a modified version of the 
spectral absorption coefficient is used: 72 


k v =P Xi SY v 


(4.20) 


The subscript ‘12’ in the line strength term S has been dropped for clarity. With this change, 
a modified version of the line-strength factor S [cm 2 atm 1 ] is used: 73 
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where T rot , T 0 . Q, and E mt are the rotational temperature, reference rotational temperature (296 
K), the total internal partition function, and the total rotational energy of the molecular quantum 
state absorbing the light, respectively. AE is the difference in energy between the upper and 
lower states of the transition. It should be noted that the final term in the equation is due to 
stimulated emission, and is negligible at temperatures below approximately 2500 K and at 
wavelengths below 2.5 pm. 74 

If the light source is spectrally narrow, the integrated absorbance for each j transition can be 
simplified: 


aj = \ kUlv j P Xl LSY v dv =P Zi LS j \Y v dv (4.22) 

If the entire transition can be resolved using a scanning technique, then the integrated 
absorbance becomes: 


a 


= P X ,LS j 


(4.23) 


The ratio of two separate absorption features, denoted by subscripts j= 1 and j= 2, with 
different rotational energies, becomes: 


a \ 


a 2 



S 2 (T 0 ) 



(4.24) 


where A E rot is the change in lower state rotational energy between the two absorption lines. 
Because many terms are cancelled in this ratio, it is possible to extract T mt . Note that it is also 
possible to calculate this ratio based on peak absorbances instead of spectrally integrated values 
to extract temperature if the line-shape function is included. This assumes that the pressure shift 
and broadening coefficients are similar for both transitions. 

Wehe et al., performed rotational and translational temperature calculations from EBO 
absorption measurements in the high-enthalpy (10 MJ/kg) Calspan 96-inch Hypersonic Shock 
Tunnel using a diode-laser sensor. A probe, which housed the critical optical components, was 
placed directly into the core of the flow. In order to improve the signal-to-noise ratio, the 
hypersonic air flow was seeded with 6% water vapor. Absorption measurements, shown as solid 
data points in Figure 4.5, correspond to an FFO absorption transition (V1+V3 band) near 1396 nm 
(M) and 1400 nm (1|). A Gaussian fit, shown as a solid line in the figure, was applied to the 
experimental data to determine the spectrally integrated absorbances and Doppler broadening of 
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each transition. The ratio of the absorbances was used to compute the rotational temperature (Eq. 
4.24) 


The Doppler broadening was used to 
calculate the translational temperature (Eq. 
4.18). The measurement uncertainty of the 
rotational temperature measurement was ± 
13 K (±2.3% of mean) whereas the 
measurement uncertainties of the 
translational temperatures based on the two 
absorption profdes (kj and k 2 ) were ± 15 K 
(± 2.7% of mean) and ± 35 K (± 6.4% of 
mean), respectively. 

These high-resolution line-shapes were 
recorded every 125 psec (8 kHz), allowing 
for time -resolved temperature 

measurements. Figure 4.6 shows the time- 
resolved rotational and translational 
temperatures extracted during the facility 1 
thermal equilibrium, which was expected ba 



Relative Frequency (cm ') 

Figure 4.5. Absorption spectra and fit through H 2 0 
transitions in Calspan 96-inch Hypersonic Shock 
Tunnel. Adapted from Ref. [75] with permission of the 
authors. 

1 time. Agreement indicates that the flow was in 
1 on equilibrium calculations reported in the work. 



Figure 4.6. Translational and rotational temperature measurements in the Calspan 96-inch Hypersonic 
Shock Tunnel. Adapted from Ref. [75] with permission of the authors. 

It can be difficult to accurately measure the line-strength due to problems in determining the 
zero absorbance intensity near the absorption feature wings, which is needed for baseline 
levels. 76 Using diode lasers, for example, the laser intensities can drift over time due to thermal, 
mechanical, and electronic effects. Liu et al. discuss how low signal levels, baseline-fitting 
errors, beam steering effects, and others can affect the estimation of the line-shapes. 76 To combat 
these problems, wavelength modulation spectroscopy (WMS) with second harmonic (2 f) 
detection can be used. Detailed theory about WMS can be found in the literature. ' Through 
the superposition of a modulation of frequency / onto the light source, the problem of excessive 
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DC noise that occurs at low absorbances can be greatly reduced. The absorbing medium converts 
some of the wavelength modulation into amplitude modulation of the laser intensity. The 
amplitude modulation occurs at the modulation frequency / and at its integral harmonics, nf. 
Then, by using phase-sensitive electronics, a selected harmonic can be isolated (e.g. n = 2). At 
sufficiently high frequencies / excess noise from the light source (e.g. laser diode) can be 
reduced. 

Goldstein et al. performed WMS-2/ rotational temperature measurements using a tunable 

01 

diode laser absorption sensor in the University of Virginia’s supersonic combustion facility. In 
this work, the facility was run with water injection to simulate vitiation effects and with 
hydrogen combustion (equivalence ratio of 0.33). Figure 4.7 shows sample experimental direct 
absorption (a) and WMS-2f/lf absorption (b) line-shapes of FTO obtained at the combustor exit. 




Time, ms 


Figure 4.7. Direct absorption (a) vs WMS-2f/lf (b) in UVa supersonic combustion experiment. 
Reproduced from Ref. [81] with permission of the authors. 


The two absorption features in Fig. 4.7(a) allow for the ratio thermometry technique to be 
used. However, the small magnitude of the neighboring feature results in a higher measurement 
uncertainty compared to the WMS-2/ technique. Both techniques give temperatures that agree 
reasonably well with CFD predictions (Table 4.1). However, the direct absorption technique has 
consistently larger measurement uncertainties compared to the WMS technique. The authors 
discuss the discrepancy in temperatures observed in the H 2 -air combustion case. 
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Table 4.1. Comparison of direct and WMS absorption calculations of rotational temperature. Reproduced 
from Ref. [81] with permission of the authors. 


Operation Mode 

Expected Value 

DA 

WMS 

11% Steam 

700-1000 K 

776 ± 10 K 

742 ± 9 K 

(Combustor Entrance) 

1 1.4± 0.2 % HO 

10.9 ± 0.1% H O 

2 

10.8 ±0.1% H O 

9% Steam 

700-1000 K 

860 ± 30 K 

831 ±9K 

(Exit Plane) 

9 ± 0.2 % H O 

9.1 ±0.2% H O 

9.1 ±0.1% H O 

12% Steam 

700-1000 k 

875 ± 50 K" 

850 ± 6 K" 

(Exit Plane) 

12 ± 0.2 % H O 

12.1 ±0.5% HO 
2 

11.5 ±0.1% HO 

9 

El -Air Combustion 9 = 0.33 

1800-2200 K 

1802 ± 94 K 

1765 ±41 K 

(Exit Plane) 

13% HO 

12.8 ±0.5% HO 

13.3 ±0.3% H O 


An alternative to wavelength modulation that has also been used with success in high-speed 
flows is ratiometric detection. ’ For this technique, the laser beam is split into a signal beam 
that propagates through the flowfield and a reference beam that does not. The photocurrents from 
the two beams are passed through a log-ratio amplifier, providing a direct measurement of the 
absorbance. If the photocurrents for the two beams are equal, the log-ratio detection removes the 
intensity modulation due to varying current on the beams, and the common-mode laser intensity 
noise, leaving only the absorption signal on the detector. The balanced variant of this technique 
includes an additional low-frequency feedback circuit that ensures that the photocurrents 
balance, zeroing the baseline of the measurement. This method, because it uses two beams rather 
than one, is optically more complex than wavelength modulation spectroscopy but electronically 
simpler, and provides similar absorbance detection limits. The output of the circuit is a signal 
that is directly proportional to the absorbance, making the analysis of the data straightforward. 

Upschulte el al. used the balanced ratiometric method to measure temperature, mole 
fraction and velocity in a supersonic combustion facility at the US Air Force Research 
Laboratory at Wright-Patterson Air Force base at Mach 2.1, using balanced ratiometric detection 
of water vapor near 1390 nm. The temperature and mole fraction measurements were made to 
within 100 K and around 0.3 % respectively, while beam steering adversely affected the velocity 
measurements. Griffiths and Houwing have also made measurements of water vapor 
temperature, mole fraction and velocity in a model scramjet engine at a flight Mach number of 
10, in the Australian National University’s T3 free-piston shock tunnel. They used two time- 
multiplexed distributed feedback diode lasers and balanced ratiometric detection at scan rates of 
up to 20 kHz. Temperatures ranging from 600 K to 1800 K were measured in the turbulent, 
combusting environment of the scramjet engine. 

4.5 Vibrational Temperature 

Similar to multi-line thermometry for rotational temperature, vibrational temperatures can be 
measured by probing absorption features at different vibrational levels. Pan and Oka measured 
vibrational temperatures of ArH + ions in a glow discharge by separately measuring the integrated 
absorbance of transitions in two vibrational levels. 86 They compared the vibrational temperature 
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to the translational kinetic temperature measured using the Doppler width of the transition and 
quantified the thermal nonequilibrium in the plasma. 

Brandt and Roth measured the absorption in the nonequilibrium flow behind a normal shock 
wave in a shock tube, in mixtures of CO and Ar, using a lead-salt IR diode tuned to absorption 
transitions lying between 2100 and 2200 cm" . Once again, the vibrational temperature was 
compared with the translational temperature to determine the degree of nonequilibrium 
downstream of the shock wave. The measured translational and vibrational temperatures 
corresponded well to the values predicted by Landau-Teller theory. 


4.6 Species Concentrations 

Similar to the thermometry techniques, species concentrations can be determined through the 
analysis of the integrated absorbance. If the temperature can be determined from one of the 
thermometry techniques, then the concentration of absorbers C, of species i can be determined 

oo 

from the integrated absorbance (Eq. 4.23) and ideal gas law: 
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(4.25) 


where R [J-mol-K" 1 ] is the gas constant. Concentration can also be extracted from analysis 
of the broadening mechanisms, which have pressure and temperature dependencies, using a 

on 

spectroscopic model that is fitted to the experimental data. 


Lin et al . , obtained CO concentration measurements using TDLAS in a Martian simulant gas 

oo 

(70% CO 2 and 30% No) that was accelerated in a shock tube. The flow immediately behind the 
6.31 km/s shock wave was in thermal nonequilibrium. Figure 4.8 shows a sample absorption 
spectra taken at 50 kHz near 2335.8 nm (4281.2 cm" 1 ) corresponding to the targeted CO 
transition. Even with excessive noise in the wings of the line-shape, a good Voigt fit was applied 
to the data, which is shown as the red solid line. The spectrally integrated absorbance was used 
to calculate the concentration. A room temperature cavity cell of CO was used to validate the 
technique. The authors report an average CO concentration of 7.46 x 10 “ molecules/cm behind 
the shock wave (Eq. 4.25). Errors in the Voigt fitting process resulted in a 1% uncertainty in the 
spectrally integrated absorbance values, which result in a 10% error in the CO concentration 
measurement. 


37 





Figure 4.8. Sample absorption measurement with Voigt fit (left) and nonequilibrium temperature 
measurements (right). Reproduced from Ref. [88] with permission of the authors and Springer Science 
and Business Media. 

4.7 Velocimetry 

Doppler-shift based velocimetry is based on the analysis of the spectral shift that occurs in an 
absorption feature when gas molecules have a component of velocity that is in the same direction 
as the excitation source. From the absorbing molecule’s frame of reference, the excitation source 
appears either red- or blue-shifted, depending on the molecule’s velocity. As a result, the spectral 
frequency at which the light beam is absorbed is either blue- or red- shifted, respectively, from 
the laboratory frame of reference. If the Doppler shift can be measured experimentally, usually 
through a spectral scan, then the gas velocity can be calculated from the following relationship: 

^Doppi er shift (26) 

^LightSource 

The velocity can be determined by passing the laser beam through a duct at two different 
angles and measuring the absorption of a spectral feature for each of the two transmitted beams, 
as shown in Figure 4. 9. 90 The Doppler shift in the peak of the absorption feature can be used to 
directly measure the velocity using Eq. 4.26. 
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Figure 4.9. (a) Typical arrangement for diode la 
absorption features for the two beams. Reproduced 



13083.6 13084.0 13084.4 13084.8 

Frequency, v [cm -1 ] 

absorption measurement, and (b) Doppler-shifted 
n Ref. [33] with permission of the authors. 


Infrared DLAS velocimetry measurements from CO, NO, and H 2 O molecules were 
performed in the S4MA, F4, and DLR piston driven HEG facilities. 91 In that work, enthalpy 
conditions ranged from ~1 MJ/kg to 20 MJ/kg and test times ranged from 100 s to 1 ms. In 
addition to velocimetry, thermometry (translational), and concentration measurements were also 
performed. Figure 4.10 shows a sample NO absorption line-shape (R 3/2) near 5366 nm (1863.6 
cm ') obtained in the DLR facility. The Doppler shift, which is used to calculate velocity (Eq. 
4.26) is visible in the figure. The integrated absorbance and Doppler width used to calculate 
concentration and translational temperature are also shown in the figure. 



Figure 4.10. NO absorption line indicating how velocimetry, concentration, and translational temperature 
are calculated. Reproduced from Ref. [91] with permission of the authors. 

Figure 4.11 shows DLAS velocimetry measurements in the F4 facility at an enthalpy 
condition of 15 MJ/kg. The velocimetry measurements are compared to CFD calculations using 
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both equilibrium and frozen assumptions. The velocimetry measurements show a closer 
agreement to the frozen calculations, indicating nonequilibrium. 
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Figure 4.11. DLAS Doppler-shift based velocimetry measurements in the F4 facility (50 MPa; 15 
MJ/kg). Reproduced from Ref. [91] with permission of the authors. 

The first oxygen velocity measurement in a high-speed flow was made by Philippe and 
Hanson using wavelength modulation detection at scan rates of 10 kHz in the A -band of 
molecular oxygen. Velocities were measured to within ±50 m/s, and agreed with theoretical 
estimates to within 5 %. More recently, velocity has been measured in a Mach 6 nozzle flow, 
using absorption in the same wavelength range, excited by a vertical cavity surface-emitting laser 
(VCSEL) emitting near 760 nm. ' The laser was able to be current-scanned over four oxygen 
transitions. Velocities of around 950 m/s were measured with an uncertainty of ±40 m/s. These 
measurements were made by propagating the laser beam through two opposing retro -reflectors, 
imposing opposite Doppler shifts as the beam propagates in the two opposing directions. This 
method allows the velocity to be measured using a single laser beam. 

4.8 Tomography 

Absorption tomography involves using multiple absorption measurements, which are path 
averaged, to reconstruct a spatial distribution (e.g. two-dimensional) of the measured property. 
Hanson suggested that new low-power fiber-based absorption sensor arrays may bring 
tomographic imaging inside aeroengine and airframe models for hypersonic ground tests. 94 
Tunable diode lasers have features that make them especially suitable for tomography. They can 
scan at very high rates, providing information at tens or hundreds of kHz, and the lasers are 
typically routed to and from the flowfield using optical fibers, simplifying the experimental 
setup. Technology from the telecommunications industry, such as fiber couplers and 
multiplexers, can be used to split the laser light into many channels and/or to combine many 
wavelengths into a single fiber, allowing multiple species to be probed with multiple colors 
along many paths through the flow. The light from the fibers is coupled into lenses, propagated 
across the flow and captured by lens-coupled multi-mode fibers, separated into wavelengths and 
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directed onto detectors. Thus, different species and spectral lines can be probed at different 
locations in a ducted flow, for example, allowing some spatially-resolved concentration and 
temperature information to be determined from the line-of-sight measurements. 


Recent advancements in TDLAS tomography include systems that scan a single laser or 
multiple lasers spatially across the flow and from different directions to allow time averaged 
tomographic reconstruction of temperature and species concentrations in turbulent combusting 
flowfields 95 or a grid of fixed optical paths transmitting across the flow, allowing a low- 
resolution but time -resolved tomographic reconstruction of concentration and temperature. 96 ’ 97 
Figure 4.12 (from Ref. 97) shows temperature and mole fraction maps obtained in a General 
Electric J85 gas turbine engine. Thirty individual laser beams — 15 in the horizontal direction and 
15 in the vertical direction, making a square mesh of 225 grid points — were used to probe the 
flow. The laser was rapidly scanned over multiple water vapor absorption lines near 1.35 pm 
with scans occurring at a rate of 50 kHz. While the smaller turbulent length (and time) scales of 
the flow are not resolved with this system, the large-scale spatial (and temporal) distribution of 
the temperature and mole fraction are resolved. 
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Figure 4.12. Temperature (left) and water vapor mole fraction (right) measured at the exit of a General 
Electric J85 engine, with 50 kHz diode laser tomography, adapted from Ref. [97]. The paper online has 
links to play animated movie clips of these data. Reprinted with permission of the authors and publisher. 

4.9 Cavity Ring-down Spectroscopy (CRDS) and Other Cavity-enhanced Methods 

Even when using sensitive detection techniques such as wavelength modulation or balanced 
ratiometric detection, it can be challenging to measure absorbances below 10' 5 . As many 
nonequilibrium flows occur at low number densities, and many important absorbing species such 
as oxygen have only very weak accessible absorption features, it becomes necessary to 
investigate even more sensitive optical systems. CRDS is a highly sensitive technique that is 
based on increasing the amount of absorbance in a medium by increasing the effective path 

go 

length (L) through the use of an optical system. For example, if an absorbing medium is placed 
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between two highly reflective dielectric mirrors, then the energy of light resonating between the 
mirrors will decrease with each pass. A small amount of light passing through one of the mirrors 
is detected and monitored as a function of time. The “ring-down time” is the time it takes for the 
light to decrease to an intensity threshold. It is also equivalent to the mean time that any photon 
will survive between the two mirrors (i.e. within the cavity). Analysis of the absorption features 
using a CRDS setup can lead to quantitative measurements similar to those described for direct 
absorption (e.g. velocity, temperature, concentration). Though it is more sensitive than 
conventional absorption spectroscopy, it is more difficult to set up and to maintain alignment. 
Scherer et al. describes the history, development and application of CRDS to pulsed supersonic 
flows." Similarly, Romanini et al., performed CRDS measurements in NOi in a supersonic slit 
jet using a narrow line, continuous wave laser source. 100 

CRDS, due to its extreme sensitivity and the requirement of locking a laser to a cavity, can 
prove difficult to employ in the often-hostile environments where nonequilibrium aerospace 
measurements need to be made. Cavity-enhanced absorption (CEA) techniques have been 
developed that provide enhanced absorption using an optical cavity for enhancing the effective 
path length without requiring complex cavity-locking mechanisms. Two of the more commonly 
used CEA variants are cavity-enhanced absorption spectroscopy (CEAS) 101 and integrated cavity 
output spectroscopy (ICOS). 102,103 These methods measure the output signal of the cavity as the 
laser is scanned across a transition, rather than measuring the ring-down at a discrete set of 
frequencies. The noise immunity of CRDS is traded for a simpler and more robust arrangement, 
while still having enhanced sensitivity compared to TDLAS. All cavity-enhanced systems are 
sensitive to degradation of the quality of the cavity mirrors, which can be a consideration when 
used in harsh environments. 

Most studies involving cavity-enhanced techniques in high-speed or nonequilibrium flows 
have concentrated on cold supersonic jets and plasmas. Bakowski et al., have used CEAS to 
make measurements of translational and rotational temperatures of No in an inductively coupled 
plasma, by exciting transitions from the A 1 E„ + (v=0) band near 686 nm. 104 Fiedler et al., have 
also expanded the CEAS method to the use of an incoherent broadband light source, and have 
performed measurements in a cold supersonic jet of azulene over a 40-nm wavelength range. 105 

4.10 Hardware Used for Absorption 

The most common excitation sources used for absorption measurements in hypersonic 
nonequilibrium flows are semiconductor diode lasers that provide tunable light in the infrared. 106 
Tunable diode lasers are relatively inexpensive compared to dye -pumped lasers and have narrow 
linewidths, which is necessary to resolve the absorption features in many common atoms and 
molecules. By adjusting the current through and the temperature of the diode, the wavelength can 
easily be controlled. Furthermore, the wavelength can be scanned at high frequencies (kHz), 
which is useful for some of the above-mentioned techniques. As a result of these advantages, the 
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diode laser is the preferred choice for absorption measurements although many other laser and 
non-laser sources are used. 

4.11 Strengths and Weaknesses of Absorption Methods 

There are many advantages of absorption methods over other optical diagnostics. In contrast 
to emission spectroscopy, absorption techniques can probe nonluminous flows, including both 
ground and excited states. Furthermore, optical setups for absorption experiments are simpler to 
set up and interpret compared to other techniques. Light sources, such as diode lasers, are 
relatively small and inexpensive compared to large pulsed laser systems needed for many 
techniques like PLIF and CARS. The frequency and scanning rate of diode lasers are also fast 
compared to Nd:YAG pumped dye lasers, allowing for high frequency temperature and 
concentration measurements to be obtained. Because of its simplicity, the absorption technique is 
also very robust. Also, unlike PLIF, there are no problems with quenching for absorption 
measurements, making them much simpler to analyze. 

The main weaknesses of absorption spectroscopy are that measurements are path-averaged 
and are normally restricted to a point normal to the flow. Multiple absorption measurements in 
combination with tomographic reconstruction techniques can be used to invert the data and 
produce spatial distributions of the desired properties. In contrast, PLIF can easily obtain two- 
dimensional measurements along a thin sheet that interrogates the flow. If there is axial 
symmetry in the flow, then Abel transforms (described in the emission spectroscopy section) can 
be used to obtain radial distributions. 

Because absorption is a path-averaged technique, uniform flow conditions across the path are 
desired for accurate measurements. Many practical flows have non-uniformities along the optical 
path, such as boundary layers, heat transfer at walls, and gas dynamic phenomena. These non- 
uniformities cause problems in the traditional signal analysis described above. In addition to the 
tomographic approaches described above, work reported in the literature describes methods to 
correct for these problems. ' Sanders et al., describes a method to deduce non- uniformity 

7 1 

along the path by probing many absorption transitions. 

5. Laser-Induced Fluorescence (LIF) 

Laser-induced fluorescence (LIF) uses a laser to probe individual species within the 
flowfield, providing information pertaining to both the thermodynamic (pressure, temperature, 
mole fraction) and fluid dynamic state (velocity) of the gas. The laser can be focused to a point 
for LIF measurements, formed into a thin sheet using a cylindrical and focusing lens in 
combination for planar LIF measurements (PLIF), or used to illuminate a volume for three- 
dimensional or stereoscopic imaging. Reviews of the laser-induced fluorescence measurement 
technique are available from Eckbreth 5 and others . 109 ' 1 10 The measurement technique works by 
inducing a transition, usually of an electron, from a lower energy state (Ef) to an excited energy 
state (£2) via stimulated absorption of one or more photons in the atomic or molecular species of 
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interest. In a two-level model assumption, the 
atom or molecule of interest in the E 2 state 
then returns to the E 7 state by transferring 
energy via spontaneous emission of a photon 
(fluorescence), or by transferring energy non- 
radiatively through a collision with another 
atom or molecule (collisional quenching). 

For an atom, such as N or O, the energy 
required to induce an absorption transition of 
an electron to the £? state from the ground (£7) 
state is equal to the energy difference between the atom with an electron occupying the excited 
electronic orbital and the atom’s ground electronic orbital configuration, respectively. For 
molecules, such as N 2 and 0 2 , the energies of the E 2 and £7 states include the energies associated 
with the vibrational and rotational motion of the molecule in addition to the energy associated 
with the molecule’s electronic configuration. Figure 5.1 shows a generalized two-level model for 
fluorescence with the stimulated absorption transition induced by a single photon. 

In Fig. 5.1, the rate at which the absorbing medium in the £7 state transitions to the £? state is 
proportional to the product of the Einstein coefficient for stimulated absorption, Bn, and the 
laser’s spectral irradiance (power per unit area per unit frequency), I v — IL V , where I is the 
irradiance (power per unit area) and L v is the laser spectral profile or laser line- shape (per unit of 
frequency ). 5 The energy of the absorbed and emitted photons is the frequency, v, multiplied by 
Planck’s constant, h. The Einstein A 21 and B 12 coefficients describe the probabilities for emission 
and absorption , 1 1 respectively, while Q 21 is the collisional quenching rate. For absorbing 
species, the line-shape function, Y v , describes the spectral width for a particular energy level. 

As detailed in Section 4.2, the line-shape function combines effects from Gaussian-shaped 
Doppler broadening and Lorentzian-shaped homogeneous broadening mechanisms. The integral 
of the product of the absorption line-shape function, Y v , and the laser’s spectral profile, L v , is 
defined by the overlap integral, G = f Y v L v dv , and describes what portion of a particular 
absorption transition is affected by the incident laser radiation. The rate constant for stimulated 
absorption, W/ 2 , which describes the rate at which species in £7 transition to £ 2 via absorption of 
a single photon 5 ' 110 is then given by: 



Figure 5.1. Two-level model of single-photon 
fluorescence. 


W 12 = B 12 IG (5.1) 

As the laser passes through a flowfield, it is absorbed at a rate corresponding to Eq. 5.1, 
inducing a transition between the £7 and £? states. Consequently, as the laser continues to 
propagate through the flowfield, the irradiance is continually diminished as a result of the 
absorption process as described by the Beer-Lambert law (Eq. 4.3). 
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In certain circumstances, the energy 
separation between the Ej and E 2 states for the 
species being examined with LIF is such that 
the frequency of a photon necessary to induce 
an absorption transition is far into the ultra- 
violet portion of the electro-magnetic 
spectrum. However, achieving frequencies far 
into the UV with conventional laser systems 
can be difficult. 111 Additionally, strong 
absorption of UV radiation by many materials 
and gases (including air) below approximately 
200 nm limits the application of laser systems capable of producing such frequencies. s One 
solution is to use two-photon LIF techniques to probe species having absorption transitions in the 
deep UV. With a two-photon technique, the frequency (and therefore energy) of each photon is 
half that required by a single photon for the same transition. Figure 5.2 shows a generalized two- 
photon fluorescence energy model for two-photon LIF. For atomic species, such as N or O, 
spontaneous emission is often observed between E 2 and an intermediate energy state, Ej. 

For the two-photon LIF process in Fig. 5.2, the two-photon absorption cross-section results in 
a rate constant different from that in Eq. 5.1 and has the form: 5,1 1 u 12-115 

W 12 = C 12 I 2 G (5.2) 

where Ci 2 relates to the two-photon absorption cross-section. Note that the dependence of the 
rate constant scales with the square of the incident laser irradiance, /. This behavior arises from 
the probability of observing the simultaneous arrival of two photons, which is the square of the 
probability for the arrival of a single photon. 111 

In both Figures 5.1 and 5.2, Q 22 is the so called quenching rate constant. It describes the rate 
at which energy is transferred through non-radiative collisions between excited atoms or 
molecules in the E 2 state and atomic or molecular collision partners of species i. This rate 
constant is computed similarly to that in Ref. 5 as: 

Q 2 i = N T 'LiX i 0s,iVs,i (5-3) 

where N r — P/k B T is the total population of the excited state, X L is the mole fraction of 
quenching species i, o si is the collision cross-section between the excited species, 5, and 

quenching species, i, and v si — ^Qk B T trans /nfj. si is the mean relative velocity 116,117 between the 

excited species and quenching species, with q s . being their reduced mass and T trans their 
translational temperature. The spontaneous emission coefficient, A, in Figures 5.1 and 5.2 is also 



Figure 5.2. Three-level model of two-photon 
fluorescence. 
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known as the Einstein A coefficient and describes the probability for spontaneous emission of a 
photon by an atom or molecule in the excited state. 111 


For the two-level model in Fig. 5.1, the rate of change of the populations Ni and N 2 , in the Ej 
and E 2 states, respectively, can be obtained from relations similar to those presented in Ref. 5: 

dJ t = - dJ t = -W 1 W 12 + N 2 (Q 2l +A 21 ) (5.4) 

N 1 + N 2 = N s = Xs f B N T (5.5) 

Equation 5.5 is a conservation law saying that the combined populations of the Ej and E 2 states is 
equal to the initial population of the excitation species’ Ej state, Ns, which is in turn the product 
of the species mole fraction, x s , the total population, N T , and the temperature-dependent 
Boltzmann fraction, f g . Equation 5.4 assumes that the laser intensity is sufficiently weak such 
that stimulated emission (W 2 /) can be neglected. 


The Boltzmann fraction, f g , describes the ratio of the number of absorbers initially 
occupying E t relative to all possible energy states at a particular temperature, T. when a system is 

1 IV 5 1 IV 118 

in thermodynamic equilibrium. This term has a general form given by: ' ' 



9 i e~ E ^ T 
1)9 f E i /kBT 


(5.6) 


where g . is the degeneracy of state j. Degeneracy refers to the number of quantized states that 

1 IV 

exist in a given energy level, Ej. The summation in the denominator is termed the partition 
function, Q as described in Eq. 3.18. 


More involved differential rate equations can be formulated for the two-photon method like 
those presented in Ref. 115. It should be noted that the models presented in Figures 5.1 and 5.2 
neglect transitions from Ej to E 2 resulting from collisions, with a rate constant of Qi 2 - This is 
usually a good assumption for large energy separations associated with the visible or UV 
transitions typically employed for LIF. Transitions from the E 2 and Ej states resulting from pre- 
dissociation ( Qp r e ) and ionization ( Q ion ), which are described in Refs. 5 and 110, have also been 
neglected. 


5.1 Linear, Steady State Solution 


If a continuous laser source is used to populate the E 2 state and detection of fluorescence 
occurs well after this source is turned on, then N 2 can be assumed to have reached its steady- state 
value. From this assumption, the left-hand- side of Eq. 5.4 is set to zero, resulting in two 
algebraic equations (Eqs. 5.4 and 5.5) for two unknowns (Ni and N 2 ). The same assumption can 
be made to determine the population, N 2 , achieved by a pulsed laser source if the time required to 
reach steady state is short compared to the duration of the pulse. This population is computed as: 
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N X s f B N T W 12 

2 W12+Q21+A21 


(5.7) 


The product of this excited state population and the spontaneous emission rate constant is 
N 2 A 21 and represents the number of transitions per unit time and per unit volume. Integrating this 
constant value with respect to time gives the total number of transitions per unit volume during 
the detection period. Substituting Eq. 5.1 for W/ 2 , and assuming that Wn is small compared to 
A21 and Q21 (which is valid for low-intensity excitation), a relation for the total number of 
photons collected via fluorescence by the detection device is: 110 

Slip = Xsf B N T B 12 IGA>t det V £77 (5.8) 

where <t> = A 21 / (Q 2 \ + A 21 ) is the fluorescence yield, t c i et is the period of detection, V is the 
volume probed by the laser source, fl is the solid angle over which detection occurs, and rj is the 
detection efficiency. The fluorescence yield, <t>, describes the fraction of de-excitation transitions 
that occur via spontaneous emission (fluorescence) relative to all de-excitation transitions (i.e. 
spontaneous emission and collisional quenching, assuming pre-dissociation and ionization are 
negligible). Neglecting the constants in Eq. 5.8, a generalized form for the fluorescence signal 
similar to that presented in Refs. 1 10 and 1 19, including its thermodynamic (x,P,T) dependencies 
and velocity dependence (arising from the Doppler effect), is: 

Slif K XsN T fB(T)B 12 IG(x s ,P,T, Xi> P ’T)t det (5.9) 

5.2 Non-Steady State Solution 

For a pulsed laser source, if N 2 is changing during the period of detection, then the entirety of 
Eq. 5.4 must be solved to obtain the time-dependent population, N 2 (t). Similarly, if the duration 
of the laser pulse is of the same order of magnitude as the time required to reach steady state, 
then the population at the end of the laser pulse, N 2 (t laser), can be obtained by solving Eq. 5.4. By 
assuming the laser intensity behaves as a Heaviside step function in time, the population N 2 (t), 
existing during laser excitation can be calculated by substituting the relation for Nj from Eq. 5.5 
into Eq. 5.4 and integrating with respect to time: 110 

N 2 (t) = Xsf B N T ^ ‘ (1 - e~ rt ), 0 <t< t laser (5.10) 

where r = W 12 + Q 21 + A 21 . The inverse of this value, r, is the characteristic time needed to 
achieve steady state. This solution assumes an initial condition which typically specifies the 
initial excited state population to be zero (N 2 (0) = 0). 

When the laser source is turned off, the only pathways to return to the ground state from the 
excited state are assumed to be through either spontaneous emission or collisional quenching. 
Therefore, for the period following laser excitation, the first term on the right-hand-side of Eq. 
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5.4 is zero. This modified version of Eq. 5.4 is then used to determine the excited-state 
population for the period after laser excitation by integrating with respect to time: 110 

N 2 (t) = Xsf B N T ^ ■ (1 - e~ rtlaser ) ■ e-Ct-tloserl/TL/F t > t laser (5.11) 


This solution assumes that the laser intensity is turned off instantaneously, with an initial 
condition given by Eq. 5.10, evaluated at t = t iaser . This solution shows that the population in 
the excited state decreases exponentially in time after the laser pulse. 

The term r LIF = (Q 21 + A 2 i) _ 1 in Eq. 5.11 is referred to as the fluorescence lifetime and 
describes the rate at which the population in a particular excited state transitions to a lower state. 
As with the derivation of Eq. 5.8, the total number of photons collected via fluorescence over the 
detection period is determined by multiplying Eqs. 5.10 and 5.11 by the spontaneous emission 
rate Einstein coefficient A 21 , accounting for the collection volume and detection system, and 
integrating with respect to time: 110 


$LIF ~ XsfB^T^B 12 IGti aser [ 1 


1 — e rt: laser 
r ^laser 



(l — e rt: laser ) 
tlaser/ T LIF 



(5.12) 


Like Eq. 5.8, this solution assumes that W 22 is small compared to A 2 i and Q 2 i. However, if 
the laser irradiance, /, is sufficiently large such that Wn and W 2 i are of the same order of 
magnitude as Q 2 i and A 22 , then both must be included in the solution. This gives r = W 12 + 
W 2 i + Q 21 + A 2 i during laser excitation. By defining the irradiance at which the fluorescence 
signal saturates as 5 ' 110 I sa t = (Q 21 + A 21 )/(B 12 + B 21 ) and rearranging r such that r — 
xlh [I / 1 sat + 1] , a more detailed formulation for the fluorescence signal, including effects from 
laser saturation, is given as: 110 
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(5.13) 


When I « I sat , Eq. 5.13 simplifies to a form similar to that of Eq. 5.10. 

5.3 Multi-Level Fluorescence Modeling 

In the preceding discussion, a simplified two-level model of fluorescence was put forth, 
which provides for an understanding of the most basic physical mechanisms and energy transfer 
processes involved. This basic description allows for the development of a general analytic 
relation between the fluorescence signal and these mechanisms and processes. However, such a 
simple description does not account for rotational or vibrational energy transfer occurring 
between the absorbing species and the surrounding gas mixture. Since each electronic energy 
level depicted in Fig. 5.1 can have vibrational and rotational fine structure for molecular species, 
such energy transfer processes result in a redistribution of the populations to various vibrational 
and rotational states. To account for these processes, analytical models including multiple energy 
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states and a system of population rate equations are used. Discussion of such modeling aspects 
can be found in Refs. 5, 109 and 1 10. 


5.4 Translational Temperature 

Translational temperature resulting from the random thermal motion of the absorbing species 
can be determined by measuring the Doppler-broadened width of the absorption line-shape 
function of Eq. 4.14. If the temperature, pressure, and mole fractions at a particular location are 
assumed to be constant, then the Boltzmann fraction and fluorescence yield in Eq. 5.9 at that 
location are also assumed to be constant. If these assumptions are valid, and if the laser line- 
shape function is known, then the fluorescence signal will be proportional to the overlap integral 
at a particular frequency. By scanning the laser’s spectral frequency across a particular transition, 
the absorption line-shape function can be obtained. 


Once the absorption line-shape is obtained, the contribution of homogeneous broadening 
mechanisms to the line-shape must be determined to allow for calculation of the Doppler- 
broadened line width. Homogeneous broadening can include components from Van der Waals 
broadening and Stark broadening effects, in addition to collisional pressure-broadening 
effects. In many cases, the Van der Waals and Stark effects can be ignored, resulting in 
homogeneous broadening from pressure alone. If the pressure is known, then the component of 
the line-shape due to homogeneous broadening can be determined. Alternately, if the pressure is 
sufficiently low, then collisional pressure broadening can be ignored. A fitting algorithm can 
then be used to determine the contribution of Doppler broadening to the line-shape. This allows 
for a more accurate temperature calculation as only one parameter, Av D , must be iterated upon to 
fit the line-shape. The translational temperature, T tmns , of the absorbing species can then be 
obtained by modifying Eq. 5.2: 


[ trans 


m 

8ln(2 )k B 



(5.14) 


120 122 123 127 

Both single-photon “ " " and two-photon “ " “ methods have been used to measure 
translational temperature. The requirement of scanning a laser’s frequency to resolve the 
Doppler-broadened line-shape can limit the time resolution of the temperature measurement such 
that only average translation temperatures can be obtained. 
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Figure 5.3 shows LIF data (dotted 
line, upper plot) from Ref. 120 obtained 
by performing a frequency scan with a 
semiconductor laser source over an argon 
metastable transition centered at 810.4 
nm. Each frequency scan over the 
absorption transition took 500 seconds to 
complete. In the experiment, argon 
plasma was generated using an 
inductively-coupled plasma torch. The 
translation temperature, T, and electron 
number density, n e , were obtained by 
fitting the Voigt line- shape function, 
which accounted for both Doppler- 
broadening and several collisional 
broadening mechanisms, to the 
experimental LIF data (solid line, upper 
plot). The lower plot in Fig. 5.3 shows 
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Figure 5.3. Experimental LIF data measured (dotted 
line, upper plot) in an Argon plasma and corresponding 
fit of Voigt line-shape (solid line, upper plot) to the 
experimental data (upper plot) and difference between 
experimental data and fit (lower plot). Image reprinted 
from Ref. 120 with permission of the authors and the 
publisher. 


the residual: the difference between the experimental LIF data and the fitted line-shape. An 
analysis of one of the scans in Ref. 120 gave T tra ns = 6800±690 K, resulting in a -10% 
measurement uncertainty. 


5.5 Rotational Temperature 


For molecules, the rotational temperature, T rot , can be obtained by exciting two absorption 
transitions and relating the signal intensities observed from each transition to the rotational 
temperature through modification of Eq. 5.9. This method was demonstrated in Ref. 128 using a 
two-photon excitation method to measure rotational temperature at a point in a cold (-300 K) 
turbulent flow. Lor each rotational transition, the population fraction, f B , is dependent on the 
rotational energy, Fj, rotational quantum number, J, and rotational temperature, T rot , of the 
absorbing species. This assumes that excitation occurs in the same vibrational state, v, such that 
the vibrational energy, G v , is constant for each probed rotational transition. If j s , G. <t>. and t det 
are assumed to be independent of the rotational state of the absorbing species, then the ratio, R, 
of fluorescence signals can be simplified to the following relation, similar to that in Refs. 1 10, 
129-132: 


S LIF,i _ g _ q B 12 ,i B i (2/i+l) ex p[ P J^/^B^rotl 15) 

S LIF,j B l 2,jP j (2//+l)exp[ , j / ^B^rot I 

where the subscripts i and j refer to the particular excited absorption transition, E is the laser 
energy such that E — at iaser I, a is the beam cross-sectional area, ti aser is the temporal pulse 
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width of the laser, and C is a constant which includes all terms that are independent of the 
rotational state of the absorbing species. Solving for T rot in Eq. 5.15 yields: 



Here, A E rot — ( Fj t — Fj ; ), and is the difference in energy between the probed rotational levels. 
To minimize the uncertainty associated with the measured rotational temperature, 8T rot , a 
propagation-of-error analysis can be performed on Eq. 5.15 by computing the derivative of the 
signal ratio, R , with respect to T rot : m,m 


8 T r pt 

Trot AE ro t R 


(5.17) 


This analysis suggests that by increasing A E rot , and thereby probing two widely separated 
rotational energy levels, the uncertainty in T rot can be proportionally decreased for a given error 
in /?. 130 


Figure 5.4 shows a mean rotational temperature map of nitric oxide measured in the wake of 
a 6.35-mm-thick flat plate in a Mach 3 turbulent flow, taken from Ref. 133. The map was 
obtained by averaging a series of images taken for excitation from rotational levels / = 8.5 and J 
= 10.5, computing the ratio of the two averaged images, and using a relation similar to that in Eq. 
5.16 to compute T rot . In this experiment, it was shown that the turbulent nature of the flowfield 
required excitation of two relatively closely spaced rotational levels (J = 8.5 and 10.5). This 
requirement arose from the nonlinear relation between R and T rot in Eq. 5.16, which can be 
heavily influenced by flowfield 

turbulence. 

Figure 5.4 shows a time-averaged 
temperature measurement obtained 
using a single laser, which probed the 
two transitions in separate wind tunnel 
runs. To obtain instantaneous 

temperature measurements in a plane 
using this method, a two-laser, two- 
camera system can be used, as in Refs. 

129-131. In Ref. 130, instantaneous Figure 5 . 4 . Rotational temperature map obtained in the 
two-line temperature imaging of OH in supersonic wake of a flat plate using two-line rotational 
a shock tube at a nominal pressure of thermometry. Image reprinted from Ref. 133 with 

40.53 kPa and temperature range of P ermission from the authors and S P rin S er Science and 

Business Media. 

1500-2950 K resulted in temperature 
measurements with errors of ~20%-25%. 
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Multi-rotational-line temperature methods also exist and can potentially provide relatively 
higher sensitivity and dynamic range T rot measurements, particularly at lower temperatures 
where rotational energy levels are closely spaced. 134 Such methods are appropriate when large 
variations in T rot are expected, requiring probing of several rotational levels. The rotational 
temperature is computed by exciting multiple rotational lines and measuring the signal. This 
measured fluorescence signal, together with the measured laser energy are then substituted into 

i in 

the following relation: 


In 


S LIF,i 

Pl2,i^i ( 2 //+ 1 ) 


(5.18) 


A plot of this logarithmic term versus F Jti , known as a Boltzmann plot, yields a linear trend with 
a slope of — ^B^rot- Refs. 134-137 used this method to characterize rotational temperatures in arc 
jet flows. An alternative multi-line approach used in Ref. 138 for measurements in a flame 
involved fitting a simulated rotational spectrum to experimental data to obtain T rot . 


5.6 Vibrational Temperature 

The vibrational temperature, T vib , can be measured using an approach similar to that used for 
rotational temperature. For T vib measurements, a two-line approach can be used. Such an 
approach is presented in Ref. 139, in which the same rotational levels, J t — Jj, are probed in two 
different vibrational levels, v t A Vj , resulting in nearly constant rotational energies, Fj, but 
different vibrational energies, G v . Assuming Xs , G • , and t det to be independent of the 

rotational and vibrational levels, the ratio of fluorescence signals, R , can be used to measure 
vibrational temperature, T vib , in a way analogous to that in Eqs. 5.15 and 5.16. 


As with rotational temperature measurements, a multi-vibrational-line temperature 
measurement method can also be used to infer T v u, by modifying the logarithmic term from Eq. 
5.18: 132 
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(5.19) 


If the rotational level, J t , is fixed and the rotational energy, Fj h is kept approximately constant, 

1 T? 

then the vibrational temperature can be computed independently of rotational temperature. 
This is done by generating a Boltzmann plot of vibrational energy, G v i , versus this logarithmic 
term, resulting in a linear relation with a slope equivalent to —1/ ( k B T vib ). 

Figure 5.5(a), taken from Ref. 132, shows a Boltzmann plot generated using a relation like 
that in Eq. 5.19 to infer vibrational temperature. The image data in the experiment was used to 
measure T V i b at a point downstream of a free-piston shock tunnel nozzle on the centerline of the 
flow. In this experiment, three absorption transitions were selected that had nearly constant 
rotational energies, Fj, but differing vibrational energies, G v . In Fig. 5(b), also taken from Ref. 
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132, the vibrational temperature of NO was observed to be constant with distance downstream of 
the nozzle, even though the nozzle was conical and the rotational temperature was shown to be 
decreasing with distance downstream, as the flow expanded. A simple computational model 

1 T? 

(shown in Fig. 5(b)) was not able to accurately predict the measured vibrational temperature. 

As shown in Fig. 5(b), T vi b = 785±30 K, giving an uncertainty of ~4%. 
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Figure 5.5. Vibrational temperature measurement obtained by probing multiple vibrational levels and 
generating (a) a Boltzmann plot and (b) vibrational temperatures plotted as a function of distance 
downstream. Images taken from Ref. 132 with permission of the authors. 

5.7 Species Concentration 

Returning to the simplified Eq. 5.8, it is possible to quantify the number density, j s IV r , using 
LIF. However, several factors complicate the interpretation of a LIF signal and its relation to 
concentration. For instance, Q 21 in Eq. 5.3 depends on the concentration and temperature of the 
absorbing species and all other constituent species in the probed mixture. If the temperature of 
the mixture were unknown, then a separate measurement of temperature would be required to 
determine both Q21 and the temperature-dependent Boltzmann fraction. The measurement system 
also requires an appropriate absolute intensity calibration to a known reference, which can be 
further complicated by the collection optics and filters used in such a measurement. The response 
of the detection system as a function of $LIF intensity and spectral frequency must also be well 
understood. The following discussion highlights some of the methods that have been used to 
provide concentration measurements that address some of these issues. 
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5.7.1 Non-quantitative flow visualization 

A relatively simple approach to 
obtaining a qualitative visualization of 
species concentration is to use laser- 
induced fluorescence to visualize the 
spatial distribution of the absorbing 
species. To first-order, the fluorescence 
signal in Eq. 5.9 is proportional to the 
relative species concentration, Xs^t , 
assuming all other terms are constant. In 
hypersonic nonequilibrium test facilities, 
the formation of reactive species at the 
high enthalpies required to simulate true 
flight conditions can influence the gas 
dynamics under study. For example, arc- 
heated facilities are often used to simulate 
true flight enthalpy conditions by 
depositing significant energy into a test 
gas and then expanding that gas to 
hypersonic Mach numbers. This process can result in a non-uniform species distribution 
resulting from the unsteady nature of the electrical arc. 

Figure 5.6 from Ref. 140 shows instantaneous nitric oxide (NO) planar FIF images in false 
color acquired in the Hypersonic Materials Environmental Test System (HyMETS) arc-heated 
facility at the NASA Fangley Research Center for a flow simulating earth’s atmosphere at 
several enthalpy conditions. The flow in this figure is from left to right. The object on the right- 
hand edge of the field of view is a 25 mm diameter silicon carbide tipped probe aligned with the 
axis of symmetry of the facility nozzle. For the lowest enthalpy condition shown (6.5 MJ/kg), the 
spatial distribution of NO within the core of the nozzle flow is irregular, suggesting that the 
energy deposition by the arc upstream of the nozzle is non-uniform. As the flow enthalpy 
increases to 18.4 MJ/kg, the distribution of NO remains irregular and the total production of NO 
appears to decrease. 

Other examples of such FIF-based flow visualization experiments are provided in Ref. 141. 
The interpretation of such flow visualization images and their relation to quantitative species 
concentration measurement is complicated. The remaining terms in Eq. 5.9, such as the 
quenching rate constant, often vary significantly in high-speed nonequilibrium flows and can be 
difficult to determine. The sections that follow detail approaches for quantitative measurements 
of concentration which address such complicating factors. 



6.5 MJ/kg 8.0 MJ/kg 9.5 MJ/kg 



11.3 MJ/kg 15.4 MJ/kg 18.4 MJ/kg 


Figure 5.6. Instantaneous nitric oxide planar LIF flow 
visualization images in false color acquired in the 
Hypersonic Materials Environmental Test System 
(HyMETS) arc -heated facility at NASA Langley 
Research Center. Flow is from left to right. Images taken 
from Ref. 140 with permission of authors. 
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5.7.2 Saturated LIF 


To avoid the complications associated with an unknown quenching rate constant, a relatively 
high laser irradiance where / » I sat can be used to probe the absorbing species. 5 In this instance, 
the stimulated absorption and emission rates dominate both collisional quenching and 
spontaneous emission rates; that is, (W 12 + W 21 ) » (Q 2 1 + A 21 ). The result of this leads to a 
modified form of Eq. 5.8 where stimulated emission is included and the fluorescence signal is of 
a form similar to that in Ref. 5: 

S l ,f = XsfBN T A 21 -^f-t det V±i 7 (5.20) 

# 12+^21 47r 

If A 2 i , B 12 , and B 2 i are known, then this method can be used to measure concentration if the 
remaining terms, such as f B , can be determined. In practice, it is difficult to achieve I » I sat 
because high powered pulsed laser beams are typically Gaussian, both spatially and temporally, 
so much of the acquired signal can originate from the lower intensity edges of the Gaussian 
beam, which excite the fluorescence linearly. 5 Further discussion of this method is found in Refs. 
5 and 109. 

5.7.3 Pre-dissociation LIF 

Pre-dissociative fluorescence involves inducing a transition to a level in the excited state, 
which can then couple to a dissociative state, resulting in dissociation of the molecule into 
smaller molecules or atoms. Pre-dissociative fluorescence concentration measurements are 
similar to saturated fluorescence measurements in that the dependence on Q 2 i in Eq. 5.8 can be 
removed, in this case by exciting the absorbing species to a pre-dissociative state where the pre- 
dissociation rate, Q pre , is faster than both Q 21 and A 21 . 5 A discussion of this mechanism can be 
found in Refs. 142 and 143 and a general discussion of the method can be found in Refs. 5 and 
109. An application of the pre-dissociative technique was presented in Ref. 144 to measure time- 
averaged OH concentrations in a supersonic hydrogen-air turbulent combusting flow. 
Concentrations as high as 5.4xl0 16 cm" 3 with uncertainties of 21% or less were reported. 144 

5.7.4 Short-Pulse LIF 

There are two methods of using short laser pulses to quantitatively determine concentration. 
The first probes the species of interest with a laser pulse having a duration much shorter than the 
time between collisions with other species. With this method, the concentration can be 
determined by relating the measured exponential fluorescence decay behavior to the initial signal 
magnitude during the short-pulse excitation period. While this decay is a function of both A 22 and 
Q 2 i, extrapolation of the exponential behavior to initial excitation provides an inferred peak 
intensity magnitude, which is assumed to be independent of collisional quenching effects. 109 This 
independence is assumed since a sufficient number of collisions, required to transfer energy non- 
radiatively from the excited state, would not yet have occurred. A discussion of this method can 
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be found in Ref. 109. Reference 145 uses a form of this short-pulse method for OH concentration 
measurements in a turbulent flame. 

A second method for measuring concentration involves using a laser pulse that is shorter than 
or comparable to the fluorescence lifetime but not shorter than the collisional timescales 
described in Ref. 5. In this approach, the fluorescence lifetime is measured directly, usually using 
a photomultiplier tube. This short-pulse LIF method can also be used to determine spontaneous 
emission 146 and collisional quenching 147 ' 148 rate constants when the thermodynamic conditions of 
the probed mixture are known. Even for unknown conditions, measurement of the fluorescence 
lifetime allows for the quenching rate to be determined and accounted for. Refs. 149 and 150 use 
a two-photon LIF technique in which fluorescence lifetime measurements are used to correct for 
collisional quenching effects in a similar manner. Figure 5.7(a), taken from Ref. 149, shows 
spectral scans of atomic nitrogen in the NASA Ames Aerodynamic Heating Facility arc jet flow 
(red data) and laboratory flow reactor (green data). Fluorescence measurements from the flow 
reactor and a krypton reference cell were used for intensity calibration in order to provide 
absolute atomic nitrogen number densities. The fluorescence lifetimes observed in the arc jet and 
flow reactor were used to correct for collisional quenching effects. Figure 5.7(b), also taken from 
Ref. 149, shows radial atomic nitrogen number density measurements for two air/argon arc jet 
runs taken 35.6 cm downstream of the nozzle exit. The reported uncertainty in number density 
was -12%. 149 




Figure 5.7. (a) Spectral scans of atomic nitrogen in the NASA Ames Aerodynamic Heating Facility arc 
jet flow and in a laboratory flow reactor and (b) measured number densities spanning the radial direction 
of the arc jet 35.6 cm downstream of the nozzle exit. Image taken from Ref. 149 with permission of the 
authors. 


A similar approach is being developed for mole fraction measurements in hypersonic 
turbulent boundary layers using naphthalene PLIF. 151452 One benefit of using naphthalene is that 
it sublimates at slightly elevated temperatures (with respect to room temperature), allowing for 
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the study of scalar transport effects in 
transitional boundary layers, turbulent 
boundary layers, and ablating surfaces. 

5.7.5 Bi-Directional Beam LIF 

If two overlapping, counter- 
propagating beams are tuned to the same 
transition, the ratio of fluorescence signals 
results in the cancellation of all terms in 
Eq. 5.9 at a point in the flow, with the 
exception of the spatially dependent laser 
irradiance. Thus, the spatially varying 
ratio of fluorescence signals can be 
equated to the spatially varying ratio of 
irradiances of the respective beams. The 
Beer-Lambert law of Eq. 4.3 can then be 
used to relate irradiance, and therefore 
fluorescence signals, to number density 

co 

via the relation: 

KM = N s (p s (v) (5.21) 

where <p s (v ) is the transition cross-section 
of the absorbing species, describing the 
absorbing species’ interaction with the 
incident irradiance on a per atom or 
molecule basis. 111 Combining Eq. 5.21 
with Eq. 4.4 allows for a relation between 
the natural logarithm of the signal ratio 
and the absorbing species number density 
to be made: 153 
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Figure 5.8. Graphical depiction of methodology used 
to measure concentration using bi-directional beam 
LIF. Top plot shows spatially varying fluorescence 
signals, middle plot shows natural logarithm of ratio of 
these signals, and bottom plot shows OH concentration 
calculated using Eq. 5.22. Image reprinted from Ref. 
153 with permission of the authors and Springer 
Science and Business Media. 


where ^ LIF, i (x) and ^ LIF, 2 (x) are the spatially varying fluorescence signal intensities of the 
forward- and backward-propagating beams, respectively. In Eq. 5.22, (p sQ is the peak transition 
cross-section, with (p s0 — <p s ( v o)- U1 


This concentration measurement technique was first demonstrated in Ref. 153. The technique 
only requires that a value for cp s 0 be known in order to make an absolute concentration 

i « 

measurement of N s . ' For total absolute concentration, however, f B must either be assumed or 
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measured, as from Eq. 5.5, N s — Xsfu^r- Figure 5.8, taken from Ref. 153, graphically depicts 
the process of obtaining a concentration of OH in a methane-air flame. In the top plot, the 
fluorescence intensities are plotted versus position across the flame. The middle plot shows the 
natural logarithm of the ratio of these two signals with position and the bottom plot shows the 
calculated absolute OH concentration using Eq. 5.22. Measurements of OH concentration were 
also made in a plane in Ref. 153 by imaging fluorescence from two counter-propagating laser 
sheets. Equation 5.22 was applied using each row of pixels to compute signal ratios. Single-shot 
imaging of a hydrogen/air/oxygen flame allowed for calculation of means and standard 
deviations of OH concentration. Using 10 images for each laser sheet direction (20 total images), 
Figure 5.8 from Ref. 153 showed a mean peak OH concentration of ~8.8xl0 14 cm 1 with a 
standard deviation of ~±1.7xl0 14 cm' 1 , which is -19.7% of the mean. 

5.7.6 Combined LIF/Rayleigh/Raman 

In reacting flows, flows undergoing transition-to-turbulence, and turbulent flows, the time- 
varying nature of flow properties makes it difficult to ascertain concentration. By making 
simultaneous time-resolved measurements of temperature and major species concentrations, the 
thermodynamic-dependent parameters in Eq. 5.9 (Boltzmann fraction, overlap integral, and 
fluorescence yield) can be determined, and a direct relation between concentration and 
fluorescence signal established. A combined fluorescence, Rayleigh scattering, and Raman 
scattering measurement technique was used for this purpose in Ref. 154. In that experiment, the 
Raman scattering was used to measure instantaneous concentrations for major species (CB, No, 
HoO, and H 2 ) in a turbulent hydrogen jet flame. This allowed for the mole fractions of the 
collision partners, Xh i n Eq. 5.3 to be calculated. The combined Rayleigh/Raman scattering 
measurements were then used to measure temperature. The temperature-dependent collisional 
cross-sections and relative velocity in Eq. 5.3 could then be calculated, and the collisional 
quenching rate, Q 2 i, determined. The temperature-dependent absorption line-shape and 
Boltzmann fraction could also be determined from the temperature measurement. The additional 
measurement of laser energy then provided for a direct relation between NO fluorescence signal 
and concentration to be established. In Ref. 155, this technique was used to make concentration 
measurements of CO in turbulent premixed and stratified CH 4 /air flames. A stated CO 
concentration accuracy and precision (single standard deviation) of 10% and 4.5%, respectively, 

155 

were given. 

5.8 Pressure 

Measurements of pressure can also be obtained using a LIF-based approach, although such 
measurements are not as common as LIF-based temperature and species concentration 
measurements (described in the preceding sections) for determination of the thermodynamic 
state of a gas. This is in part due to the implicit relationship between fluorescence signal and 
pressure in Eq. 5.9, as opposed to the explicit relationships with temperature and concentration. 
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One method by which pressure can be measured is by tuning to the center frequency of an 
absorption transition and relating the measured fluorescence signal to pressure. For such 
measurements, the selected absorption transition must be carefully chosen such that the 
temperature dependence (i.e. f B ) can be removed over the range of temperatures expected for the 
experiment. When this is done, the fluorescence signal is inversely proportional to the static 
pressure if certain additional simplifying assumptions are made. This approach has been used in 
low-quenching ( Q 21 « A 21 ) flows 133 or quenching-dominated (Q n » A 21 ) flows where the 
pressure broadening is much larger than the Doppler broadening. 156 ' 157 

Often for quenching-dominated LIF-based pressure measurements, however, the dependence 
of fluorescence signal on pressure when tuned to the center of an absorption transition is nearly 
removed, as N T (and hence pressure) can be factored out of both the numerator and denominator 
in Eq. 5.9. An alternative approach is to infer pressure from the absorption lineshape function. 
This can be accomplished if a separate measure of temperature can be made so that the 
absorption lineshape dependence on temperature can be isolated from its dependence on 
pressure. Common methods include those which relate the width of the absorption lineshape to 
pressure, 158 those which resolve the entire absorption lineshape function, 159 and those which 
relate the ratio of fluorescence signal from the wing of the lineshape to the signal from the center 
of the lineshape. 160 " 162 

5.9 Doppler-Based Velocimetry 

The Doppler effect can be used to determine flow velocities for atomic and molecular 
species. The translational motion of the absorbing species in the direction of the excitation 
laser’s propagation, described by a velocity component U, results in a shift of the absorption 

co 

line-shape function away from its transition center frequency, v 0 , according tor 

Av = % 0 (5.23) 

This velocity-dependent frequency shift of the absorption line-shape function is implicit in 
the overlap integral, G, and thus, its effect on fluorescence signal can be seen through the 
dependence of G on the velocity component, U, in Eq. 5.9. When a component of translational 
motion of the absorbing species opposes the laser’s direction of propagation, corresponding to 
the negative (— ) solution of Eq. 5.23, the incident laser radiation appears to be at a higher 
frequency from the perspective of the gas. Hence, as the laser’s frequency is scanned over the 
absorption transition, Y v , the measured intensity of fluorescence corresponding to this profile is 
shifted toward a lower frequency, or red-shifted. The converse is true when motion is in the same 
direction as the laser’s propagation, corresponding to the positive (+) solution of Eq. 5.23, where 
the absorption profile is shifted toward a higher frequency, or blue-shifted. 
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In one implementation, velocity component measurements can be made by scanning the laser 
over an absorption transition in both the measurement volume and a reference cell. In this case, 
the Doppler shift between the absorption profiles is used to compute an average velocity 
according to Eq. 5.23. Examples of such measurements include those in an arcjet, ’ ’ ' 

supersonic underexpanded jets, 158 ’ 159 ’ 164 ' 166 shock tunnel, 167 168 and non-reacting supersonic flow 
with a rearward-facing step. 169171 Such measurements require that the flowfield be relatively 
steady since shot-to-shot fluctuations in fluorescence intensity from thermodynamic and fluid 
dynamic variations can affect the measured line-shape. The velocity measurements can also be 
complicated by a frequency shift in the line-shape function resulting from collisional effects 
(pressure). In compressible flows, the pressure field can vary significantly, and hence the 
collisional shift in the line-shape function can likewise vary. Absorption of laser energy can also 
result in an apparent frequency shift in the line-shape function. Recall that the line-shape 
function must be inferred from the overlap integral, G, which represents the convolution of the 
absorption line-shape and laser spectral line-shape. If absorption is significant, then irradiance 
will vary spatially according to Eq. 4.3, resulting in a spatially varying overlap integral. An 
analysis of how absorption affects the line-shape frequency shift is presented in Ref. 168. If flow 
symmetry is assumed, then the frequency shift due to pressure can be estimated. 167 Flow 
symmetry can also be used in the application of Eq. 4.3 to correct for the frequency shift 
resulting from absorption. “ Alternatively, if two counter-propagating laser beams are used, 
these frequency shift effects can be completely removed. The use of counter-propagating beams 
results in two excitation peaks, separated in frequency by twice the velocity-induced Doppler 
shift. The need for a reference cell measurement to ascertain velocity can also be removed, as a 

165 172 173 

counter-propagating beam approach is self-referencing. ’ ’ 
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Figure 5.9 in Ref. 168 shows a two-component velocity measurement about a heat shield 
model taken in a hypersonic free-piston shock tunnel. In this experiment, laser sheets were 

directed in both the radial (vertical) and 
axial (horizontal) directions so that 
measurements of the Doppler-shifted 
absorption profiles for the respective 
directions could be obtained and 
compared with measurements from a 
static reference cell to infer velocities. 
An estimation of velocity errors incurred 
from the frequency shifts due to 
absorption and collisional effects was 
also performed. In Ref. 168, three 
separate absorption transitions were 
probed to measure two different velocity 
components, with one transition being 
used for both components. This resulted 
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Figure 5.9. Two-component 
using fluorescence-based 
technique. Vectors denote flow direction and color map 
represents magnitude. Image reprinted from Ref. 168 
with permission of the authors and Springer Science and 
Business Media. 
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in a measured axial freestream velocity of 2394+68 m/s and a measured radial velocity of 53+50 
m/s, giving respective uncertainties of -2.8% and -94.3%. 168 

Another form of the fluorescence-based Doppler velocity measurement is a fixed frequency 
method, which can allow for an instantaneous velocity component measurement. With the fixed 
frequency technique, a narrow linewidth laser is tuned off the absorption profile peak to a point 
where the slope of the profile is maximum, as described in Refs. 160 and 174. Assuming that the 
absorption profile is approximately linear in the region of maximum slope, the measured signal 
intensity can be related to the Doppler shift of the profile. This fixed frequency Doppler 
velocimetry technique has been applied to a free jet, 175 supersonic underexpanded jet, 176 and 
reacting supersonic flow. The stated random and systematic errors in Ref. 176, when added in 
quadrature, gave a total uncertainty of -12%. In Ref. 177, the stated lowest time-averaged and 
single-shot uncertainties achieved were -3% and -15%, respectively, for a 1600 m/s velocity 
range. 

5.10 Flow-Tagging Velocimetry 

Another technique by which velocity can be measured using fluorescence is flow-tagging 
velocimetry. Fluorescence-based flow-tagging velocimetry is a time-of-flight technique that 
involves laser excitation — or tagging — of the gas along a line, series of lines, or grid pattern. 
With this form of velocimetry, the species of interest in the gas absorbs the incident radiation 
from a laser source, which induces either of the following: 1) fluorescence, 2) a reaction that 
forms a product that then emits a photon via fluorescence, or 3) a reaction that forms a product 
that can then be probed with another laser source to induce fluorescence. Images of the 
fluorescence pattern are acquired at two time delays, with velocity computed by measuring the 
displacement of the tagged molecules between images. Typically, a line or series of lines can be 
used to measure a single-component of velocity while a crossed grid pattern can be used to 
measure two-components. Two general fluorescence-based methods of flow-tagging velocimetry 
are discussed here; one that requires a single laser source and another that requires multiple laser 
sources. The main advantage of flow-tagging velocimetry, as compared to most Doppler-based 
methods (which are time averaged), is that it can make instantaneous (single-shot) measurements 
with fast time resolutions (as short as a few hundred nanoseconds). A disadvantage, however, is 
that flow-tagging velocimetry cannot provide full velocity field information. A broader 
discussion of molecular-tagging velocimetry, which relies on molecular tracers for flow-tagging, 
is provided in Ref. 178. 
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5.10.1 Single-laser methods 

The first method involves either direct or indirect excitation of fluorescence with a single 
laser source. The first application of this method to a gaseous flow involved excitation of 
phosphorescence of biacetyl molecules, as described in Ref. 179. If this method is used, the 
fluorescence lifetime of the tagged molecules must be long enough so that advection provides for 
measurable displacements with tagged regions having signal intensities that are still above the 
detection limit of the imaging system at the time the second image is acquired. Typical 
experiments involve capturing a reference image acquired during the tagging process, or a 
relatively short time thereafter. If a single-framing camera is used, such as in Ref. 180, a single 
reference image or set of reference images is acquired. The timing of the single-frame camera is 
then delayed and a subsequent image or set of images is then acquired. If a two-camera system or 
dual-framing camera is used, such as in Ref. 181, the delayed image is acquired in sequence after 
the reference image. The velocity is computed by measuring the displacement of the tagged 
molecules that occurs in the time between when the reference and delayed images were acquired. 
The form of this technique relying on direct excitation of fluorescence for flow tagging has been 
applied to the study of supersonic jets, ' hypersonic boundary layers, ’ ’ ’ and arcjet 

flowfields. 122 

Figure 5.10, taken from Ref. 180, shows images from tagging a single line of nitric oxide 
using direct excitation of fluorescence within a flat plate hypersonic laminar boundary layer. In 
this Figure, the left-most image corresponds to the reference image, while the remaining images, 
from left to right, correspond to delayed images taken at 250 ns, 500 ns, and 750 ns after tagging, 
respectively. Measurements of freestream velocity, spatially averaged from a point just above the 
velocity boundary layer (3 mm) to 15 mm above the flat plate, resulted in a mean of 3,035±100 

i on 

m/s at 90% confidence, giving an uncertainty of 3.3% of the mean. Single-shot uncertainty 
estimates for a 3,000 m/s freestream flow 
and for camera delay settings of 250 ns, 

500 ns, and 750 ns were 4.6%, 3.5%, and 
3.5%, respectively. 180 

An indirect excitation scheme, as 
described in Refs. 191 - 193 , relies on 
photodissociation of molecular nitrogen 
for flow tagging. The technique uses a 
femtosecond laser pulse to dissociate 
molecular nitrogen into two nitrogen 
atoms, which then recombine after a 
collision, forming molecular nitrogen in 
an intermediate state. A subsequent 
collision brings the molecular nitrogen to 



Figure 5.10. Single-line excitation of nitric oxide 
fluorescence used to study hypersonic boundary layer 
flow over a flat plate. Images, from left to right, 
correspond to camera delay settings of 0 ns, 250 ns, 
500 ns, and 750 ns. Image reprinted from Ref. 180 
with permission of the authors. 
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an excited electronic B state, which then emits a photon via fluorescence upon transitioning to 
the excited electronic A state. Ref. 193 provides a description of this process. One benefit of this 
indirect technique, known as Femtosecond Laser Electronic Excitation Tagging (FLEET), is that 
the recombination rate of dissociated atomic nitrogen allows for a much longer fluorescence 
lifetime. This would allow displacements to be measured over greater time scales, providing for 
accurate measurements of velocity in low-speed flow regions, such as in a hypersonic wake flow. 
Additionally, the technique relies on molecular nitrogen for tagging, which is present in most 
hypersonic facilities. 

5.10.2 Multi-laser methods 

A second method of fluorescence-based flow-tagging velocimetry involves writing a line, 
series of lines, or grid pattern into the flowfield by one of several laser-based mechanisms. This 
pattern can then be interrogated, or read, by subsequent laser pulses to induce fluorescence, 
allowing for the determination of velocity through measurement of the displacement of the 
pattern. Such techniques usually involve two or three different lasers and are therefore more time 
consuming to set up and more difficult to execute. 

One mechanism by which a pattern can be written into the flowfield is via ionization of the 
absorbing species, known as Laser Enhanced Ionization (LEI) flow tagging. The tagging process 
is accomplished by promoting the species (such as sodium in Refs. 194-196) to a higher energy 
state via laser excitation near the ionization limit. Collisions then result in the ionization of the 
species, with the tagging pattern corresponding to the ionized regions. Subsequent laser pulses 
are used to induce fluorescence of the absorbing species in regions that have not been photo- 
ionized. Supersonic measurements of velocity in a shock tube were performed using LEI flow- 
tagging in Refs. 194 and 195, and hypersonic velocity measurements in an expansion tube were 
performed in Ref. 196. 

A second mechanism that can be used to write a pattern into the flowfield is by vibrational 
excitation of molecular oxygen via Raman pumping. This is followed by reading the pattern of 
the vibrationally-excited oxygen by inducing fluorescence. The technique, known as Raman 
Excitation and Laser-Induced Electronic Fluorescence (RELIEF), is advantageous as it relies 
on the flow-tagging of oxygen which is a common working gas of most hypersonic facilities. 
This technique has been used to characterize turbulence in a free jet 198 and underexpanded jet. 199 
The RELIEF technique, however, is limited to temperatures below 750 K. Above this 
temperature, a significant fraction of oxygen molecules are vibrationally excited, making it 
difficult to distinguish the tagged molecules from the background. 200 

Yet another mechanism involves using one laser to photo-dissociate a molecular species. 
This results in the formation of a product species for which a second laser can be used to read the 
location of the written pattern by exciting laser-induced fluorescence in the product species. The 
formation of the product species typically occurs through one or more reactions. A list of partner 
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species used in the writing and reading process include: H 2 0-0H , 201-207 N 2 0-NO , 208 O 9 -O 3, 209,210 
N 2 /0 2 -N0 , 21 1-214 and N0 2 -N0 . 184,215-222 

Figures 10(a) and 10(b), taken from Ref. 219, show images of vibrationally excited NO 
fluorescence in a grid pattern formed via photodissociation of N0 2 in a mixture containing 6.3% 
N0 2 in N 2 . The images were taken in a supersonic underexpanded jet 400 ns (Fig. 5.11(a)) and 
800 ns (Fig. 5.11b) after the pattern was written into the flow with a two-dimensional array of 
355 nm beams. Two components of velocity were obtained by relating the displacement of the 
grid in the left image to a grid imaged in a stationary gas. The upper half of Fig. 5.11(c), taken 
from Ref. 219, shows measured streamwise velocities compared with computation, shown in the 
lower half of Fig. 5.11(c). The use of two pulsed dye lasers in this experiment permitted the 
excitation of both a low (7 = 1.5) and high (7 = 8.5) rotational level within the same vibrationally 
excited (v = 1) NO state. The fluorescence images corresponding to the low-7 (Fig. 5.11(a)) and 
high-7 (Fig. 5.11(b)) rotational levels also allowed for measurement of the rotational temperature 
by using a calculation similar to that presented in Eq. 5.16. The top half of Fig. 5.11(d), from 
Ref. 219, shows measured rotational temperature compared with computation, shown in the 
lower half of Fig. 5.11(d). Stated root-mean- square (RMS) uncertainties in the velocity 
measurement were ~5% with high signal-to-noise.' The stated RMS uncertainties in rotational 
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temperature ranged from 9% to 35% prior to the Mach disk. 

N0 2 -NO flow tagging velocimetry has also been used to study boundary layer transition on a 
flat plate in a Mach 10 wind tunnel. - " ’ ~ Pure N0 2 was seeded into the boundary layer through a 
spanwise slot located downstream of the sharp leading edge. The angle of attack of the flat plate 
was 20 degrees, reducing the edge Mach number to about 4.2. Parallel focused beams of 355 nm 
light oriented normal to the surface of the flat plate (and in a plane parallel with the streamwise 
direction of the flow) dissociated the N0 2 , creating NO. After a 40 ns delay, the NO was probed 
by a 226 nm laser, at which time the camera acquired an image. One microsecond later, a second 
226 nm laser probed the NO and the second image was acquired with the same camera. The two 
images were compared to compute velocity profiles. Figure 5.12 shows the resulting velocity 
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Figure 5.11. Flow-tagging images of vibrationally excited NO fluorescence obtained (a) 400 ns and (b) 
800 ns after photodissociation of NCF. The two-dimensional grid pattern allows for calculation of two 
velocity components. Measured (c, top) and computed (c, bottom) streamwise velocity maps. Measured 
(d, top) and computed (d, bottom) rotational temperature map. Image reprinted from Ref. 219 with 
permission of the authors and the publisher. 
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profiles for two cases, one with no trip (i.e. tripping element, a protuberance designed to trip the 
flow from laminar to turbulent) and one with a 1-mm tall cylindrical trip, where the 
measurements were made downstream of, and on the centerline of, the trip. (The laminar 
boundary layer thickness was also approximately 1 mm thick.) The figure shows mean profiles 
(black, top two charts) as well as profiles of the fluctuating component of velocity, u’ (red, 
bottom two charts). For the case of no trip, the velocity profiles compare well with a 
compressible Blasius solution (shown in light grey). When the trip is present, the measured mean 
profiles depart from the laminar solution, showing a profile that is more full than the laminar 
profile near the plate surface with a pronounced velocity deficit near the edge of the profile. The 
fluctuating streamwise velocity increases by a factor of three, up to 250 m/s, between the two 
cases, with the fluctuations highest in the wake of the trip. The single-shot measurement 

precision was 15-25 m/s, which was 1-2% of the maximum velocity in the boundary layer. The 

221 

accuracy was estimated to be 5-15 m/s. 

5.11 Advantages and Limitations of Laser-induced Fluorescence 

The time scales associated with inducing fluorescence via laser excitation are typically a few 
hundred nanoseconds, which is much shorter than hypersonic flow time scales, therefore 
providing sufficient temporal resolution for high-speed transition-to-turbulence measurements. 
The availability of pulsed lasers capable of picosecond and femtosecond excitation allow for 
fluorescence measurements with time scales much less than those associated with collision and 
reaction time scales. Recently developed kHz- and MHz-rate pulsed laser systems have allowed 
image sequences consisting of tens to thousands 223 of images to be acquired, providing time- 
resolved information pertaining to high-speed fluid dynamic behavior. Both fluorescence-based 
velocimetry and visualization ' “ experiments in hypersonic flow fields have been 
performed with these types of laser systems. The spatial resolution of a LIF technique is also 
sufficient for many applications, and higher than other techniques such as Raman or CARS, with 
laser sheet thicknesses typically in the range of 0.1 to 1 mm and magnifications of tens of pixels 
per mm, depending on the experimental setup. Fluorescence-based measurements are more 
sensitive than other techniques (for example Raman spectroscopy) with sensitivity on the order 
of parts-per-million or better. 5 
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Figure 5.12. Streamwise velocity profiles on a flat plate in a Mach 10 facility for the case of no trip (a) 
and a k = 1 mm tall, 4 mm diameter cylinder trip (b) and the fluctuating streamwise velocity for no trip 
(c) and a k = 1 mm tall, 4 mm diameter cylinder trip (d). Images reprinted from Ref. 222 with permission 
of the authors. 


Additionally, a wide range of species including intermediate combustion species can be 
probed using fluorescence techniques. Reference 109 provides an extensive listing of many 
species that have been detected using LIF and other methods. Another advantage of LIF is that it 
is readily extended to planar or volumetric measurement (see Refs. 226 and 227). 


Several factors complicate acquisition and interpretation of LIF signals, complicating 
quantitative measurements. Quenching of the fluorescence, which prevents easy quantification of 
signal intensities, was discussed extensively above. Absorption of laser energy as the laser light 
passes through the flowfield can limit the effectiveness of fluorescence-based measurements, as 
the energy will decrease in an exponential manner over a given spatial path length according to 


66 


the Beer-Lambert relations in Eqs. 4.3 and 4.4. This makes quantitative measurements difficult, 
as the laser energy at a particular location may not be easily determined. Absorption can be 
significant when the concentration of the absorbing ground state population is high, the transition 
cross-section is relatively large, the Einstein B coefficient for stimulated absorption is relatively 
large, and/or the path length through which the laser radiation passes is relatively long. To avoid 
strong absorption, a transition may be selected for which the population is small based on 
analysis of the Boltzmann fraction, as was done for PLIF visualization measurements in a 
hypersonic shock tunnel described in Ref. 228 . Absorption can also limit measurement 
capabilities when fluorescence from the probed volume is re-absorbed by the species of interest. 
This effect, known as radiative trapping , occurs when fluorescence emission at frequencies 
readily absorbed by highly populated states must pass through gas containing these potential 
absorbers before reaching the imaging system. Measurements in a non-uniform or turbulent 
mixture are especially susceptible to errors associated with absorption and radiative trapping 
effects, as the absorption coefficient is a spatially- and temporally-varying property. A discussion 
of these issues, and some methods used to circumvent them, is provided in Ref. 5. 

Another disadvantage of LIF is that it usually probes only a single species at a time, 
compared to Raman or CARS, which can interrogate many species simultaneously. However, the 
fact that PLIF can measure spatial distributions of species can compensate for the single- species 
capability in some applications. 

Consideration must also be given to the optical access of the test facility. Since much of the 
work described in this discussion requires fluorescence excitation using laser frequencies in the 
UV portion of the electromagnetic spectrum, the window material used in hypersonic facilities 
must be capable of transmitting such frequencies with minimal absorption. Also, typically two or 
three windows are required for LIF or PLIF applications. The laser (beam or sheet) is typically 
brought in from one window and observed through another window at right angles to the first 
window. A third window can allow the laser to leave the test section, reducing scattered light, 
and allowing the quantification of absorption of the laser beam/sheet in some applications. These 
windows are typically relatively large compared to those required for CARS or diode laser 
absorption measurements. When short-pulse LIF experiments are performed, the inverse 
relationship between the pulse temporal width and spectral width (which for a Gaussian pulse is 
r ia Se r = 0.44/Av) 111 can result in distortion of the temporal characteristics of the pulse. This is 
because the refractive properties of the optical windows of the test facility affect the speed with 
which the frequency components of the short pulse pass through the window material. Reference 
111 discusses issues related to short pulse excitation and interaction with optical components. 
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6. Rayleigh and Raman Scattering 

6.1 Introduction 

When a light beam passes through a 
gaseous medium, it can interact with the 
gas molecules or particles in the gas, 
thereby scattering light away from the path 
of the incident beam. Elastic scattering 
occurs if no energy is gained or lost to the 
medium. If energy is either absorbed or 
lost by the medium, the scattering is 
inelastic. Light scatter from particles that have a diameter, d, on the order of or larger than the 
light wavelength, A, is termed Mie scattering. Scattering for which d« A is known as 
spontaneous Rayleigh scattering if elastic and spontaneous Raman scattering if inelastic. These 
processes are shown schematically in Figs. 6.1 and 6.2. Represented on the energy level diagram 
in Fig. 6.1, a photon excites a molecule from an originating state, 1, to a ‘virtual’ state, 2, from 
which the scattered photon is emitted. The virtual state, represented by a dashed line, is not an 
actual resonant state of the molecule. Instead, it indicates a non-resonant, short-lived state in 
which the electron distribution of the molecule is distorted. This virtual state immediately 
relaxes to the originating state (in the case of Rayleigh scattering) or another state (in the case of 
Raman scattering). Relaxation to a higher lying (e.g. vibrational) state than the originating state 
is termed Stokes Raman scattering. In this case, the molecule absorbs a quantum of energy 
through this process. In anti-Stokes Raman scattering the molecule imparts a quantum of energy 
to the scattered photon so that the scattered photon has higher energy than the incident photon. 
In this case, the originating state must not have been a ground state. 

Figure 6.2 shows notional Raman/Rayleigh 
spectra, not drawn to scale. Rayleigh scattering 
is shown at the laser’s wavelength. Discrete pure 
rotational Raman lines, associated with rotational 
quanta imparted to or subtracted from the 
incident laser frequency, are shown on opposite 
sides of the Rayleigh peak. Vibrational Raman 
bands are located further away, spectrally shifted 
towards the red (Stokes) and the blue (anti-Stokes). The vibrational Raman bands show 
rotational fine structure. Raman scattering is much weaker than Rayleigh scattering. Raman 
scattering is typically three orders of magnitude smaller than Rayleigh scattering for most gases 
of interest to supersonic and hypersonic flows. 5 
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Figure 6.2. Notional Raman/Rayleigh spectra. 
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Figure 6.1. Energy level diagram indicating incident 
radiation, Rayleigh scattering and Raman scattering. 
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6.2 Theory of spontaneous Rayleigh and Raman scattering 

Comprehensive reviews of the theory and application of Rayleigh and Raman scattering have 
been given by others. ’ ' ’ Herein we provide a brief introduction and overview. Later 

sections will show how this theory is applied to measure thermodynamic properties. Consider an 
electromagnetic wave incident upon a molecule, perturbing the molecule’s electron cloud and 
making it oscillate at the same frequency as the incident wave. These oscillations cause a 
periodic charge separation within the molecule, known as an induced dipole moment. Oscillating 
dipole moments act like antennas, emitting radiation. As described in more detail by Baldwin,” 
if these oscillations are in-phase, the emission adds constructively, producing a coherent beam. 
If out-of-phase (for example, in the direction orthogonal to the beam), the emitted light interferes 
destructively and the radiation cancels. For a monochromatic plane wave passing through a gas, 
the constructive interference occurs only in the forward direction. The resulting emitted coherent 
light is perfectly in-phase with and indistinguishable from the incident wave. In a gas composed 
of a finite number of molecules, the destructive interference at other angles is not fully complete 
because of statistical variations in the number of particles located in different wavelength- sized 
volumes of the gas. 23 T That is, there are not exactly the same number of particles in each 

wavelength sized volume, which would be required to cancel out the radiation perfectly. This 
statistical variation in the number density then leads to Rayleigh and Raman scattering. Rayleigh 
and Raman increase significantly in intensity at shorter wavelengths: both scale approximately as 
1/A, 4 . 5 

The radiant intensity, 1°, which is the scattered power per unit solid angle, is proportional to 
the square of the induced dipole moment. The induced dipole moment, p, is given by: 5 

p = s 0 a E (6.1) 

where s 0 is the permittivity of free space, a is the molecular polarizability and E is the incident 
electric field given by E — E 0 cos (m 0 t), where E 0 is the amplitude of the electric field, co 0 is the 
frequency of the laser light and t is time. 

The polarizability of a molecule depends on its internal structure and varies with time during 
vibrational oscillations at the natural frequency of the molecule, a) v , vibrating in the direction of 
its normal spatial coordinate, Q. The polarizability can be approximated with a Taylor series 
expansion: 

a = a ° + i^) 0 dQ (6.2) 

where the small physical displacement, clQ, of the atoms about their equilibrium positions 
(denoted by the subscript 0) during vibrations is: 

dQ = Q 0 cos (co v t) (6.3) 
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Combining Eqs. (6. l)-(6.3): 


[ a ° + © 0 cos (<°vt) Zo% co s(<u 0 t) (6.4) 

Expanding and using a trigonometric identity: 5 

V = a 0 s 0 E^cos(m 0 t) +|(0) s 0 Q 0 £^[cos(m 0 - <u„)t + cos(m 0 + m v )t] (6.5) 


The first term on the right side of Eq. (6.5) describes Rayleigh scattering at a frequency 
corresponding to the incident laser’s wavelength. The second term indicates Raman scattering 
which is shifted from the Rayleigh scattering by ±co v resulting in upshifted (anti-Stokes) and 
downshifted (Stokes) Raman scattering. While significant theory has been developed to describe 
the physics of Raman and Rayleigh scattering, the strengths of the scattering for different gases 
is generally measured and reported as a temperature-independent differential cross section: 5 



which can be rearranged as: 

< 6 - 7 > 


where the subscript zz refers to a polarization in the z (vertical) direction caused by an incident 
electric field oriented in the z direction, N is the number density of the gas and I is the laser 
irradiance. 5 The differential cross sections vary by process (Rayleigh vs. Raman), by molecule, 
and vary with the laser wavelength but are independent of pressure and temperature. 


Shift, V 


Rayleigh scattering from different 
molecules cannot be distinguished 
spectrally, so it is not usually used to 
detect individual species. Under 
circumstances where the composition is 
fixed or known, or limited to vary under 
controlled conditions, ' Rayleigh 
scattering can be used to measure the gas 
density, p. As discussed further below, the 
Rayleigh scattering cross section varies 
from molecule to molecule, with some 
hydrocarbon species having cross sections 
more than an order of magnitude larger 

than N 2 . The cross section for a mixture of gases is equal to the mole-fraction-weighted average 
of the individual cross sections. If the composition can be estimated or measured, for example 



Figure 6.3. Notional sketch of spectrally dispersed 
Rayleigh scattered light (solid line) which has been 
Doppler shifted from the incident laser light (dashed). 
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by Raman scattering, then the effective cross section for the gas mixture can be determined, 
allowing the density to be determined from the measured Rayleigh scattering. Under conditions 
where the pressure and composition are known or can be determined and where the perfect gas 
law applies, the gas temperature can be inferred from the measured density. 5 Such measurements 
can suffer from interference from Mie scattering from particles and laser scattered light. 5 To 
avoid Mie interferences, the gases can be filtered to remove particles. Careful experimental 
design and blackening of surfaces can minimize scattered laser light. 


Spectral analysis of Rayleigh scattering yields additional parameters that can be measured. 
Figure 6.3 shows a schematic of the dispersed Rayleigh scattered light (solid line), separated 
from the incident laser light (dashed line). The spectral shift of the Rayleigh from the laser is 
caused by the Doppler shift of the gas relative to the incident light beam and depends on the 
detection angle. The width of the Rayleigh spectrum depends, in part, on the gas temperature. 
The amplitude of the scattering, as described above, depends on the gas density (and 
composition). Thus, in principle, by spectrally dispersing the Rayleigh scattered light, the 
temperature, velocity, and density of a gas can be measured simultaneously. These effects are 
detailed below with examples given. 


Revisiting Eq. (6.5), notice that for Raman scattering to occur, (j^j must be nonzero. This 

happens when a molecule exhibits a change in its polarizability with vibrational displacement. 
For example, as the atoms in the N 2 molecule vibrate, they become less polarizable as the atoms 
approach each other (charges exhibiting more force because they are closer to each other) and 
more polarizable when further apart (charges more weakly interacting). This results in a nonzero 


(— ) . Such vibrations are Raman active. On the other hand, some vibrational modes do not 

Va<2/ 0 

change the molecule’s polarizability, for example the asymmetric stretch of CO 2 has the 


following vibration pattern: ()-> <-C 0~>. Such modes have = 0 and are termed Raman 

inactive. Further quantum mechanical selection rules and molecular structure considerations 
(atomic mass, bond length, moment of inertia, etc.) determine the shape and structure of the 
Raman spectra. 5 Since the frequency of Raman spectra depend strongly on the individual 
species’ molecular structure (each vibrational resonance occurring at a different energy), Raman 
spectra from different molecules appear spectrally separated when they are dispersed, typically 
using a grating-based spectrometer equipped with a camera to acquire the spectra. Because 
Raman scattering is species-specific, it can be used to measure individual species concentrations, 
where in Eq. (6.6), N is the species being detected and the differential cross section is that of the 
individual species. 


6.3 Advantages and Disadvantages of Rayleigh and Raman Scattering 

Rayleigh and Raman scattering have many inherent advantages compared to other 
measurement techniques. Both use just a single laser and the laser need not be resonant with any 
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particular molecular resonance as in absorption or laser induced fluorescence (LIF). This allows 
high-powered, fixed frequency lasers to be used for Rayleigh or Raman. Since they are linear 
laser techniques involving only a single excitation beam, Raman and Rayleigh are relatively easy 
to set up and understand, and the data are, in general, easier to analyze than nonlinear techniques 
like coherent anti-Stokes Raman spectroscopy (CARS). Neither Raman nor Rayleigh is sensitive 
to collisional quenching, a phenomenon which complicates the interpretation of LIF signals. 
Absolute intensity calibration of both techniques is relatively straight forward and easily 
performed. Raman and Rayleigh can be performed simultaneously using the same laser to 
measure species concentrations (Raman), temperature (Raman and/or Rayleigh), density (Raman 
and/or Rayleigh) and velocity (Rayleigh). 

The major disadvantage of Raman and Rayleigh scattering is the low signal intensity. The 
low signal from Raman scattering generally prevents minor species (less than a few percent by 
mole fraction) from being detected using this technique. Because the signals are so low, large 
(low /-number) collection optics are usually used. Ideally, the detection optics need to be placed 
close to the measurement region and, in ducted flows, large windows are required. Hypersonic 
and combusting flows can sometimes be luminous. Natural luminosity, like spontaneous 
Rayleigh and Raman scattering, emits light in all directions. In order to improve the signal-to- 
noise ratio, temporal, spatial, and spectral filtering can be used to collect the desired radiation 
and block unwanted natural luminosity. For example, gated detection of pulsed signals (or lock- 
in detection of continuous signals) can be used. 

6.4 Translational Temperature, Velocity, and Density Measurements 

The gas temperature and velocity can be determined from the Rayleigh spectrum by 
resolving the Doppler broadening and Doppler shift, respectively. The Rayleigh scattering 
linewidth is typically in the range of 1-6 GHz (0.03-0.17 cm' 1 ) for supersonic and hypersonic 
flow experiments in the range of a few hundred to a few thousand kelvins, depending on the 
temperature and angle of incidence of detection and collection of the light (see Ref. 229 for 
detailed information about the Rayleigh scattering lineshape and angular dependence). A high- 
resolution laser can be used so that the broadening caused by the laser’s lineshape is negligible 
compared to this Doppler broadening. For example, injection seeded, pulsed Nd:YAG lasers 
typically have linewidths of about 0.1 GHz (0.004 cm' 1 ) while continuous sources can have 
linewidths that are orders of magnitude smaller. To spectrally resolve Rayleigh scattering, a 
high-spectral-resolution instrument is required. Two methods, described below and also in the 
introduction, are typically used: gas vapor cells and Fabry-Perot etalons. 

A gaseous cell filled with iodine vapor can be used in combination with excitation in the 
visible wavelength range, where L has many absorption lines, by placing a low-pressure gas cell 
containing crystalline iodine and L vapor in front of a detector or camera. The diameter of the 
cell should be larger than the diameter of the collection lens. The cell length, pressure and 
temperature of the I 2 cell are chosen to control the L gas concentration and spectral line shape 
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and therefore the absorption magnitude and profile. Attenuation of transmitted light by a factor 
of 10 5 can be achieved using high-resolution cw lasers, though it is more difficult to achieve 
attenuations greater than -10 3 with off-the-shelf injection-seeded, pulsed Nd:YAG lasers. 236 The 
absorptions can be used to reject spurious scattered laser light, while passing the pressure- and 
Doppler-broadened and Doppler-shifted Rayleigh scattering, allowing background-scatter-free 
density measurements. These sharp absorption features can also be used for temperature and 
velocity measurements by using the steep edge of the absorption spectral profile to provide 
spectral resolution. Since the absorption spectrum of the gas filter is well known, if the laser 
frequency is scanned across the filter, then the Doppler-broadened and Doppler- shifted Rayleigh- 
scattered light will transmit through the filter, being acquired by the camera. Each pixel on the 
camera will then have obtained a convolution of the Rayleigh-scattered light with the absorption 
spectrum. These spectra can then be deconvolved to determine the gas temperature through the 
thermal broadening and gas velocity through the observed Doppler shift. Miles et al. 
demonstrated this approach in a Mach 2 pressure-matched jet flow, shown in Fig. 6.4. 
Variations of this technique can provide both time-averaged and single-shot measurements. For 
single-shot measurements, the dynamic range of velocities that can be measured is limited, 
although by increasing the gas cell buffer gas pressure, the dynamic range has been increased by 
Elliot et al,. who made single-shot measurements with -9% uncertainty over a velocity range 
from 200 to 600 m/s. 237 ’ 238 
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Figure 6.4. Temperature (left), pressure (center) and velocity (right) measured in a Mach 2 supersonic 
jet flow using Rayleigh scattering observed through a gas vapor cell using a frequency-scanned, 
injection-seeded Nd:YAG laser. 229 © IOP Publishing. Reproduced with permission of the authors and of 
IOP Publishing. All rights reserved. 

By using a Fabry-Perot etalon, single-shot Rayleigh spectra can be obtained, allowing 
instantaneous (in -100 ns) and simultaneous measurement of temperature, density and velocity. 
A Fabry-Perot etalon consists of two planar, reflective surfaces that cause interference, 
dispersing the transmitted light spatially, so that in the focal plane of the etalon, the spectrum of 
the light is separated spatially and can be resolved spectrally. Two different strategies have been 
developed to acquire these spectra. The first uses a CCD camera in the focal plane of the etalon 
to capture the dispersed spectrum. " The second method uses spatial masks or mirrors to 
direct the Rayleigh-scattered light to single-point detectors such as photo-multiplier tubes 
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(PMTs). 243 ' 244 Typically, CCD cameras read out more slowly than single-point detectors such as 
PMTs, so using PMTs generally results in higher-speed detection. Measurement rates up to 32 
kHz have been reported using this approach. 244 However, using a CCD camera offers several 
measurement advantages, described below. 

Figure 6.5 shows one such example of CCD-based detection of Rayleigh-scattered light from 
Bivolaru et al. A pulsed Nd:YAG laser is focused into a heated Mach 1.6 jet flow. The 
Rayleigh-scattered light is collected at right angles by a lens system that down-collimates the 
collected light and passes it through a solid etalon. The etalon-processed light is then focused on 
an electron-multiplying charged-coupled device camera (EMCCD) where the interference 
fringes are realized. This image is called an interferogram. An optical fiber (OF) directs some of 
the spectrally narrow laser light into the optical path so that a circular reference fringe will 
appear in the interferogram, as in Fig. 6.6(a). The pair of oval patterns located in the white 
rectangle show Doppler-shifted, Doppler-broadened, Rayleigh-scattered light collected by the 
lenses and dispersed by the etalon. A similar pattern is shown on the left side of the 
interferogram. These four oval patterns originate from four different spatial locations a few mm 
apart in the flow, thus resulting in four simultaneous, spatially-separated measurements. The gas 
velocity, V, and the Doppler shift, A / , are related by: A / = ((k s -k 0 ) ■ V)/A, where A is the 
wavelength of the incident light, k a is the wave vector of the incident light and k s is the wave 
vector of the collected light. (The wave vector points in the direction of light propagation and has 
a magnitude of 2n/X.) In the vector diagram shown in Fig. 6.5, k 0 i is the wave vector of the 
incident light and k s i is the wave vector of the collected light. The observed Doppler shift 
measures the velocity component Vi in the direction defined by k s i - k 0 i, which bisects the angle 
/ 3 between the laser and collection wave vectors. Similarly, a mirror M r reflects this incident 
beam back through the lens, L?, and into the measurement volume with wave vector k 0 2 resulting 
in collected Rayleigh signal having a wave vector k S 2 . This signal is sensitive to the velocity 
component VT The reflected beam was slightly misaligned in the downward direction so that 
the two measurements would be spatially separated on the interferogram, shown in Fig. 6.6(a). 
Furthermore, the geometry of this experiment was constructed so that the Vi would be 
orthogonal to the jet axis, measuring a radial velocity component while V 2 was parallel to the jet 
axis, measuring the axial velocity component. In Fig. 6.6(a), four spatially separated 
measurements of axial velocity are collected in the bottom half of the interferogram while four 
measurements of radial velocity are collected in the top half. Figure 6.6(b) shows the two boxed 
peaks after they have been processed to linearize the interferograms and to bin the data into 
single spectra. 240 ' 241 These spectra were fitted with Gaussians to determine the Doppler shift and 
broadening associated with Rayleigh scattering, relative to the reference peaks. Subsequent 
work by this team and others has been able to simultaneously determine the gas velocity, 
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temperature and density from similar spectra.' ’ 

In Ref. 241, Bivolaru et al. reported velocity measurements with a precision of ~40 m/s in a 
flow with -1200 m/s, or about 3% of the maximum velocity. The dynamic range of the 
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instrument was -3000 m/s. (The dynamic range is mainly determined by the thickness of the 
etalon, which sets the free spectral range, or fringe-to-fringe spectral separation.) Thus, 
expressed as a percentage of the dynamic range of the instrument, the measurement precision is 
-1%. To obtain more precise measurements at lower temperatures, the experiment can be 
designed to use higher spectral resolution. 

The use of the camera-based approach has some advantages over the PMT-based approach. 
First, it can tolerate (and in fact, benefits from) scattering from stationary surfaces in the flow. 
Such scattering, as long as it is not too large, provides a reference frequency to determine the 
Doppler shift. Second, having a laser frequency reference in each interferogram makes the 
system insensitive to variations in the etalon transmission spectrum or the laser wavelength. In 
typical PMT-based experiments, the etalon must be temperature controlled, vibration isolated (if 
it is an air spaced etalon) and the laser frequency must be carefully controlled. Uncontrolled 
drift in either the etalon or laser would result directly in a systematic error in the PMT-based 
approach while it is automatically corrected in the CCD-based approach. A third advantage of 
the CCD-based approach is that it is somewhat more tolerant of scattering from clusters and 
particles in the flow, lessening the need for gases to be filtered. Scattered light from particles 
typically appears in CCD-based interferograms as spatially distinct, circular artifacts. These can 
sometimes be removed by image processing. However, if the scattering is too large either 
spatially or in intensity, it can corrupt the measurement; even in this case, other spatial locations 
in the flow may yield measurements from the same interferogram. As shown in Fig. 6.6, the 
CCD-based detection allows multiple spatial points to be measured simultaneously, allowing 
measurements at adjacent spatial locations to be correlated. Furthermore, with CCD-based 
detection, it is straightforward to measure multiple velocity components with the same 
instrument, also shown in Fig. 6.6. Finally, CCD-based detection has often been combined with 
pulsed-laser excitation, which has two benefits: it makes the measurement instantaneous, 
occurring in -10 ns, and the signal intensity is much higher, allowing measurement in much 
lower density flows, such as atmospheric -pressure flames and low-pressure gas flows. 
Instantaneous measurements at 1/5 of atmospheric density or lower are possible with this 
approach. 241 
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Figure 6.5. Two-component interferometric Rayleigh scattering system from Ref. 241. Mirrors are 
denoted by M, dichroic mirrors by DM and lenses by L. PBS is a polarizing beam splitter while NBF is a 
narrowband filter. IR is an iris and FPC is a focal plane camera used to monitor the beam alignment. In 
this experiment, the dichroic mirrors were required to filter out other laser beams associated with a dual- 
pump CARS measurement that was being performed simultaneously. Reprinted with permission of the 
authors. 


1.2 



Figure 6.6. (a) Fabry-Perot interferogram of Rayleigh-scattered light obtained by laser beams from two 
different directions as well as laser-light, resulting in the circular pattern, (b) the linearized Rayleigh- 
scattered spectrum obtained from the boxed region in (a), showing best fits to the reference and Doppler- 
shifted light. 241 Reprinted with permission of the authors. 

6.5 Rotational and Vibrational Temperature Measurements 

Raman scattering is sensitive to individual molecular rotational and vibrational transitions. 
Since the population of these lines depends on (actually defines) the temperature, it is possible to 
measure rotational and vibrational temperatures from spontaneous Raman spectra. Figure 6.7 
shows a typical Raman spectrum obtained in a hydrocarbon-air flame at high pressure. A 
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pulse- stretcher was used to extend the duration of the 500-mJ, 532-nm pulse by nearly a factor of 
10, thereby lowering the peak laser power to avoid laser-induced breakdown while maintaining 
high pulse energy to yield a sufficiently high signal-to-noise ratio. This spectrum shows 
many temperature and concentration dependent features. The relative heights and shapes of the 
different bands depend on temperature and gas concentration. Modeling, calibration and analysis 
of such Raman spectra can yield rotational and vibrational temperatures as well as 
concentrations. 



Figure 6.7. Five-hundred-pulse-average Raman spectrum in a high-pressure CH 4 air flame. The 
excitation laser wavelength was 532 nm. A 532 nm filter blocks Mie, Rayleigh and spurious laser 
scattering as well as some of the low-rotational-quantum-number rotational Raman lines. A subframe 
burst gating (SBG) technique was also used to subtract background emission from this spectrum.' 48 
Reprinted with permission of the authors and publishers. Combustion science and technology by 
TAYLOR & FRANCIS INC. Reproduced with permission of TAYLOR & FRANCIS INC. in the format 
reuse in a book/textbook via Copyright Clearance Center. 

Several different strategies for temperature measurement based on spontaneous Raman 
scattering exist. Rotational temperatures can be measured from pure-rotational Raman scattering 
either using high-resolution or low-resolution detection, as indicated in Fig. 6.8. 249 The 
advantage of low-resolution detection is that the spectrum can be acquired simultaneously with 
the same instrument used to acquire multi-species spectra like that shown in Fig. 6.7. A 
disadvantage of this low-resolution technique is that many different species have similar pure 
rotational Raman shifts, so they overlap in the same spectral region close to the excitation laser. 
Such interferences from different species can lead to measurement errors. This technique works 
over a wide temperature range, including at room temperature. Alternately, rotational 
temperatures could be determined by resolving the rotational- vibrational Raman scattering, 

OCA OC1 

typically of N 2 , as is often done for CARS. When flows are in rotational-vibrational 
equilibrium, it is more common to measure the temperature using vibrational bands, as described 
below, because they result in higher signal-to-noise ratios and consequently, more precise 
temperature measurements. 


77 



Envelope (low resolution) 


T=500K A 


A 1 

A 


/ 

/ f 

J.,.. 

\ 

\ 

\ 


510 


520 


530 


540 


550 


510 


520 


530 


540 


550 


Vibrational temperatures 
can be measured from the 
relative heights of different 
vibrational Raman bands. A 
commonly used method 
compares the ratio of the 
Stokes to anti-Stokes (S-AS) 
vibrational Raman bands of 
N 2 . 5 252 254 n 2 is often used 
because it is present in high 
concentration in many 
supersonic, hypersonic and 
combusting flows, resulting 
in adequate signal-to-noise 
ratio. Also, N 2 is used 
because it is relatively well 
understood and well resolved 
spectrally from other species. 

Recall from Fig. 6.1 that the 
anti-Stokes scattering 

originates from an excited 

vibrational state whereas the Stokes scattering can originate from the ground vibrational state. 
Thus, the integrated intensity of the anti-Stokes N 2 spectrum at 473 nm in Fig. 6.7 can be 
compared to the Stokes scattering at 607 nm. This ratio, which is a monotonic function of 
temperature, can be used to determine the temperature. Since this method integrates 
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Figure 6.8. Theoretical calculation of the pure rotational Raman 
spectra of N 2 at two different temperatures and two different spectral 
resolutions, adapted from Ref. 249. © IOP Publishing. Reproduced 
with permission of the authors and of IOP Publishing. All rights 
reserved. 


the signal intensity in the two different bands, it improves the signal-to-noise ratio, allowing 
single-shot determination of temperature. Single shot precisions of -20% of the measured 
temperature have been obtained at several points along a hne using this method. ' This 
technique becomes insensitive to temperature below about 700 K because the low population in 
the excited vibrational state causes low signal-to-noise ratio. 5 


A second method of measuring the vibrational temperature is to spectrally resolve the 
different vibrational bands of a molecule such as N 2 , though the individual rotational lines need 
not be resolved. The relative heights of the vibrational levels can be plotted on a Boltzmann plot 
or can be fitted spectrally to determine the vibrational temperature.' For example, Sharma and 
coworkers' ~ used a KrF excimer laser operating at 248 nm to generate spontaneous Raman 
scattering in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The 
EAST facility, normally operated as a shock tube, was fitted with a two-dimensional converging- 
diverging nozzle so that a shock reflection generated a high pressure (100 atm), high temperature 
(5600 K) reservoir of N 2 gas that expanded through the nozzle. The sudden expansion through 
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the nozzle results in vibrational-rotational nonequilibrium, which was monitored with Raman 
scattering. Measurements were performed at different distances downstream to study the 
evolution of the vibrational relaxation. Figure 6.9 shows a sample spectrum generated and 
spectrally fit to determine the rotational and vibrational temperatures. The vibrational 
temperatures were then plotted versus distance downstream in the flow so that different 
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vibrational relaxation models could be tested.” Though not explicitly stated in the paper, the 
temperature measurement precision appears to be about 5%. 



Figure 6.9. An experimental Raman spectrum of N 2 , fitted for rotational and vibrational 
temperature (left), and resulting vibrational temperatures measured as a function of distance 
downstream of the nozzle (right), adapted from Ref. 255 with permission of the authors. 

6.6 Species Concentration Measurement 

Rayleigh scattering has been used in some specialized experiments for determining species 
concentration. For example, the Rayleigh scattering cross section for propane is 13.5 times 
larger than for air, allowing mixing to be quantified in propane/air jets. 256 However, Raman 
scattering is much more commonly used to measure species concentrations. Sandia National 
Laboratories have extensively used Raman scattering for major-species concentration 
measurements (along with Rayleigh and LIF measurements for temperature and minor-species 
concentration measurements, respectively). 155 This system uses a series of four frequency- 
doubled Nd:YAG lasers to produce temporally- stretched 1.8 Joule pulses which are focused to a 
0.22 mm (l/e“) diameter spot size. A separate laser excites LIF of CO. An imaging system 
directs the collected Raman/Rayleigh/LIF scattered light into the optical analysis system shown 
in the bottom-right panel of Fig. 6.10. The system contains two mechanical chopper wheels (one 
“slow” and one “fast”) that are synchronized with the laser to transmit the signals while rejecting 
flow luminosity. It also has a transmission grating and multiple lenses, beam splitters and cooled 
(low noise) CCD cameras. This system images a 6-mm-long segment of the Rayleigh, Raman 
and LIF probe volume, allowing concentrations and temperature to be obtained along a line. The 
spatial resolution for the Raman measurements was -10 pixels / mm along the 6-mm-long probe 
volume, though the actual spatial resolution was slightly worse due to optical distortions from 
the flame. 155 
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A sample of the resulting concentration measurements are shown in Fig. 6.10. Raman has 
been used to measure N 2 , 0 2 , CH 4 , C0 2 , H 2 0 and FI 2 while CO was measured with two-photon 
LIF excited at 230.1 nm. Temperature was measured from the intensity of the Rayleigh- 
scattered light, using the ideal gas law to convert density to temperature as described above. The 
data are graphed versus temperature to show the correlation of species with temperature and to 
allow comparison with an unstrained flame calculation, shown in blue. The individual (red) 
data points correspond to different single pulses of the laser. Thus, these measurements are 
‘single shot’ having been obtained with flow freezing (-100 ns) time resolution. Data such as 
these have been acquired at many locations in a variety of different flames, allowing the 
development of a large database of turbulent flames. The accuracy of the Raman concentration 
measurements varies from 2% to 10% of the measured concentration, while the precision (based 
on one standard deviation) varies from 0.7% to 7.5% depending on the species. The accuracy of 
the temperature measurement, based on Rayleigh scattering signal intensity, was reported to be 
2% with a l-o precision of 0.75%. 






0 08 


0.10 


unstrained flame calculation 


<fr = 0.73 


0 08 


0 06 


fshl 


0 06 


0 04 


0.02 


D C 2 


0.65 


0.06 


0.05 
0.04 
x 8 0.03 
0.02 
0.01 
n nn 


0.020 


0.00 

0.20 


0.15 


X s 0 10 


0.05 


0.05 


0.00 

0.80 


0.00 


0.60 


1200 1500 
T (K) 


T (K) 


0.015 


0.00 

0.20 


0.010 


300 600 900 1200 1500 1800 2100 

T (K) 


Slow 

wheel 


Holographic 

9ra,in9 jjt* 

| 

m 


amera ; 

J Fast 


Fast 

wheel 


Rayleigh 

camera 


Figure 6.10. Scatter plots of mole fraction and equivalence ratio ((j>) data at one location in an 
atmospheric pressure, methane-air flame compared with a calculation, and a schematic of the optical 
components associated with the collection and analysis of light from the Sandia Raman/Rayleigh/LIF 
imaging system (bottom right). Adapted from Ref. 155 with permission of the authors and The 
Combustion Institute. 
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7. Coherent Anti-Stokes Raman Spectroscopy 

7.1 Introduction 


Comprehensive reviews of the theory and application of Coherent anti-Stokes Raman 
spectroscopy, often called coherent anti-Stokes Raman scattering, or CARS, in a gas have been 
given by others. 5 ’ 109 A brief introduction and overview are provided herein. CARS is a non- 
linear optical process in which three laser beams interact with the gas generating a fourth, laser- 
like signal beam. The energy level diagram of this process (each arrow represents a change in 
state due to photon absorption, an up arrow, or emission, a down arrow, and the length of the 
arrow is proportional to energy change or, equivalently, photon frequency) is shown in Fig. 7.1 


The pump beam and the lower 
frequency Stokes beam interact with the 
gas, a pump frequency photon is absorbed 
and a Stokes frequency photon is emitted 
coherently with the Stokes beam, and the 
gas is excited to a higher energy state via a 
state that in most CARS setups is a virtual 
state (although it could also be a real state). 
The difference between the frequencies of 
the upper and lower state, m T 



Signal 


-’pump 


Figure 7.1. Energy level diagram showing the CARS 
process. 

( n stokes^ 's called the Raman shift. A probe beam photon is coherently scattered from this excited 
state, shifted up in frequency by the Raman shift to form part of the signal beam, and the gas 
molecule returns to its original state. Total momentum as well as energy is conserved. 
Therefore, the momentum of the scattered photons equals that of the incident photons, leading to 
the following equation: 


kpump "b kprobe 


— k 


Stokes 


+ k 


signal 


(7.1) 


Since the k’s are the wave vectors (length proportional to photon frequency with direction the 
same as that of the beam), this equation allows the direction of the signal beam to be found from 
the directions of the pump, Stokes, and probe beams. 


A physical interpretation of this process 
is that the interaction of the pump and 
Stokes beams establishes an optical fringe 
pattern in the gas. If the pump and Stokes 
frequencies are the same, as with a related 
technique called degenerate four-wave 
mixing (DFWM), this fringe pattern is 
stationary. In CARS, the frequencies are 
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Figure 7.2. Planar BOXCARS beam geometry. 
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different and the fringe pattern moves across the interaction region, modulating the intensity at 
any particular point at the beat frequency, i.e., the Raman shift frequency. This moving fringe 
pattern excites a polarization response in the gas which acts as a moving grating from which a 
probe beam photon is scattered, in a manner similar to Bragg scattering, to form the signal beam 
photon. Since the grating is moving, the frequency of the signal photon is shifted relative to the 
probe photon by the Raman shift frequency. (An analogous effect occurs in an acousto-optical 
modulator where acoustical waves in a solid material, such as glass, form a moving grating of 
varying index of refraction, and incident laser light is coherently or Bragg scattered from this 
volumetric “grating” with its frequency shifted by the acoustical frequency. ) 

If the pump, probe, and Stokes lasers are all single frequencies, then the signal is also at a 
single frequency, as indicated in Fig. 7.1. If the pump and/or Stokes lasers are broadband lasers, 
while the probe is single frequency, then the signal is broadband also, and contains a spectrum 
that reflects the variation of CARS susceptibility of the molecules in the probe volume as a 
function of Raman shift. In many broadband CARS setups, the probe and pump frequencies are 

c 0<C 

the same (derived from the same laser source). In the dual-pump CARS technique, ' pump 
and probe frequencies are different and these two laser beams have interchangeable roles (each 
beam performs the role of pump in one CARS process and probe in a second CARS process). 
The signal is thus the coherent superposition of the signal from the two processes, generated over 
two different ranges of Raman shift. When the desired Raman shift is small, as in rotational 
CARS (i.e., when the transitions probed are pure rotational transitions, which differ by relatively 
small quanta of energy and are thus relatively closely spaced spectrally), the pump and Stokes 
beams can be derived from the same broadband laser while the probe beam is single frequency 
(dual-broadband CARS). 5 

Since CARS is a non-linear process (see Section 1.2) it requires high irradiance levels. 
Consequently, CARS signal in a gas is typically generated by focusing pump, Stokes, and probe 
beams, and overlapping them at their common focus. If the beams are initially separate and 
crossed at the common focus, the arrangement is called BOXCARS; if they lie in a plane (Fig. 
7.2), it is called planar BOXCARS. In this arrangement the measurement volume is the small 
region of overlap of all three beams. 

7.2 CARS Theory 5 ’ 109 ' 257 

CARS, like Raman and Rayleigh scattering, arises because of time-varying polarization 
induced in the gas in the presence of electromagnetic radiation. CARS, specifically, arises due to 
the third order susceptibility, Xcars> f° r which the induced polarization is described as follows: 

{^signal) ~ ^oX CARSp(S^pump)P^^Stokes)P^^probe) (7.2) 
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The £”s are the complex electrical field amplitudes. This equation can be substituted into the 
wave equation relating the electrical field to the induced polarization and solved by integration 
along the direction of the signal beam. The CARS irradiance is thus: 


1 signal ^pump^Stokes^probe 


\Xcars\ 2 L 2 


(7.3) 


L is the length of the measurement volume, i.e., the length of the region along which the pump, 
Stokes, and probe beams all overlap. Rigorous calculations of the CARS susceptibility require 
quantum mechanical treatments, but classical derivations are simpler to understand. The gas is 
modeled as a simple harmonic oscillator in which the variation of the normal coordinate, Q 
(intermolecular spacing in the case of vibrational states), with time is described by a second 
order linear ordinary differential equation, with a time-dependent forcing function proportional 
to the average over an optical cycle of the square of the electrical field. (Thus, the forcing term 
oscillates at the Raman shift frequency.) The constant of proportionality in the forcing term is 


proportional to , and consequently to the Raman scattering cross section, where a is the 

optical polarizability of the molecule. The CARS susceptibility is found by solving for the time 
dependence of the normal coordinate, combining with Eqs. (6.1) and (6.2) for the polarization 
(from the chapter on Rayleigh and Raman scattering), and comparing with Eq. (7.2): 


XCARS ^ 
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(7.4) 


The susceptibility is proportional to the molecule number density, N, and the differential Raman 
cross section, and reaches a peak where the detuning, Am = 0) v - (oj pump — O) st0 kes ) , 

approaches zero. The damping coefficient, T, determines the line width. A quantum-mechanical 
treatment yields the following equation: 

= 2/ + Xnr Where K, oc NA, (f§) I . (7.5) 


The index j refers to a particular transition and A ; - is the fractional difference in the population of 
the two states of the gas between which a transition is taking place (the two real levels in the 
energy level diagram of Fig. 7.1); Xnr * s the non-resonant susceptibility, a nearly constant term 
for the CARS process in which all the states in Fig. 7.1, other than the ground state, are virtual 
states. The (real part of the) line shape in Eqn. (7.5) is a Lorentz function with full-width at half 
maximum equal to Tj. This line shape is the same as the (complex conjugate of the) line shape 
for the classical solution when the detuning is small relative to the vibrational frequency. The 
line width depends inversely upon the lifetimes of the two states involved in the transition, which 
depend upon the rate of collisions between molecules (collisional broadening). Light emitted by 
molecules in motion is observed at a slightly different frequency by a stationary observer by an 
amount proportional to the velocity component towards the observer. This effect, when averaged 
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over many molecules travelling in different directions, results in Doppler broadening. Since the 
distribution of molecular velocity components about the bulk mean in a gas in equilibrium is 
Gaussian, the line shape becomes a Voigt profile, a convolution between Gaussian and 
Lorentzian functions. At pressures much above atmospheric, additional line narrowing effects 
occur. 

Since the CARS signal is proportional to \Xcars \ 2 > a CARS spectrum reflects the 
populations of the molecular energy states involved in the transitions. Since these populations, in 
equilibrium, are related via the Boltzmann equation to the temperature, CARS can measure 
temperature. Integrated CARS signal is also strongly dependent upon number density and, in 
principle, could be used to measure density. This is not typically done because of experimental 
difficulties in maintaining a consistent geometry of the laser beams at the beam intersection. 
Small uncontrolled motions of the beams due to refraction in inhomogeneous gas fields, 
movements of the optical system, variation in the quality of the laser beams, etc., cause changes 
in signal intensity, and calibrations fail. However, where two gas species are resonant in a 
spectrum, the ratio of the population of one species to another may be found from the shape of 
the spectrum. Where only one species present is resonant, but Xnr is known, the fraction of that 
species may be found from the relative amplitude of the resonant signal to the non-resonant 
“background”; however, Xnr depends on the number density of all species present. If all species 
but one are resonant in the spectrum then the composition can be fully determined by reference 
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to the non-resonant background.” Use of CARS as a diagnostic tool requires comparison of 
experimental and theoretical spectra. 260 Calculation of theoretical spectra is quite complex and 
numerically time-consuming: after computation of theoretical susceptibility, spectra must be 
convolved with laser line shapes and instrument probe function. 261 

Nitrogen is a very useful species for 
measurement of temperature since it is 
usually present in fuel-air combustion or 
hypersonic flows, and is readily probed by 
CARS. The band head of the No Q branch is 
located at a Raman shift of 2330 cm 1 (units 
of inverse wavelength, proportional to 
frequency divided by the speed of light), 
and is readily accessible using available 
lasers. The Q branch occurs as a result of 
transitions between adjacent vibrational 
states (Av=+1, where v is the vibrational 
quantum number) with no change in the 
rotational state (AJ=0, where J is the 
rotational quantum number). Figure 7.3 shows a portion of the No Q-branch spectrum near the 
band head. Transitions are between v=0 and v=land each peak corresponds to a different J value. 



Figure 7.3. CARS N 2 Q-branch susceptibility and 
spectrum near band head at 1 atm., 1400 K. 
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Since rotation of the molecules affects the energy associated with a given vibrational state, the 
energy difference between v=0 and v=l (at given rotational level) is a function of the rotational 
level. Two spectra are shown in this figure - the first is the CARS susceptibility and the second 
is a computed typical broadband CARS spectrum formed by convolution of typical laser line 
shapes and a typical instrument function with the susceptibility. 261 As may be seen, line shapes 
are typically not fully resolved in experiments. 

Figure 7.4 shows the sensitivity of the signal spectrum to (a) temperature and (b) 
concentration. The calculations are for (a) air and (b) either air (79% N 2 ) or 10% N 2 , 21% O 2 , 
69% H 2 O, with a pressure of 1 atmosphere. The CARS signal strength varies strongly with 
temperature, partly through the effect of temperature on density (AO via the N dependence of the 
CARS susceptibility. By using the Sandia CARSFT code, 260 the pressure and temperature 
dependence was determined to be proportional to p 2 j~ 3 - 5 at moderate to high pressures. 262 
However, the shape of the spectrum also changes. At low temperature, only one vibrational band 
is present, associated with v=0 to v=l transitions, and is relatively narrow (because fewer J levels 
are occupied). At the higher temperatures, “hot” bands are present also, associated with v=l to 
v=2 and v=2 to v=3 transitions, and more J levels are populated. As the fraction of N 2 is reduced, 
the amplitude of the band structure becomes smaller relative to the non-resonant background, 
and sensitivity to N 2 is lost for concentrations less than a few percent. This lack of sensitivity 
below a few percent is typical and prevents measurement of minor species, although a number of 
CARS techniques, for example resonance CARS (e.g., Ref. 263), have been developed to 
circumvent this limitation. 




Figure 7.4. Effects of (a) temperature and (b) concentration on the CARS N 2 Q-branch spectrum. 

7.3 Ultrafast CARS 

Ultra-fast techniques have been surveyed by Roy et al. 264 The CARS technique as described 
above is a steady-state process involving the mixing of four beams simultaneously present. For 
given laser energy available in the pump, probe, and Stokes laser beams, the energy in the signal 
beam is proportional to t , where t is the time over which the measurement is made. 
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Nanosecond (ns= 10 9 s) CARS is performed with Q-switch pulsed lasers where the pulse energy 
is distributed over 5 ns to 10 ns, at a pulse rate on the order of 10 Hz. There are many 10’s or 
100’s of thousands of Raman frequency cycles during each laser pulse, so that CARS is 
effectively a steady-state process. However, with the advent of femtosecond (fs= 10 15 s) lasers 
(e.g., titanium-sapphire), very short pulse widths of order 10 14 seconds are possible at high pulse 
repetition rate (1 kHz or more), and the pulse width is less than the period of a Raman cycle. 
Due to the t scaling, strong signal energy is available with modest pulse energy. 

The theory of fs CARS is different from ns CARS since a steady state problem is replaced 
with an impulsively initiated, time-dependent one. 265 First, a polarization grating is established 
in the gas by the interaction of pump and Stokes beams. This grating evolves in time through 
rotational- vibrational relaxation and then, after some delay (typically up to several hundred 
picoseconds in measurement applications), a signal is generated by scattering the probe pulse off 
the grating. By making measurements at different probe delays the relaxation of the grating, 
which carries with it an imprint of the vibrational and rotational state of the probed molecules, 
may be observed. “Chirping” techniques have been developed where the probe beam is 
broadened both temporally and spectrally in such a way that the frequency of the probe varies 
with time. The temporal relaxation is thus mapped into frequency space, obtained in a single 
laser pulse, and may be analyzed with an optical spectrometer. 266 

Picosecond (ps = 10" 12 s) CARS employs mode-locked solid state lasers with pulse lengths 
on the order of 10 ps and has characteristics of both ns and fs CARS. As with fs CARS, 
generation of the grating by the pump and Stokes beams, and scattering of the probe to form the 
signal, are typically separate steps. However, the probe pulse is still relatively long compared 
with the Raman cycle period, and the signal contains an optical spectrum similar to that of ns 
CARS. 267 An advantage of ps CARS over ns CARS is that the ps pulses required to obtain 
adequate signal-to-noise ratio CARS spectra are better suited to pass through commercially 
available optical fibers without damaging the fibers.' Both fs and ps CARS have the potential 
advantages over ns CARS of much higher data (pulse) rates and the absence, when the probe is 
delayed relative to the pump-Stokes laser pulses, of non-resonant background effects. Another 
advantage of fs CARS is simpler theoretical modeling, since few or no collisions occur in the 
time of the measurement (negligible collisional effects). 264 

7.4 Advantages and Disadvantages of CARS for probing Supersonic, Hypersonic and 

Nonequilibrium Flow 

In hypersonic propulsion systems (scramjets), flow velocities are supersonic in inlets and 
nozzles, and may be supersonic or a mixture of subsonic and supersonic in the combustors. In 
combustors, pressures are roughly one atmosphere while temperatures are similar to those of low 
speed combustion. Thermal nonequilibrium (not typically present in low speed combustion) as 
well as chemical nonequilibrium can be present due to the short flow through time (on the order 
of milliseconds). CARS has several advantages in this application. The signal comes as a laser 
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beam, which means that it can be collected through a relatively small aperture and may be 
separated from non-coherent interferences by spatial filtering. This is particularly useful in an 
engine combustor where there may be emission from the gases and thermal radiation from the 
wall, and where, for structural reasons, it may not be possible to incorporate large windows. 
Typical CARS measurements are spatially and temporally resolved, with ~1.5 mm long and 50 
pm diameter measurement volume, and a 10 ns time scale. CARS is able to non-intrusively 
measure local temperature and composition. Since the Q-branch spectrum reflects the rotational- 
vibrational state of the molecule, the populations of these states may be directly determined, 
which is useful when thermal equilibrium (and therefore a single “temperature”) does not exist. 

CARS has been used much less frequently in hypersonic freestreams, where temperatures 
and pressures can be very low. Due to the previously-mentioned p 2 T~ 3 - 5 scaling of the CARS 
signal, the signal to noise ratio may be low. At low temperatures only the vibrational ground 
state and far fewer rotational states of molecules are populated. Vibrational CARS (e.g., of the 
No Q-branch) then depends upon measuring high resolution spectra of the rotational structure of 
the “cold” band near the band head, which depends upon the pressure-dependent collisional 
effects on line shape, as well as temperature. Measurements of both pressure and temperature, at 
pressures down to about 0.1 atm., have been made in an underexpanded supersonic jet using this 
approach. 269 If the v>0 bands are populated at low rotational temperature (vibrational 
nonequilibrium), this can be easily measured. Alternatively, pure rotational CARS (Av=0, 
AJ=±2) is sometimes used at low temperatures since the lines have greater separation than in the 
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Q-branch, but no information on the vibrational state is obtained.” 

In nonequilibrium flows CARS has a significant advantage that it can simultaneously 
monitor multiple individual state populations. Thus non-Boltzmann distributions can be probed, 
as in the examples below. Using vibrational broadband CARS, separate rotational and 
vibrational temperatures can be determined from CARS spectra. 

CARS has some disadvantages compared to other measurement techniques. It typically 
measures at only a single point, and so does not typically provide simultaneous information at 
multiple spatial locations. Instead, the probe volume is typically scanned around the flow 
allowing time-averaged spatial properties to be measured. The data rate is low in ns systems 
(order 10 Hz) so that in pulsed hypersonic facilities only one single measurement is obtained per 
facility run. Also, its -1.5 mm long probe volume can be too large in certain applications where 
it is desired to resolve small length scales such as shock waves, shear layers or turbulent eddies. 
Many of these limitations are being overcome with fs CARS where ID” and 2D” imaging 
have been demonstrated recently with smaller probe volumes and data acquisition rates in the 
kHz regime have been achieved. 264 CARS requires optical access on two sides of the flow, 
which may limit application in some facilities. The experimental setup for CARS, involving two 
or three pulsed laser systems along with large spectrometers ), is relatively complicated and time 
consuming to set up. Furthermore, the nonlinear nature of the theory of CARS complicates the 
interpretation and analysis of the resulting spectra. However, with the combination of accurate 
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and precise temperature and multi-species measurement capabilities, CARS is used in many 
supersonic, hypersonic, and nonequilibrium flow applications. 

7.5 Mole Fraction and Nonequilibrium Temperature Measurement 

Extensive dual-pump CARS measurements have been made in a dual-mode scramjet burning 
hydrogen.” The experiments were conducted in a facility which provided electrically heated 
clean air to the scramjet combustor via a Mach 2 nozzle. The CARS lasers, which were located 
outside the scramjet lab, consisted of an injection-seeded Nd:YAG, frequency doubled to 532 
nm, an in-house broad band dye laser (Stokes laser) centered around 603 nm with FWHM of 10 
nm, and a commercial narrow-band dye laser centered around 550.5 nm. Beams were relayed to 
the experiment via a translation system that could move the measurement volume and through 
special slotted windows in the scramjet. Beams were focused and crossed in the scramjet in a 
planar BOXCARS arrangement (Fig. 7.2) to form the measurement volume. The signal beam 
was transmitted out of the scramjet, recollimated, separated from one of the pump beams, 
relayed to a 1 m spectrometer, then imaged onto a cooled CCD array with 1340x100 pixels. The 
measured spectra had background subtracted and were normalized by a CARS spectrum in 
argon, which has no resonances, to remove the spectrum of the broad -band laser. The resulting 
spectra were fitted to theory, using the Sandia CARSFIT code 260 to generate the theoretical 
spectra and an in-house fitting code 261 to do the fitting, resulting in temperature standard 
deviations of -3%. The fitted parameters included vibrational temperature of No, a single 
rotational temperature for all resonant species, and mole fractions of No, O 2 , and tF. 

Vibrational nonequilibrium was observed in the flow at the exit of the nozzle and in the 
freestream of the scramjet flow, which could be fitted to around 1200 K in vibrational 
temperature of N 2 and 1000 K in vibrational temperature of O 2 , while the facility air stagnation 
temperature was 1200 K and the nozzle exit static temperature was around 680 K.” The low 
rate of relaxation of N 2 vibrational energy relative to rotational (which rapidly comes into 
equilibrium with translational), and the higher but still relatively low rate for O 2 is consistent 
with known rates in the literature. 276 Measurements in which a varying amount of steam was 
added to the air flow are consistent with the known effect of water vapor to promote 
equilibrium.” Figure 7.5 shows the effect of water vapor in the dual-pump CARS spectrum; 
notice the large effect of steam on the hot band of the N 2 spectrum and a lesser effect on the hot 
band of the O 2 . The intensity of these bands is proportional to the population difference between 
the v=l and v=2 vibrational states (roughly proportional to the population of v=l since the 
population of v=2 is small). 
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Figure 7.5. Experimental dual-pump CARS spectra in a Mach 2 flow with varying amounts of steam 
added to the air. The N 2 Q-branch is on the right and the 0 2 on the left. Reprinted from Ref. [275] with 
permission of the authors. 


Figure 7.6 shows some typical (averaged) measured spectra and fits to theory in the scramjet. 
Spectrum (a) was in the freestream air of the flow and is shown with a fitted rotational 
temperature of 923 K and N 2 vibrational temperature of 1 133 K. No and 0 2 Q-branch resonances 
may be seen. Spectrum (b) is in the combustion plume and shows H 2 rotational (S) lines as well 
as the aforementioned Q-branches, with fitted rotational and vibrational temperatures of 1588 K 
and 1766 K, respectively. Figure 7.7 shows contour maps of vibrational temperature in the 
combustor operated in the “scram” mode. The wireframe represents the corners of the flow path; 
a single hydrogen injector is located at the downstream end of the small ramp visible on the top 
surface of the flow path. Flow enters at Mach 2 from the test facility nozzle and is from left to 
right. The development of the combustion may be seen; combustion is initiated on the top of the 
plume of hydrogen, near the fuel injector, and wraps around and engulfs the plume further 
downstream. 
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Figure 7.6. Typical fits of theory to experimental dual-pump CARS spectra in a dual-mode scramjet: (a) 
freestream, (b) combustion plume. Reprinted from Ref. [275] with permission of the authors. 
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Figure 7.7. Contour plots of CARS-measured mean vibrational temperature in a dual-mode scramjet. 
Reprinted from Ref. [273] with permission of the authors. 

Simultaneous with the temperature measurements, species concentrations were also 
determined from CARS data. The concentration information was derived from the relative 
intensity of the different resonant features in the CARS spectra, comparing to the intensity of the 

tcq 

nonresonant background, and considering that the mole fractions must sum to one." Figure 7.8 
shows CARS mole fraction measurements obtained at the same conditions as Fig. 7.7. The mole 
fraction is uniform air in the first plane at the left of each of the figures. Evidence of the cold FE 
fuel jet is seen in the temperature map (Fig. 7.7) and all three mole fraction maps in the second 
plane, which is located just downstream of the fuel injector (Fig. 7.8). As the fuel jet spreads 
spatially and is consumed by combustion, the presence of No in the center of the downstream 
fuel plume shows evidence of fuel-air mixing. However, the CF mole fraction does not track 
exactly with the No because it is reacting with the Ho. Not until the final plane at the right side of 
each figure, when all the Ho is consumed, does the Oo penetrate to the center of the duct, in the 
wake of the fuel plume. 
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Figure 7.8. Contour plots of CARS-measured mean mole fractions of N 2 (top), 0 2 (middle) and H 2 

273 

(bottom) in a dual-mode scramjet at the same conditions as Fig. 7.7. Reprinted from Ref. [273] with 
permission of the authors. 

7.6 Measurement of non-thermalized population distributions using CARS 

Nonequilibrium flows sometimes result in non-Boltzmann distributions, for example among 
rotational or vibrational states. Montello et al. used picosecond CARS to probe non-thermal 
vibrational distributions in No generated in low pressure (100 Torr) N 2 and air pulsed discharges. 
They used a 150 ps duration, frequency-doubled Nd:YAG laser to generate the two green pump 
beams. The Nd:YAG laser also pumped a red broadband dye laser for the Stokes beam. These 
frequencies were chosen to probe the Q-branch Raman spectrum of No. The CARS signal was 
dispersed on a spectrometer, allowing instantaneous quantification of the N 2 vibrational 
distribution function (VDF). By varying the time between the CARS acquisition and the pulsing 
of the plasma discharge, the time evolution of the VDF could be studied. 


91 







Figure 7.9(a) shows vibrational No CARS spectra obtained at different times after the 
beginning of the -100 ns duration pulsed discharge. The narrowness of the bands indicates a 
relatively low rotational temperature, compared to the VDF. The populations of the excited 
vibrational levels are observed to increase with time (though the rotational distribution stays 
relatively constant). Vibrational levels from 0-9 were observed to have measurable population. 
When the vibrational populations are graphed with respect to vibrational quantum number (Fig. 
7.9(b)) a thermalized population distribution would appear as a nearly straight line on this semi- 
log plot. Instead, departure from the Boltzmann distribution is observed, as evident from the 
non-linear trends in the data. A theoretical model captures much of the nonequilibrium physics 
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occurring in this plasma." 
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Figure 7.9. (a) N 2 Q-branch CARS spectrum obtained in a pulsed discharge operating with 100 Torr of 
N 2 shown for three different measurement times after the beginning of the -100 ns discharge pulse, (b) 
Experimental and predicted vibrational energy level populations based on (a), showing non-thermal 
vibrational distributions which evolve with time. Reprinted from Ref. [278] with permission of the 
authors. © IOP Publishing. Reproduced by permission of IOP Publishing. All rights reserved. 


7.7 CARS summary 

In summary, CARS is a powerful technique for simultaneously and quantitatively measuring 
multiple flow parameters in gas flows with thermal and/or chemical nonequilibrium. If the gas is 
in thermal equilibrium, CARS can be used to measure temperature; if not it can provide 
information on rotational and vibrational (for rotational-vibrational CARS) state populations, 
allowing separate determinations of rotational and vibrational temperature or rotational and 
vibrational distribution functions in non-thermalized flows. It also can be used to quantify the 
relative concentration of the species probed. It is spatially and temporally resolved, with a short 
measurement time (sub-ps to 10 ns) and a measurement volume that is typically on the order of 1 
mm long. Nanosecond CARS is limited by available lasers to measurement rates of about 10 Hz, 
but femtosecond CARS data rates are often 1 kHz or more. Nanosecond CARS is a well- 
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established technique, picosecond and femtosecond CARS in gases is quite recent, due to the 
advent of suitable lasers, and techniques are still under development. The primary limitations of 
state-of-the art CARS systems include its complexity (which has meant relatively few 
applications to large-scale test facilities), the fact that measurements are pointwise and not planar 
or volumetric, and (for nanosecond CARS) the low data rate. Other issues include sensitivity to 
optical misalignment due to vibrations or refractive index variations along the beam path, and 
lack of spectral and collision-rate modeling for some species. 

8. Other Molecular-Based Measurement Techniques 

8.1 Introduction 

Many different surface-based, probe-based and particle -based instruments have been used to 
study nonequilibrium hypersonic flows. However, this manuscript focuses only on molecular- 
based measurement techniques. Optical emission spectroscopy, absorption spectroscopy, Raman 
scattering, Rayleigh scattering, laser-induced fluorescence (LIF) and coherent anti-Stokes Raman 
spectroscopy (CARS) are some of the most commonly used techniques, and for this reason have 
been detailed in the preceding sections. This section briefly describes additional molecular- 
based measurement techniques that have been used, or that show some promise for use, in 
studying nonequilibrium hypersonic flows. 

8.2 Interferometry 

Interferometry uses interference of light to measure small changes in optical path length. 
Variation in gas flow properties along the optical path cause the index of refraction to vary, 
which generates these path length changes relative to a reference path. The two paths are then 
aligned onto a one- or two-dimensional detector where interference (or a so-called 
interfero gram) is observed. This method is very sensitive because it can measure optical path 
differences that are fractions of a wavelength of the light from which they are generated. 

In fluid mechanical applications, interferometry is usually used for flow visualization or to 
measure the gas density. The gas density is related to the index of refraction through the 
Gladstone-Dale equation: 


n-l = Kp (8.1) 

where n is the index of refraction, p is the density and K is the Gladstone-Dale constant, which 
depends on the gas and the wavelength of light. The Springer Handbook of Experimental Fluid 
Mechanics recently reviewed applications of various types of interferometry to fluid 
mechanics.” These include Mach Zehnder interferometry, holographic interferometry, shearing 
interferometry and others. The different interferometric methods can be sensitive to the density 
or the gradient in density. Interferograms typically have light and dark fringes, which 
correspond to interference in the light caused by index of refraction variations introduced to one 


93 



of the paths. Sometimes, to enable post-processing, the optics are purposely misaligned 
producing a series of mostly-parallel sinusoidal fringes that act as a spatial carrier frequency 
allowing Fourier transforms to be used in post processing. Such interferograms are called finite 
fringe interferograms (as opposed to infinite fringe interferograms wherein the two paths are 
precisely aligned). 

As a path-averaged measurement technique, most interferometry techniques require 
experimental methods and/or analysis to determine spatially-resolved information. If 
measurements along multiple lines of sight can be obtained, tomographic reconstruction can be 
used to deconvolve spatial information. Or, in axisymmetric flows, an Abel inversion can be 
used, as described previously. However, in the general case of 3D turbulent flows, these 
approaches are challenged by the short spatial and temporal scales of the experiments, which 
require many simultaneous high-resolution views from different angles. While time-averaged 
measurements in such turbulent flow fields are possible, it is more difficult to use these methods 
to resolve small temporal or spatial scales. 

Interferometry is particularly useful in luminous high-enthalpy flow facilities, where spectral, 
spatial and temporal discrimination can be used to obtain the interferograms while rejecting flow 
luminosity. Interferometry was used extensively to study nonequilibrium hypersonic flows 
several decades ago. See, for example, the use of Mach-Zehnder interferometry to study the 
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effect of nonequilibrium flow on the shock stand-off distance and post-shock density field.” 
More recently, holographic interferometry has been used in expansion tube facilities which are 
simply too luminous to use other techniques such as laser-induced-fluorescence. In one 
experiment, the 532 nm and 355 nm outputs from an Nd:YAG laser were used to produce two 
simultaneous finite fringe interferograms in air passing over a cylinder at 12 km/sec in the XI 
superorbital expansion tube at the University of Queensland, Australia. Both holographic finite 
fringe interferograms were acquired on the same photographic plate. After reconstruction, a 
Fourier transform method was used to process the interferograms. In this experiment the fringe 
shift is sensitive to both the density and the electron concentration, necessitating the two 
measurements to determine the two properties. Density and electron concentration maps were 
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determined from the interferograms.” 

8.3 Background Oriented Schlieren (BOS) 

A relatively new measurement technique that is sensitive to density (through index of 
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refraction gradients) is known as Background Oriented Schlieren (BOS).” ’ ” BOS uses a CCD 
camera equipped with a lens to image a flow placed in front of a background consisting of a 
random dot pattern (or other” ’” random pattern). A reference image of the background pattern 
is acquired with the flow absent and then a second image is acquired with the flow present. 
Refractive index gradients cause the dots to move slightly from one image to the other. The pair 
of dot pattern images is processed using particle-image velocimetry (PIV) software to recover 
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the path-integrated light deflections (pixel shifts) caused by the refractive index (density) 
gradients. 

To obtain the highest possible spatial resolution, both the dot pattern and the flow 
disturbance should be in crisp focus. This suggests the use of small lens apertures to obtain a 
large depth of field. To obtain maximum sensitivity, the flow disturbance should be located as 
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close to the camera lens as possible. Since it is difficult to achieve a long enough depth of 
field to optimize both of these requirements, in practice, the object is usually put closer to the 
background dot pattern, for example from 0.5 to 0.875 of the way from the camera to the 
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background pattern. ’ Also affecting spatial resolution are the camera pixel size, the 
resolution and quality of the background pattern used and the interrogation window size used in 
image processing. 

In a two-dimensional flow, the pixel shifts can be used to directly compute the density 
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gradient and thus the density field.” In an axisymmetric flow, the density gradient varies along 
the optical path and additional processing of the data is required to make quantitative 
measurements. In such axisymmetric cases, an Abel or Radon inversion can be used to compute 
the density field, for example in the supersonic flow over a cone at zero degrees angle of 
attack.” In non-axisymmetric flows, tomography using multiple camera views is required to 
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obtain density information.” The multiple views can either be obtained using a single camera 
that observes the flow from different angles at different times, resulting in time-averaged 
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measurements,” or by using multiple cameras to take many simultaneous images, providing 
time resolution, which is a current research topic. Since BOS measures gradients and is thus a 
relative measurement technique, the density in one point in the flow must be known or estimated 
to obtain the density field. BOS has the advantages of being a low-cost technique and being able 
to provide both qualitative flow visualization and quantitative measurements for large fields of 
view. One disadvantage is that BOS does not produce real-time visualization and measurement 
results since the raw images require processing. Also, as a tomographic technique using a 
limited number of views, it has much lower spatial resolution than laser-sheet based techniques 
such as planar laser induced fluorescence, at least in unsteady flows where it is impractical to 
have a large number of views. 

8.4 Electron-Beam Fluorescence 

In low-density (< 10 16 molecules/cm 3 ) hypersonic flows, electron beam fluorescence is an 
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attractive measurement approach. An electron beam can be directed into the flow to excite a 
line, or swept through the flow to stimulate fluorescence in a plane. Multiple species, such as 
NO and N 2 , can be excited and detected at different wavelengths. The fluorescence signal, if 
quenching is negligible, can be related to species concentration. Decades ago, density and 
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temperature were measured in hypersonic boundary layers using electron beam fluorescence.” 
More recently, electron-beam fluorescence has been used to distinguish rotational and vibrational 
temperatures by spectrally resolving the emitted fluorescence, as well as to measure the gas 
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velocity from time-of-flight measurements, using a technique similar to that shown in Fig. 
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5.10.“ However, at higher pressures, several problems complicate the use of electron-beam 
fluorescence. First, quenching becomes significant, complicating interpretation of concentration 
data and minimizing fluorescence lifetimes used to obtain flow tagging velocity. Second, the 
beam begins to disperse at higher densities, preventing measurement in a line or plane. Finally, 
excitation from scattered electrons becomes significant at higher pressures, which complicates 
interpretation of the images. Nonetheless, the electron beam fluorescence method has been 
shown to be useful in a wide variety of low-pressure hypersonic flows and even has potential for 
measurements in flight. 290 288 

8.5 Degenerate Four-Wave Mixing (DFWM) 

Like CARS, Degenerate Four-Wave Mixing (DFWM) is a nonlinear, four-wave mixing 
measurement technique that involves the crossing of three laser beams to generate a fourth beam 
that is detected. DFWM’s three input laser beams are all of the same color (degenerate) and are 
usually tuned to an absorption in the atom or molecule being probed. The signal beam from 
DFWM measurements is typically analyzed in the spectral domain using scanned or broadband 
lasers to measure a spectrum of signal intensity versus wavelength. DFWM has been used to 
measure gas concentration, temperature and velocity. 5 Compared to LIF, it has an advantage for 
measuring species concentrations because it is relatively insensitive to collisional quenching, 
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owing to the fact that it is an absorption-based technique (not depending on fluorescence). “ 
Since each beam is resonant, it is also relatively sensitive and can be used for detecting minor 
species. 5 Since the DFWM signal is a laser-like beam, it is suitable for studying luminous flows 
since background radiation can be spatially filtered with an aperture. DFWM is usually 
implemented as a single-point measurement, though it has been extended to line and 2D imaging 
applications. 5 The main drawbacks of DFWM are that it is more difficult to align and implement 
than LIF, it has much lower spatial resolution than LIF (because of its elongated probe volume), 
it is sensitive to absorption of the laser beams, it is sensitive to saturation of the molecular 
transitions, and it requires very high beam quality to obtain 2D images. 5 DFWM temperature 
measurements in a flame have reported accuracies and precisions more than a factor of two 
worse than CARS' ’ while also having a larger probe volume. 

8.6 Laser-Induced Thermal Acoustics (LITA) 

Another four-wave mixing technique, laser-induced thermal acoustics (LITA), uses two 
pump beams to generate acoustic waves that scatter a probe laser into a signal beam that is 
detected.' The two pumps and probe can be the same wavelengths, in which case it is also 
DFWM. Or, the probe beam can be a different wavelength than the pumps, which reduces light 
interferences from the incident, high-powered pump beams. Whereas the signal beam from the 
DFWM techniques described in the previous section is usually analyzed in the spectral domain, 
LITA usually uses the time domain (signal beam intensity versus time) to measure gas velocity, 
speed of sound (which can be converted to temperature if the gas composition is known) 293 294 
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and gas pressure.” In the time domain, the signal intensity is characterized by a damped 
oscillation, which can have beat frequencies containing the velocity and temperature 
information. Velocity and temperature accuracies and precisions on the order of 1 m/s and 1 K, 
respectively, have been reported in air flows. 294,296 These accuracies and precisions are the best 
reported in this manuscript for these parameters. LITA is usually implemented as a single-point, 
single-velocity-component measurement. The main limitations of LITA are its relative 
complexity and that the measurements are obtained over relatively long probe volumes - on the 
order of 1 to 2 cm. While this is much longer than most turbulent length scales occurring in 
hypersonic flows, it can be sufficient to probe laminar flows if the probe volume is oriented with 
its long axis oriented perpendicular to the steepest gradients in the flow. 

9. Conclusions 

This manuscript has described numerous molecular-based measurement techniques that are 
potentially applicable to the investigation of hypersonic nonequilibrium flows. One might ask, 
“Which technique is best?” or “Which one should I use?” The answer to the first question is 
that no one is best. The different techniques measure different properties in different ways and 
have relative merits. In addition, flow conditions and facility constraints may prevent 
application of specific measurement techniques. The answer to the second, more important 
question, depends on many factors. It depends on the measurement requirements as well as the 
past experience of the research team. It depends on the available equipment, the budget for new 
equipment and the time frame available to make the measurement. Some techniques, such as 
CARS, require a year or years of training to gain enough experience to apply the technique 
expertly. Dual-pump CARS requires home -built equipment (i.e. a broadband dye laser) as well 
as two commercially available lasers, a spectrometer, a low noise CCD camera, and other 
equipment. It is therefore relatively time consuming and expensive to set up. If a researcher or 
research team has experience with similar technology (for example Nd:YAG lasers, dye lasers, 
CCD cameras), then the work will progress faster. Less complicated (easier to set up and to 
understand) methods might yield usable results in a shorter time frame. These practical issues are 
certainly important, but a major consideration in planning an experiment is tailoring the chosen 
technique to the measurement requirements, the characteristics of the facility and the expertise of 
the scientist or engineer performing the experiment. 

A measurement campaign should begin by answering the following questions, for example 
by interviewing the final “customer” of the data: 

• What parameter(s) need to be measured? 

• Must multiple parameters be obtained simultaneously to determine correlations? 

• What spatial resolution is required? 

• Is imaging required or are single -point or line measurements sufficient? 

• What temporal resolution is required? (e.g. time required for a single measurement) 
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• Do measurements need to be time resolved? (e.g. a continuous sequence of data) 

• What accuracy is needed? 

• What precision is needed? 

• What quantity of data is required? 

• When is the data needed? Is instant (real-time) data required? 

• Where in the flow are measurements required? (inflow, exit, near walls, etc.) 

• What type of optical access is available? 

• Can (toxic) seed gases be introduced? Will they influence the properties being measured? 

• What is the ordered priority of the above requirements? 

This list of requirements must be matched with available measurement technologies such as 
those reported in the sections above. There is rarely a perfect solution for this exercise. Instead, 
there is a compromise between the needs of the customer and the capabilities of the 
instrumentation team and the chosen technique. Goals should then be established before the 
research actually starts. This exercise is fruitful for the measurement scientist because it often 
results in a need for a new measurement technique that can be developed if existing techniques 
are not suitable 
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